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FOREWORD 

This document is the final report of a research 
program performed by Lockheed Missiles Space 
Company, Inc., Huntsville Research & Engineering 
Center, while under contract to NASA-Marshall Space 
Flight Center, Contract NAS8-25712. The report 
summarizes a research effort accomplished between 
29 June 1971 and 15 March 1972. Tasks completed 
before this period were reported in an interim report, 
LMSC-HREC D225157, dated June 1971. The NASA 
technical monitor of the contract was Mr. Frank Bugg, 
S&E-AERO-DDS. 
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SUMMARY 

The marker-and-cell numerical technique was applied to the study o I 
axisymmetric and two-dimensional flow of liquid in containers under low 
gravity conditions. The purpose of the study was to provide the capability 
for numerically simulating liquid propellant motion in partially filled con- 
tainers during a docking maneuver in orbit. A computer program to provide 
this capability for axisymmetric and two-dimensional flow was completed 
and computations were made for a number of hypothetical flow conditions. 

In order to extend the numerical simulation capability to more realistic 
flow conditions, a research effort was undertaken to develop a three- 
dimensional marker-and-cell computational technique. For this initial effort 
container boundaries were limited to rectangular shapes. A pilot computer 
program was successfully developed as a result of this research effort. The 
computer program requires 64K core storage with four drum areas for tem- 
porary storage. Computations were made for several test cases with reason- 
able results obtained. This pilot program can be more fully developed to 
include such features as the capability for treating containers with curved 
boundaries. 
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Section 1 
INTRODUCTION 

The Space Shuttle and Space Station programs as currently planned 
have generated much interest in recent years in the dynamic behavior of 
large amounts of liquids stored in containers under low -gravity conditions. 

Of particular concern is the motion of the liquids, and the resulting forces 
that are exerted on the container walls during a docking maneuver in orbit. 

Until recently there have been no satisfactory analytical means for 
calculating the transient flow characteristics for liquids in containers. To 
develop such a means a research effort was undertaken to adapt the Marker- 
and-Cell (MAC) numerical technique (Ref. 1) to the axisymmetric and two- 
dimensional flow of liquids in containers. The capability for treating curved 
container boundaries was included in this adaptation. The resulting computer 
program will satisfactorily predict flow properties, including forces and 
moments on container walls. Results of this research effort are described 
in Ref. 2. 

The success of the axisymmetric and two-dimensional liquid dynamics 
computations described in the preceding paragraph encouraged further effort 
to develop a three-dimensional liquid dynamics numerical technique. The 
research effort described in Ref. 2 was therefore extended to include this 
additional effort. Results of the extended research effort are described in 
this document. 
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Section 2 
FORMULATION 

2 . 1 GOVERNING EQUATIONS 

The differential equations which govern the transient flow of a viscous 
incompressible fluid are 


J 


= - (v -V) v - + P F 2 v + g 


(1) 


D 2 V ♦ v = 0 


( 2 ) 


where 


v 

<P 

P 

P 

t> 


velocity vector 
£ , 


3 n = P r © 8 8ure function 
H 


pressure 
mass density 

kinematic viscosity coefficient 


g = equivalent gravitational acceleration vector, 


»* 


In Cartesian coordinates, Eqs. (l)and (2) may be written as 


9u 9u _ 9uv . 9uw 

9t “ * 9x 9y 5z 


9v 

9t 


9w 
9 t 


9uv _ 9 v 
9x " 9y 


9vw 

9z 


9uw _ 9vw _ 9w‘ 
9x 9y 9z 


1 W ■ - / 9 2 w 9 2 w + 9 2 w\ 

"W 17 8. 2 / 


( 3 ) 


f ’ 1 
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D 0 §u , 9v + Bw a 0 
u n 8x + 9y 9 * 


whoro u, v,w and g x , g y , g g arc the components of vectors v and g 

in the x, y, z-directions, respectively. The stress tensor of the system is 


a =-(/?+ Zu -|— 

xx r 9x 

A v 

a =-</?+ Zv r~ 

yy 9y 




_ ,,/9u , 9 v \ 

V \»y 9*/ 

_ .. /9v , 9w\ 

yz~ \9z 9y / 

_ _ /9w , 9u\ 

^zx ” ^ \9x 9z/ 


In this study, the fluid is considered to be at rest initially; i.e., v(x, y, z, o) 
p (x, y, z, o) = 0. The boundary conditions of the fluid are 


v n = a t = 0 


cr = cr = 0 
n t 


(at a rigid boundary) 
(at a free surface) 


where a variable with a subscript n or t denotes the normal or tangential 
component of the variable, respectively. 
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2.2 FINITE DIFFERENCE FORMULATION 

Cubic meshes are employed in writing the finite difference equations 
of the formulated problem. The velocity components and the pressure of a 
fluid are specified at 'ho boundaries and at the center of a coil, respectively 
As shown in Fig. I, Eqs. (3) and (4) of cell (i, j,k) may bo expressed in the 
following form: 


1 , n+1 ,, . 

Tt 'Vij.k’^+ij.k* 


* - Tt [< u i + l,j,k > 2 - ‘“i.j./] 'Ty [ (uv) i+|,j+ik ' (u V) i+i.j-l.k 


k [ (uw) i+i,j,k+i- <uw >i+i. j.k-i] • “’i.j.k 1 


+ v [-^i <u i+f.j.k tu i-i,j,k- 2u i+i,j.k> 


+ -Z2 (u i+i,j+l.k + u i+i.j-l.k ' 2u i+ij.k> 

6y 


1 


+ 6z Z (Ui+ 2»J* k+l + “i+i j* k_1 ’ 2U i+Uk> I + g x 


JL / n+1 . 

fit V i, j+i', k " i, j+i,k' 


k [ (uv) i+i,j +i.k • (uv) i-i,j + i,k] • ir [ <v i,j+i.k> 2 • (v i. j .k ) ’ 


~k [ (vw >i.j+i.k+i • (v ' w) i,j+i,k-i] • -^‘n.j+i.k- ^t.j,k> 
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♦ " |’t z<v i+>.j+i.k + Vi -i,j+i,k- 

+ ^7< v l,J+i,k + v l.j-i,k ' 2v i,j+i,k ) 

+ “2 (v i.j+i,k+l + v i,j+i,k-l ' 2v i,j+i,k ) ] + 8 y 


tt (w i n j‘k+r w i,j.k+i ) 


- to [ <uw) i+i. j.k+i - (uw) i-'i,j.k+i] ■ 7y [ (vw) i.j+i, 

■ (vw) i,j-i.k+i] - to [ <w i,j,k+i> 2 - lw i.j.k' 2 ] - to <n,j.k+i 


“ <Pi 4 J + V 7. 


i» j*k 


1 + w 


i - 2w, 


2 ' w i+l,j,k+i T w i-l,j,k+£ “”i,j,k+i 


+ (w i,j+l.k + i + w i,j-l,k+|- 2w l.j.k + i» 
+ -h (w i.j,k + -} + w i.j.k-i • 2w i,j,k + i>] + 


D i,j,k 3 to <u i+i,j.k ' “i-i.j.k* + 6y (v i,j+{,k ' v i,j-i,k ) 


+ -r— (w 4 1 4,1 - W. . , x) 
6z i,j#k+a 
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whoro 


<u i,J.k> n VS,J,kM,J.k 


* v i,j.k* " v i,j+i.k v i.j-a.k 


(w i,j.k> “ w i.j.k+i W i,j,k-J 


(uv ’i+'i. j+a.k *3 <U l+a.j+l.k + “i+i. j.k* <v i+l, j+a.k + v l,j+a. k* 
<vw >i. j+i, k+i 4 (v i, j+i, k+1 + v i, j+i, k> (w i,j+l,k+i + w i. j, k+i> 
<WU, i+i.j.k+i 4 < w t+l. j.K+i + W i,j.k+i )(u i+i.j,k+l + u t+i,j,k>- 


Note that the superscript n+1 refers to values at time (n+l)6t. Where there 
is no superscript it is understood that reference is made to values at time 


n 6 1. 


The boundary conditions at a free surface car. be written as 


Zlf , \ 

^i.j.k = fix U i+i» j» k " u i-i,j,k' 


(o a 0) 

' XX 


9 i,j,k“ fiy ^ V i,j+i,k “ v i, j-i,k^ 


(a . =0) 

yy 


Zlf , \ 

^i» j» k 3 fiz W i» j* k+i * W i,j,k -r 


(a = 0) 

' zz 
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t\y ^ u l+4,jtl,k U i+1» J» k^ * 6x^ V i+l ( Jt»»k V i*j+i*k^ 


<"xy ° ®> 


(lOd) 


6z ^ V i* j+a, k+l ’ v l,J+4,k> + 6y ^ w i, j+i,kl*i ’ w i.J»k+^ 0 


(a o 0) 

' ya ' 


(lOe) 


~k (w i+i,j,k+i “ w i, j,k+i> + £ (u i+ij,k+i ■ u i+*,j.k> s 0 


<v B 0) 


(I0f) 
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Section 3 

NUMERICAL SCHEMES 

3.1 THE MAC COMPUTING TECHNIQUE 

Similar to Eqa.(8a), (8b) and (Be), tho u, v and w velocity components 
of colls (i- 1, j, k), 1, k) and (ir 4 »k-t), respectively, can bo found. Sub- 

j^*{' j[ 

stituting these two sets of equations into Eq. (9) and assuming y j,, s t , a 
tentative pressure field (in terms of ( P ) of the fluid at time (n+1)* so V i, 


! 

I 

I 

I 

I 

I 

I 


^^i+l.j.k + n-l.j.k' 2V ’i,j,k )+ { *2 fe ’i,j+l,k' ! ^i. j- l.k ‘ ^i.J.k* 

+ ^K,j i k t .^i,i,k-|- 2 »i, J , k l=Y- Q i,i,k t 6 i ,i,k (U > 


where 


‘i.j.k 


* [ <u l+l,j,k )2 + <u i-l,j,k )2 " 2 (u i.j.k > ] 

+ pr [ (v i.j+i.k> z + (v i.j-i,k > 2 • 2 (v i.j,k> 2 ] 

+ [<»i.j,k +1 > 2 + '*i,j.k-i> 2 - 2 Kj.k> 2 ] 

+ «^r[ (uv) i+i.j+i.k + (uv) i-ij-i.k • (uv) i+i,H.k 
• (uv, i-i.j+i.k] 


** 
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• (uw) i-i.j,k+{] 

+ ife ( (vw) l,j+i.k+i +(vw »i,j-ik-K (vw) i.j+i.k-i 


A 


<H 


I 

I 

1 

! 


I 

I 


• <vw) i,j-ik + i 


and 


D i* j.k ” V 


[■^T (D i+l,j,k + D i-l,j,k ' 2D i.j.k ) + {y 2 <D i,j+l.k 
+ D i,j-l,k ' 2D i,j,k* + 77 (D i,j,k+l + °i,j.k-l 


(12) 


• 2 D i.j.k» 


To find the correct pressure field of the fluid at time (n+1) fit, an iterative 
process is used. Using Eq. (11) the pressure of cell (i» j»k) a-t the comple 
tion of (h+ 1 ) th iteration is expressed as 


(13) 


V 7 


h+1 
i. j»k 


L* w iiM.k + ^i.j.k’ 


■y 1 5 ' ~2 

fix* 6 y fi z 


♦ 77 <V-i.k> + 77 <j.«* + <V.k-.> 

- s i.j.k] ' 


n h 

“^i.j.k , 
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whore a is an over-relaxation parameter and 


‘i,j,k 


j i K „ q 4* D 
fit Q i, jik + D i, jik 


( 15 ) 


The value of a can be taken between 0 and 1, and it is introduced for speeding 
up the iteration process. Hence, the MAC method calculates a set of velocities 
at time (n+1) fit from the velocities and pressure at time nfit. Then, using 
the incompressibility property and the boundary conditions, an iterative pro- 
cess is used to compute the pressure of the fluid at time (n+1) fit. The pro- 
cedure is repeated in the next computing cycle. 

The major steps of a MAC computing cycle are: 

Step 1: Use Eqs. (8a), (8b) and (8c) to calculate a set of velocities 

u n _*i l . , , v n+1 , i . and wf* 1 w ,i , then compute S. . . 

i+2»J»k i, J+2*k i,j*K+a l »J* K 

from Eqs, (9), (12), (13) and (15), 


,h+l 


Step 2; Iterate the pressure field, j until the following 
equation is satisfied throughout the entire flow field 


I* 


h+1 

kiJi 


<p' 


<flk+l + (p ® . . I 
k, j,k ^i,j,k 


i.i. k | 


< c 


<P 


where is a small constant which is chosen to pro- 
vide the necessary accuracy of a solution. 

Step 3: Displace the marker particles according to their local 

velocities. 


Step 4: Adjust the boundary velocities and pressure of the newly 

obtained flow field. 


1 
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3.2 BOUNDARY CONDITIONS 

To facilitate the conatruction of the numerical scheme and direction of 
the execution sequence, the following four types of cells are defined for identi- 
fying the status of a cell in a computing cycle: 

1 . Empty cell (E): A cell having no marker particle, 

2. Boundary cell (B): An empty cell whose boundary face 
forms a portion of the wall of a container, 

3. Surface cell (S): A cell having marker particles and 
neighboring with at least one empty cell, and 

4. Full cell (F): A cell having marker particles and not 
neighboring with any empty cell. 

Since there are many possible arrangements of empty cells around a surface 
cell, a surface cell is further classified to 63 cases (see Fig. 2). 


The velocities of a boundary cell are calculated to satisfy the conditions 
of a rigid smooth wall. The velocities of a surface cell are computed to satisfy 
the incompressibility property and the free surface condition of the fluid. Equa- 
tions for computing the velocities of a surface cell are given below, and the 
appropriate equations to be used in each case are listed in Table 1. 


u i+i» j*k 


jl ) ( 16 ) 

2 


W 


i* j* k+i = w ; 


ty . OX / 

Vij,k-67 (v i.j+ik- V i.H,k 6z W i,j,k+J w i,j,k- 

. . . ‘ lx (u i+i, j,k ‘ U i-i, " 6y (V i. j+ 2 .k V i.j-i*k ) 

jf 2 


V i.j+ik 


U i-£,j,k 


W i i k - 1 
Jl K £ 


=5 V 


i.j-i.k • “l-i.j.k* -& (w l.j.k+i • w i.J.k-i 


(18) 


■ Vi.j.k + l7 ,v i.j + l.k- v i.H.k )+ 6t ,w i,i.kH- w U.K-^ < 19 > 

■i >J .k+t*fe‘"l+Uk ' ‘ ’M-U* <2< ” 


= W 


11 


l(X;KHttO - HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC D225632 


Tabic 1 

EQUATIONS TO BE USED IN EACH CASE 


Case 8 


Equations to be Used 


1, 5, 11, 15,49, 59 

2,7, 10,50,55,58,63 

4, 14,52,62 

8, 13, 56,61 

16,21,26,31,48,53 

32,37,42,47 

19 


(22) , (23) and (24) 

(23) , (24) and (25) 

(22) , (23) and (27) 

(23) , (25) and (27) 
(22), (24) and 26) 

(24) , (25) and (26) 
(22), (26) and (27) 

(25) , (26) and (27) 
(22) and (17) 

(25) and (17) 

(22) and (20) 

(25) and (20) 

(22) and (18) 

(23) and (17) 

(25) and (18) 

(23) and (20) 

(22) and (21) 

(26) and (17) 

(25) and (21) 

(26) and (20) 
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v i. j - ; , k ' v i , j+i , k + li (u i+i, j. k • v i , j. k> + & < w i, j. k+r w i. i. k- 

Ul) 

u i+‘,J.k = VJ.j.k 


(22) 

i • j ^ a i k t » J " a * k 


(23) 

W i. j.k+i = W i,j,k-£ 


(24) 

U i-£» j» k i+2» J>k 


(25) 

V i,j-£,k = V i,j+ik 


(26) 

W i,j.k-£ = W i,j,k+£ 


(27) 


In a three-dimensional MAC formulation, there are also 12 possible configu 
rations involving two empty cells neighboring with a free surface (Fig. 3). The 
velocity component between these two empty cells needs to be considered in 
order to preserve the no-shear stress conditior at a free surface. The following 
equations are used for calculating these velocities: 


Case Equation to be Used 

t \ 

1 u i+£, j> k+1 = U i+i, j,k " fix (w i4l, j,k+£ " W i, j,k+£' 

2 U i+£»j.k-l * U i+£. j.k + fix * W i+l, j,k-£ ’ W i,j,k-£* 

3 u i+£, j+1, k “ U i+£, j,k ■ fix * V i+l, j+£, k “ V i, j+£,k* 

4 u i+i,j-i,k * u i + i,j,k + ^ (v i+i,j-ik- v i.j-i > k> 
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V i-l,j+i,k s V i, j+i, k + 6y (u i-i, j+l,k " u i-i, j.k* 

^ iKI 

v i+l, j+i, k = v i, j+i, k " fiy (u i+i j+i, k ' U i+i, j, k* 


7 

8 

9 

10 

11 

12 


V i, j+i, k+1 
V i, j+i, k-1 

W i-1, j, k+i 
w i+l, j.k+i 
W i» j-1, k+i 
W i, j + 1 * k+ i 


6z ( 

V i, j+i, k " 6y lw i, j+l,k+i 
* V i, j+i, k + 6y ^ w i, j+1 , k- i " 

s W i, j,k+i + |l (u i-i, j,k+l 
~ w i,j,k+i 6z ^ u i+i, j, k+1 

* + (v i,j-i.k+l 

= w i, j, k+| ‘ Tz (v i, j+|, k+1 


W i,j,k+| ) 

W i,j.k-i’ 

• “i-ij.k 1 
U i+1. j. k* 
- v i.j-i.k> 
• v i.H.k» 


The pressure of a boundary cell is considered to be equal to the pres- 
sure of the neighboring full or surface cell. In general, a surface cell does 
not carry a pressure unless it neighbors with one and only one empty cell. 
As shown in Fig. 2, there are six cases of a surface cell which can have a 
pressure. These boundary pressures are defined below: 


^i,j,k ~ fix ^ U i+i, j,k 
- t 

^i, j, k “ fiy ' V i» j+i, k 

«’i, j ,k = i£ lw i.j.k+i 


u. i . ) (for cases 1 and 4) 

j# k 

v. . i , ) (for cases 16 and 32) 
i, j“ a* K 

w. . , i) (for cases 2 and 8) 


t* 


it 
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3.3 DISPLACING OF MARKER PARTICLES 

A volumo-volocity woighting scheme is usod to calculate the velocity 
of a marker particle. Figure 4 shows tho mechanics of the scheme, and 
Table 2 gives the formulas for calculating the velocity components of a 
marker particle in cell (K, J, I) or (i, j,k). At the completion of a computing 
cycle each marker particle will be displaced by a distance v fit. The cell 
status and the boundary velocities and pressure are then readjusted in ac- 
cordance with the new flow field. 

3.4 FORCES AND MOMENTS 

The dynamic loads exerted by a moving liquid on its container are com- 
puted from the following surface integrals 


and 


F 


M 



p n dS 


prxndS 


where F and M are the force and moment vectors, respectively. Notations 
~r and n are the position vector and unit normal of a fluid particle on the 
surface, respectively. In the MAC formulation, the integrals are evaluated 
by summing up the quantities contributed by all cells neighboring with the 
container wall. The loads are calculated at the end of one or several com- 
puting cycles about the X, Y, Z- coordinate system shown in Fig. 9. 




15 
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Section 4 

THE COMPUTER PROGRAMS 


4.1 USER'S GUIDE OF THE LHMAC2 PROGRAM 


I 

I 

I 

I 


f 


A detailed discussion of the Lockheed -Huntsville two-dimensional and 
axisymmetric MAC computer program (LHMAC2) is presented in Ref, 2, The 
program has been modified to include the capability of simulating the transient 
flow of a liquid in a shallow container. Variables to be used in preparing the 
data deck of the LHMAC2 program are defined in Table 3, and the sequence 
and format of the data cards are shown below: 

Data 


Set 

Format 

Variables 

1 

215, 3F8.3, 
215, F8.3 

IBAR, JBAR, DR, DZ, DT, IPHM, PC, ALP 

2 

1 2A6 

NAME 

3 

4F4.1, 8F8.3 

BCB, BCR, BCT, BCL, A, B, C, NU, EPS, 
GR, GZ, VSCALE 

4 

4F10.3, 4110 

T, TWPLT, TWPRT, TWFIN, LPR, NPRT3 

5 

1015, 2F8.3 

TYPE, LI, L2, L3, L4, L5, L6, L7, NXB, 
NYB, UL, UR 

6 

1615 

NSGMTS, NJCI, NJC2, LQUDHT, NPRT2, 
ISUR, ICYCLE, IPLOT, NGLVL 

7 

8F10.3 

(RCOORD(I), 1=1, NSGMTS) 

8 

8F10.3 

(ZCOORD(I), 1=1, NSGMTS) 

9 

8F10.3 

DEPS, VEPS, DBETA, SIGNVN, STH, STR, 
STZ, DS, COFST, RHO, THCKNS 

10 

8F10.3 

(GLVLTT(I), 1= 1, NGLVL2)* 

11 

8F10.3 

(GRT(I), 1=1, NGLVL1)* 

12 

8F10.3 

(GZT(I), 1=1, NGLVL 1) 

13 

1615 

(JHYB(I), 1=1, NJCELL)** 


* NGLVL2 = NGLVL + 1, NGLVL1 = NGLVL + NGLVL 
** NJCELL - NJCI + NJC2 
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Data 

Sot 

Format 


Variables 

14 

1615 

(J,HYR(J), I 

1, N.ICKU,) 

15 

1615 

(NHYU(I), I 

- 1, Nil CELL) 

16 

215, 6FH.3 

NX, NY, XC, YC, XD, YD, UO, VO 


Table 3 

Variables to Be Used in Preparing the Data Deck of the 

LHMAC2 Program 

Description 

0 , 0 , 0 . 0 , 0.0 

0.0 

1 . 0 , 1 . 0 , 1 . 0 , 1.0 
surface tension coefficient 
iteration step 
0.25 

mesh size in the r and z-directions, respectively 
0.4 

time step 

0.0002 

initial and final points of time intervals of the 
equivalent gravitational acceleration 

0.0, 0.0 

equivalent gravitational acceleration in the r and 
z-directions, respectively 

number of interior cells in the r and z-directions, 
respectively (see Fig. 5) 

0 

4 

0, flat free surface; 1, curved initial free 
surface; 2, surface tension effect included 


1 

i r 


1 

I 

I 


Variable(s) 

A, B, C 
ALP 

BCB, BCR, BCT , BCL 

COFST 

DBETA 

DEPS 

DR, DZ 

DS 

DT 

EPS 

GLVLTT(I) 

GR, GZ 
GRT(I), GZT(I) 

IBAR, JBAR 

IPHM 

ICYCLE 

ISUR 
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\ 


fj 


i 

i 


Variables) 

IP LOT 

.mYD(I), UTYB(I), 
NJIYB(I) 

LQunriT 

LPU 

LI, L2, .... L7 

NAME 

NGLVL 

NJC1 , NJC2 
NPRT2, NPRT3 
NU 

NXB, NYB 

NX, NY 

NSGMTS 

PC 

RCOORD(I), ZCOORD(I) 

RHO 
SIGNVN 
STH, STR, STZ 
T, TWFIN 

THCKNS 

TWPLT, TWPRT 
TYPE 

UL, UR, UO, VO 

VEPS 

VSCALE 


Dnucrlpl Ion 

I 

used for defining ITS colln (nee Fig. 5) 

liquid height (hop Fig. B) 

0, no velocity plot; 1, having velocity ploto 

0 , 0 , . . . , 0 

title of a problem 

number of equivalent gravitational acceleration 
intervals 

for defining HS cells (see Fig. 5) 

0, 0 

kinematic viscosity coefficient 

0, 0 

number of marker particles per cell in the r 
and z -directions, respectively 

number of segments 

0, axisymmetric container; 1, channel -type 
container 

for describing container geometry (see Fig. 5) 
Limitation: ZCOORD(I) cannot be specified to 
cross a cell in the z -direction 

mass density 

- 1.0 

for defining a curved free surface (see Fig. 5) 

initial and final time points of a problem, 
respectively 

1.0 

time increments for plots and printout, 
respectively 

0, container with flat ends; 2, container with 
ellipsoidal ends 

0 . 0 , 0 . 0 , 0 . 0 , 0.0 

0.04 

scale factor of velocity plots 
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I 


1 


Var in bin (n) 
(XC, YC), (XD, YD) 


PnHc r Ipt ion 

coordinator! in meters of the lower loft and upper 
right, corner it of a container having fiat endn, 
respect ivoly 


Note that multiple runn can ho made, and execution in terminated by netting 
IDAR o=l 


U The data duck of Sample Problem I in shown below (aloo boo Fig. 6) 


in 

"ip 0,1 «7« 

0 , 107*5 o,oifapf oo 

c 

• 

c 

c 


PROPFLLANT DYNAMICS — CASE ! 

AX 1 SYMMETRIC. H «s 

0.76 M 

1.0 1 

c 

c 

• 

c 

• 

c 

• 

c 

• 

c 

c 

• 

0.0 

t #C -6 

0.000? 

0,6 

0,06P* 

0,P^ 

1 .01 


1 

P 






1 6 

* n 

c 

c 

< 

4 1 

7 


c 

• 

c 

P#0 

6,0 

9,0 

\p.n 

1*5.0 

in.o 

0,0 





0,0 

0*0 

0,37*5 

1 .0 

P.0 

c 

• 

n 

"*P» o 

7?,0 





0,00 

0,04 

1 .0 

-1,0 

4*0 

c 

• 

CVj 

0,00003fl 000,0 

1 .0 




0,0 

0,0077* 

0,71 

1 . 1F6 




- 7,0 

0,0 

c 

• 

c 

0*0 

0.0 

11 *4 

1 1 #4 

11.4 

1 1 .4 

0*0 

0.0 

1 

P 7 

4 « 




1 

to 17 

16 17 




n 

7 3 

P. P 




c; 

* 





n 






-1 







o*o 

0 


1 1 *4 


16 . 7*5 

4.0 

3.0 


lfl.O 

6.0 

0.4 


4.2 USER’S GUIDE OF THE LHMAC3 PROGRAM 


The LHMAC3 program needs approximately, but no more than, 64K 
core space and four drum areas for temporary storage. The program can be 
used to simulate the flow in a rectangular container requiring up to 5000 cells. 
The output of the program will provide the following information: 
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• Plata of 3-D flow and 2 D velocity fioldn of r>. transient flow, 

• Prnnnurn dlntribution on tho container wall , and 

• Dynnmic londn induced by thn moving liquid, 

A brief block diagram which a how a tho organisation of tho LHMAC3 program 
in given In Chart l. 

The aequence und format of preparing the data deck of tho LMMAG3 
program are given below, and variables to be uaed are defined in Tabled. 


Data 


Set 

Format 

Variubleu 

1 

1615 

ITYPE, IBAR, JBAR, KBAR 

2 

1 2A6 

TITLE 

3 

1615 

LNTH1 , LNTH2, LNTH7 

4 

1615 

NMPPUX, NMPPUY, NMPPUZ 

5 

1615 

(IOPT(I), 1= 1,16), (IPLT(I), 1= 1,16), 
(IPRT(I), I = 1,16) 

6 

1615 

NGRT, LHT, NVPLT, (NSEGV(I), 1= 1, 3), 
(J PLANE (I), 1=1,3) 

7 

16F5.1 

(XV(J,I), 1 = 1, NSEGV(J))) Repeat J = NVPLT 

8 

16F5.1 

(ZV(J,I), 1=1, NSEGV(J))) time8 

9 

8F10.4 

(BDRY(I), I - 1,6) 

|L 

10 

8F10.4 

(GRT(I), 1= 1, Jl) 

11 

8F10.4 

(GRX(I), 1=1, J2)* 

12 

8F10.4 

(GRY(I), 1= 1, J2) 

13 

8F10.4 

(GRZ(I), 1=1, J2) 

14 

8F10.4 

DT, DBETA, DX, DY, DZ, EPSA, EPSD, EPSP 
EPSV, RHO, RNU, VSCALE, WALL 

15 

8F10.4 

TIN, TPLT, TPRT, TCOMP, TFIN, TCPU 


*J1 = NCRT +1, J2 = NGRT + NGRT 
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Set up a 

Problem 



Assign Marker Particles 

. 



Compute Velocities 


Adjust Velocities of 
Boundary Cells 


Iterate the Pressure Field 


Displace Marker Particles 


Check Cell Status 


Adjust Pressure and 
Velocities of Boundary Cells 


Calculate Dynamic Loads 
and Generate SC 4020 Plots 


Chart 1 - Program Organisation of the LHMAC3 Program 
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Tabic 4 

Variables to be Used in Preparing the Data Dock of 
the LKMAC3 Program 


Variable(s) 

BDRY(I) 

DBETA 

DT 

DX, DY, DZ 

EPSA, EPSD, EPSP, 
EPSV 

GRT(I) 

GRX(I), GRY(I), GRZ(I) 

IBAR, JBAR, KBAR 

IOPT(I) 

IPLT(I) 

IPRT(I) 

ITYPE 

JPLANE(I) 

I,HT 

LNTH1, LNTH2, LNTH3 
LNTH4, .... LNTH7 
NGRT 

NMPPUX, NMPPUY, 
NMPPUZ 

NSEGV(I) 

NVPLT 

RHO 

RNU 

TCOMP 


Description 

1.0 
0.125 
time step 

mesh size in the x, y and z -directions, 
respectively 

0.5, 0.25, 0.0001, 0.04 

initial and final points of time intervals of 
the equivalent gravitational acceleration 

equivalent gravitational acceleration in the 
x, y and z -directions, respectively 

number of interior cells in the x, y and 
z -directions, respectively 

0 (or blank) 

1 for I = 1, 2, blank for all others 

1 for I = 1, 2, 3, 4, 5, blank for all others 
1 

velocity projection on Y = Jth cell plane 
liquid height* 

equal to IBAR, JBAR and KBAR, respectively 
blank 

number of equivalent gravitational acceleration 
intervals 

number of marker particles per cell in the x, 
y and z-directions, respectively 

4 

number of velocity plots per time point (max. 3) 
mass density 

kinematic viscosity coefficient 

time intervals for computing dynamic loads and 
pressure distribution on container wall 


Measured as number of cells 
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Variables 

TGPU 

TIN, TFIN 

TPLT, TPRT 

TITLE 

VSGALE 

WALL 

XV (I, J), ZV(I,J) 


Description 

number of CPU time allowed for executing this 
run. It is used to save SC -4020 plots in caso 
of max. time being reached. 

initial and final time points of a problom, 
respectively 

time intervals for plots and print -out, 
respectively 

title of a problem 

scale factor of velocity plots 

1.0 

coordinates of container geometry in the x and 
z-directions, respectively (measured as number 
of cells). 


Note that multiple runs can be made, and execution is terminated by 
setting ITYPE a -1. Selection of time step DT is suggested to satisfy the 
following two conditions: 

4 | v | fit < rnin (fix, fiy, 6z) 

and 


4v6t < 


(fix) 2 (fiy) 2 + (fi v) 2 (fiz) 2 + (6z) 2 (fix) 2 
(fix) 2 + (fiy) 2 + (fiz) 2 


Example: The data deck of Sample Problem II, Case 3, is shown below. 





1 16 n 14 

PROPFLLAn r RLORH1NO IN A RFOTAN6ULAR TANK 


15 

R 

14 




9 

9 





0 






1 

1 





1 

1 

1 1 

1 



9 

R 

1 4 

4 4 

6 

cc 

n .o 

16.0 

0*0 




0.0 

0,0 

C 

• 

c 

c 

• 

< 

•— 1 




1 .o 


1 .0 

1 .0 

1 .o 

1 .0 

0.0 


O.F 

C 

• 

IT 



e; .o 


1 .n 

1 .0 

1 .0 


o.n 


o.o 

0.0 

0.0 


o.o 


o.o 

0.0 

0.0 



n. 126 

0.?5 

0.25 

0.26 

< 

c 

• 

c 


1 ooo.o 

O. 000001 

0.25 

1 .0 

0,0 


0.0626 

0.0626 

0.0625 

1 .0 
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Section 5 
EXAMPLES 

Sample Problem I : Flow of a Viscous Incompressible Fluid in an Axisym- 

metric Container 

The container geometry and liquid height of the problem is shown in 
Fig. 7. Material properties, time increment, iteration step and other param- 
eters are given in the example of Section 4.1. Using the LHMAC2 program 
the first second of the flow under the following g-level is simulated: 



g a = 11.4 m/sec 6 (0<t _<_0.7 sec) 

= 0 (t>0.7 sec) 

Fig. 8 shows the flow and velocity fields of the fluid at selected times. 

Sample Problem II: Flow of a Viscous Incompressible Fluid in a Rectangular 

Container under an Arbitrary g-level 

The flow of a viscous incompressible fluid in a rectangular container 
is used as a sample problem of the LHMAC3 program. Many cases are 
studied to show the influence of viscosity and g-level on a flow. Figure 9 
shows the geometry of the problem. Parameters which are used in model- 
ing this problem are given in the example of Section 4.2. Some of the cases 
which have been investigated are listed in Table 5. The flow and velocity 
fields of the fluid are shown in Fig. 10. 
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Table 5 

VISCOSITY AND g-LEVEL USED IN VARIOUS CASES OF SAMPLE PROBLEM II 
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Section 6 

CONCLUSIONS AND RECOMMENDATIONS 

This research effort resulted in the development of analytical tools for 
he study of propellant motion in tanks during a docking maneuver in f'pace. 

The axisymmetric flow computer program is capable of accurate simulation 
of propellant dynamics in realistic shaped containers for docking loads aligned 
parallel to the tank axis. For off-axis loads the two-dimensional program 
provides a reasonable indication of propellant motion and the forces and 
moments on the tank wall. 

The pilot 3-D program can be used for the study of propellant dynamics 
in a rectangular shaped container and can be developed further to study prob- 
lems of a more general nature. 

Propellant dynamics problems to be encountered in the future space 
shuttle flights will be mostly of a 3-D nature. The lack of symmetry is due 
either to the tank geometry or the acceleration vector. Development of a 
computer program for investigating these problems is needed. Recommen- 
dations of continued research effort to extend the LHMAC3 program for 
studying the dynamics of liquid propellant in commonly used tanks are to: 

• Include the capabilities of handling curved boundaries and 
arbitrary initial free surface into the LHMAC3 program. 

• Refine the current 3-D MAC computing technique for analyzing 
flows having waved free surface. 

• Develop a new numerical scheme of including the surface 
tension effect for flows in a low-g field. 

• Study the roles of viscosity and tank geometry in flows under 
a wide range of g-levels. 
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and the Velocity Fields of a Liquid 
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Fig. 3 - Twelve Possible Cases of Two Empty Cells Neighboring 
with a Free Surface 
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NOTE: Length is measured as number of cells 



Fig. 5 - Variables to be Used in Preparing the Data Deck 
of the LHMAC2 Program 
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Fig. 6 - Ceil Status of Semple Problem latts 1 sec 
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t s 0.25 sec 





t a 0.5 sec 


Sample Problem 1 — Flow and Velocity Fields at Selected Times 
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Fig. 8 
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t = 1,0 sec 


Fig. 8 - (Concluded) 
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Fig. 10 - Sample Problem II 


Flow and Velocity Fields at Selected Times 
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t = 0.3125 sec (Case 2) 


Fig. 10 (Continued) 


t = 0.5 sec (Case 3) 


Fig. 10 - (Continued) 
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t = 0.34375 (Case 4) 


Fig. 10 - (Continued) 
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t = 0,375 sec (Case 7) 


Fig. 10 


(Continued) 
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y a 1.125 m t = 0.5 sec (Case 1) 



Fig. 10 - (Continued) 
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y a 1.625m t = 0.375 sec (Case 7) 

Fig. 10 - (Concluded) 
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Appendix A 


LISTING OF LHMAC2 PROGRAM 


I 


Lifting of LHMAC2 Program 
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*RUN,//T 

* m^g ,n i-orKHrrn-HuMT^vii tr ?n mat program itapf gen 00027 2) 

• A r >G « T l-HMAGP • T . GAV/rOG . LOCK HEEPMAC , '’PROGRAM 
, print no ihmagp, 

• for, | G rn \ , \ 


r 

r |,or KHH f/HUNTGV II LIT pf) MAT UROGRAM ( I. MMAC?, 1 >OK CORK. GRACf ) 


COMMON/I 1 / A 1 IS PR* rv* PT, IMAR, IPHM, JPAP, l ,r • , 

1 Al(IPOn>« API I J'OO) , API »P0(I> , A4(li"O0), AAI )pOO)« Af I ) POO ) « 

p A7MPOO), AGIIPPO)* AO ( I POI * i « AHUlPnoj, Atl(IPnO), AJPIJPOO). 

'» A 1 "» (1 POO ) 

f HU | N/AI I NCf I A 1 ,P ) * I AP , U ) « I AM » V ) » < A4 ,UT IU > • I A' > « VT 1 1 I * I AO«R'-i! > * 

K I A 7 , THl TA ) « I AO* 0 ) « I AO t KT ) « I A 1 O « I ) * I A ) 1 « H ) - I A 1 ? , M ) , I A | 3 , U ) 

INTEGER A | , AO, P , 00 
iMTPr.ro tv or ( pp ) 

PAT A TYUF/1f»M O fM PAPER , pO»AH / 

1 FORM A T I P 1 r -> • PPH. p , P IP* PB.P ) 

PALL f PI 'M ( 0 ,TVPP) 

jn ppAp 1 , JMAP, JpAR, PR, D7» PT , IPHM, PC, ALP 
tnpu trap*' 5 
JPPt: JPAR+.3 
PALL rtKOUT 
IFIfnAR) po,?o,pn 

PO CALL PROCI I P? « JP? , P ,U , V ,UT I L , VI I L , Pb I , THETA , D ,K|- , G, H »P, Q > 

GO TP IP 
OP CALL emPJOB 
STOP 
FNO 

• FOR, IP cn->,c,n? 

SUBROUT I NE PROG I 1 P? , JP2 ,F,U,V,UT1L»VTIL,PS1 ,THETA,D«NE ,G»H,P«G) 
COMMON/L!/ ALP, PR, DZ, DT , IBAR, 1PHM, JBAR, PC, 

1 A | ( 1 POO ) , A?I1?00), A3<1?00), A4I1200), A^llrOO), A6II200), 

P A7 ( 1 POO ) , ABIIPOO), A9I1P00), A10I1200), A11I1200), A12U200), 

•7 A | 7 ( 1PP0) 

D IMFMS l ON F ( l PP , JP? ) , U I I P2 , JPP ) , V I I P2 * JP? ) ,UT IL ( I P2 , JP2 ) ♦ 

1 VTILI I P? , JP2 ) ,PS I I IP?, JP?) , THETA I 1P2, JP?) ,D( IP 2 , JP?) ,KF< IP?,JP2) « 
p C'’PP,JP 2 ), H(|P?,JP 2 ), P I I P2 , JP? ) , OI IP?, JP?) 

DIM,. jS I ON BCLTI60), BCRTI60), NAME I 1 2 ) , R I 66 ) « RIP I 66 )' RM0RPI66)* 

1 ROS('A), RPORM < 66 ) , RR 1 I 66 ) , PRP < 66 ) 

DIMENSION GAMMA I BO), XBDRY I BO), YBDRY < BO), ALPHA I BO ) , DL ( 80 ) , 

I DBI 1 P 8 ), FNI12RH FXIlpB), FYI12B) 

DIMENSION I SFGRR ( 3B ) < ISEGRZI36), IDUMPl?), JHYBI30), LHYBI29), 

I nhvBI?), VNIRO), lEMP(pO), JEMPI20), GLVLTT I ? I ) , GRT I AO ) ,GZT I AO ) 
COMMON/L?/ RC00RDI3B), ZC00RDI3R), XP I 10000), XI, Y1 
COMA.ON/L3/ SPIA 8 P), DS, STC« STH, STR, ST? , STR2, ST2?, STR22, 

1 ^T-?PP, NSP, LSFG(F), LSFGSIB), SMC I 0,0) 

COMMON/L 4 / SPTI2), DUMP! R, 60), NPRT3, CYCLE, 1NDSMP, LERROR, 1 SUP 

PFA( MU 

INTEGER A|, AP, F, PC 

INTEGER CYCLE, TYPE, BND , EMP, FUL , LPB , OB, SUR , HYB, HSUR 
DATA BND, FUL, SUR, EMP, OB, HYB , HSUR , LPB/ 1,2,3,4,6,6,7,10000/ 

2 FORMAT! I H 1 , 1 PA6) 

7 FORMAT ( 4F4 • I »RFP»3) 

4 FORMAT ( 4F! 3,4 1 1 0 ) 

c; FORMAT! ' IBAR='I3,» JBAR=*I3,» DR=»FB.3,* DZ-'FR. 7 ,» PT=* 

1 FR*e,» | PHM= M 7 , • PC=» t’’ »• ALP=*FB. C ) 

6 FORMAT I • PCB~R“T-L= »4< FB, 1 , IX) * IX, *A= • FB • 3 , 4X , ' B- 'FB«3« 6X , • C = • 


»> 


I 


A- 1 
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07= »E9. 3) 
Ffl • 3 * 

TYPE*: i 12) 


XC = ' F8« 3* OX * » YC= 'FB«3* * 

n.3) 

system* ) 


XD = • 


GWn • 1PF lO.r. 


I Ffl • 3/ * NU=* F9.3,4X* *FPS"»Ffl.O* r 'X* 'GR^'F9.3. • 

7 FORMAT ( * T«» Ffl. 3, • TWPLT"* Ffl. 3* • TWPRT" • 

t ♦ TWFIN- *FA*3» * LPR=* IP* 

p FORMAT <21R,AFP.3) 

9 FORM AT !* NX - *11*' NY" M3** 

1 FB*3»OX* ' YD= »F£l.3« • 1,10" * F fl # 3 » ' V0=*Ffl*3) 

in FORMAT! *0 NX=0* 1 OX* IB* ' PARTICLFS IN 

11 FORMAT! *0 PART I CLF STORAGF FXHAUSTFO') 

1? FORMAT! t Hi * 1 2A6* 4X< ' T" ' 1PF12.B* ' CYCLE” * I • 

1 3 1 FORMAT ! AX * M *4X* »w»*6X« »F( I * J) »flX* *U! I « J) *BX* »V! I « J) '7X. *PSI ( I « J) • 

1 By, • THETA! 1 * J) *AX« • D( I « J> • ) 

14 FORMAT ! 3X* 2 ! FX « I 3 ) * I 9 ♦ EX* R ( 2X « 1 PF 1 ? *F ) ) 

1R FORMAT! 1M K ) 

lft FORMAT!/// • NSGMTS = ♦ IP* BX* ' NJC1 = ' 13* ^X* • NJC2 = I* 

- * 13 , fx« • NPRT2 = • IP**!/* *DEPS = ' Ffl«2«BX» ♦ VEPS 

rp.p// * 1SUR =' I?* AX* ' STH = • F6*2* 

• STZ =' F6 • 2 * BX* ' DS *• F6*3* SX ♦ 

« 1CYCLF = • 13* 4X * • I PLOT = • 12* 4X* ' COFST = • F9.6// 

« PHO = • F 1 0 .» 4 * BX* ' THCKNS = ' Ffl *2 * AX* • NGLVL = • 13/// 

• VALUES OF RCOORO ( I ) * ZCOORO ( 1 ) * GLVLTT ( I ) * 

7P!l)« LHYR(!) AND NHYB ( 1) •> 

1 -7 FORMAT!// !1MA)) 

FORMAT! PFIO.3) 

FORMAT!' T- * 1 Pfll 2#B* ' CYC.LFs'Ifl*' ITFRst 17) 

FORMAT! 1 0 I R « ?Ffl • 3 ) 

FORMAT! '0 CYLINDRICAL COORDS ALLOWS NO INFLOW') 

FORMAT! 'O NO ORST OR I/O ROUND AR I FS * ) 

TYPF = • 12 « ' LI * • IF' 


1 

2 

3 

4 
«=; 

6 


• LQUDHT 
FS.2* AX* • DBFTA 

* <^TR =• F6.2* AX* 


AX* 

: • 

AX* 


GRT( I ) * GZT ( I ) * JHY 


1 « 
lO 
2n 
21 
22 
23 


FORMAT! • 


• L2= ' I fl * ' L3= • I fl * ♦ L4= • I 5* 


1 

2 


Lfl= • Ifl* ' 

I iRs • FP • 3 ) 


L6='Ifl*' 1.7* • I fl » » NXB='Ifl*' NYB=»I5*' UL=»F8.3* 


24 

2R 


• BX* I 4* F 1 0 »B ) 


FORMAT! // ! 1 ftFB *3)1 

FORMAT!/// ' **** DIMFNSION CHECK **** 

PFAO 2* NAMP 

READ 3* F)CB,BCR«BCT*BCL*A,B*C*NU*FPS*GR*G7.VSCALE 
RFAO 4* T* TWPLT* TWPRT, TWFIN, LPR* NPRT 3 

RpAn 20* TYPF*Ll*L2*L3*L4«Lfl*L6*L7«NX8.NYB*UL.UR 

1F(P<" .FO. O) R0L=l*O 
2* MAMF 

R * I BAR « JBAR*DR *DZ *DT * I PHM * PC * ALP 
6 « BCR « BCR * BC T * BCL « A * B * C * NU ♦ EPS » GR * G7 
7* T* TWPLT* TWPRT* TWF I N »LPR » TYPE 
IF ! TvPF .NF* 2) GO TO 27 
READ lfl* NSGMTS * N JC 1 * NJC2* 


PP I NT 
PRINT 
PRINT 
PRINT 


LQUDHT* NPRT2 * 1SUR* I CYCLE* I PLOT* 


1 NC| VL 

MGl.VI 1 = 2*MG| VL 

MOL V! 2 = NO! VI + 1 

1NDO s 1 

MJCFI l = NJC 1 + NJC 2 

PFAn IB* (PCOORD! I 1 , 1 = 1 * NSGMTS ) 

RFAn lfl* ( 7C00RD ( I ) * I = 1 * NSGMTS ) 

RFAO lfl* DFPS « \/EPS * DBFTA* SIGNVN* STH* STH* STZ* DS* 


PHO* 

THCKNS 

prAn 

1 

( Cl. VL TT < I) * t = 1 * NGLVL2 ) 

RFAn 

1 

(OPT! M * 1 = 1 *NC-l.VL 1 1 

PFAn 

Iflt 

( C-7T (I ) * I = 1 »NOLVL 1 ) 

pfAo 

♦ 

( JHYR! t) , 1 = 1 . N JC=T.L ) 

RFAr 


« LHYP ( I ) * I = 1 * NJCFLL ) 


COFST* 


A-'i 
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* 


PF An IP* ( MHYO ( t ) * I - 1 * N JCFI.I- ) 

PRINT tR« M r >GMT e >* NJC1* NJrp, LOUDHT « NPRT?* DFPS» VFPR* DPI-TA* 

| IRuRi 5>TH • STR * RT2 « D r > * ICYCLF* I Pl-OT* COFGT* PHO, THCKNS • NGI-VU 
print * i proopnil ) * f n 1 «n c .gmtr 1 

PP f N't 'Vl « ( 70OORn ( | ) « I - 1 « N r -RMT r . ) 

PRINT PR* ( 01 VI..TT I I) ♦ I * I «NGLVLP ) 

PRINT PR* < CRT ( I ) « I n I «N0LVL l ) 

PRINT PR* < P7T ( I ) * J " 1 *NRUVI 1 > 

Pfi ronM^TI// (JRFIP.AM 
TOAGr n P 

PRINT IT, t JHYP ( I ) * ! n 1 »N J0F|. L ) 

PRINT IT, (| HVPI I ) * T “ I «Nv.Tri L> 

PRINT IT, (MHVP( I ) * I a 1 *Njr r l ,U > 
nRFTA a PBFTA/P7 
I rnRAP = R 

p-r comttnnf 
y=P 0 

Y = nP#FLOAT( 1-PO) 
no Pn 1=1*1 PP 
p jp< t ) =v+y 

RRP ( T )=1 • /P 1 P ( I ) 

P( T)ePTP( I )-.**Y 
pni(T)si,/P(l) 

•7 = 4. WP( I ) 

PMOPP ( T ) = ( 7-Y ) / ( 7+Y ) 

PPOPM (!)=!« /PMOPP ( ! ) 

pn y= y+ v 

IP 1 = tpap+i 
JP 1 = JRAR +1 
t PHM- I PHM + 1 
POQ = 1 , /HP 
pppp_pnR*Rnp 
Pn7= 1 ./P7 
pn7T=pn7*Rn7 
ppon- 7 =np*pn 7 
n7onn=n7*ppp 
pnDn-7=pnP*Pt)7 
OOF 1 -t,*NU*R0P 
OOF7 = t,*N( l*P07 

>>(=( 1 , + ALP)/(T.*(P n P7+P n 70) ) 
nTonp=nT«-pnp 

nTonT=nT*POT 

TPLT = TWPL.T 
TPPT = TWPPT 
T'HPl. T = TPLT + T 
TWPPT = TPPT + T 

rvoLf* = P 

T = T + A, AAAOAl 

C TO *FT UP THF MAROINFS OF THF GFOMFTRY OF THF PROB 
YRs IPAP*00 
yt= jp AP* n7 
VF = I PAP 
YF = JPAP 

jppp r a 1 OPO/AMAyO ( I BAP* JRAR ) 

KI1 s P* I PAP 

IF(Pr ,FO. n» I RPPO = 1 0?T/AYAyn (Ml ♦ JPAR ) 
pppr) | = p| nAT(IPPPO) 

IPActp = |PrtP»IPPor 

A-3 
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JR A FTP n JRAP# JRPPF 
1 V|. ~ MOf>R - TRA r TP)/R 

| VP rr IVL. + TRAOTR 
|VR rr (lORR - JRAGTR)/? 

IYT - !vn + JPARTR 
no ■’n J^I.JPP 
no RO I r ) , 1PP 
r( f « j) ^ a 
U ( ! ♦ J ) n 0,0 
V< I « J> " o,o 
IITTI(T.J) - 0.0 
VTILfl.J) - o.o 
P e 1 ( t 4 J) B 0.0 
THPTAf 1 » J) = 0,0 

-»o n( i « j) = o.o 

r to fpfc. try fipoy ffllo and hyp cfllf 

no Ri jsi.JP? 

F ( 1 ♦ J ) = 1 

01 F<TP?.J1 = 1 

TF(TyPF .FO. ?) GO TO rr 

no i = ?. i pi 

F( T . 1 ) = 1 

R? F< 1 . jPO) = 1 

•v* FONT 1 Ml IF 

t F ( TvP p «NP. ?) GO TO R7 
T FMP ( 1 1 a TRAP 
JFMP ( 1 ) = 1 

Ml = ? 

on rc m=i«mjcfll 

N? = I P 1 - LHYR(N) - NHYR(N) 

1 F(N^ .LT. 1 ) GO TO RP 
t F MP (Ml ) = N? 

JFMP (Ml) = JHYP(M) + 1 

N 1 = N 1 + 1 

RF FONT! Ml IF 

IFMPfMl )rr 1 PAP 
JFMP (Ml) = JPP 
MPFWPF a M 1 
TNHFMP a O 
I PR a I PAP + R 
JH = IP?/? 

IH1 a 1H + 1 

NGGFTH = nggmt c /p 
JH = JO?/? 

COFR = F.?PO?*nP 
IFCPO .FO. n GOFRsTHCKNG 
VH = JRAO/? 
yH = t F AR/? 
y m a 1 P A P 

1F(PF ,FO. O) yHaO , 0 
no RF K|= 1 .MJFFLL 
J = JHYR(M) + 1 

N 1 = LHYP(M) + 1 

N? a LHYR(N) + NHYP ( i 
F( IPo. J) = FMP 
1 F ( P<* .FO. 1) F ( t . J ) =FMP 
no 1 =N 1 ♦ M? 

1 F ( PF ,rn. 1) F ( IP1-1 . J) =HYF J 
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06 F( ! « .n ~ HYP 

rjpo - cjp#*p 

c.T 7P - 6T7**P 
GTPOo n 6TPP/GT77 
r.T7Po n ^T7?/ r ~TP? 

OTC r 6TH + 6TZ 
77 rONTTMIJF 

y n ft,'>H |c .OBT(AR^(rtR'*yf 5 ) 

Y a O • 77#6QPT ( APS ( G7 ) *YT ) 

DROll r. AMAV1 I V * Y «l J|.. ♦! IP) 
pn 7n ,J= 1 < JP? 

PfLT( J1 =nri. 

POe.( j>sl #o 
TP. RroT ( .M = n CP 
K=1 

mp=o 

IFITvPF.^O.l) GO TO inn 

print ">? 

C. TO AGBlGN MAPKFR PART I CLF6 AND I N I T I AL I 7F U(I«J> AND V<1«J> 
40 RFAp p, NX # NY # XC *YC « XD # YD « l lO«V r 
IFINV.FO.O) GO TO flO 

PRINT 9i NX « NY # XC# YC » XD# YD# UP# VO 
XTF = 1 * /NX 
YTF=1 ,/NY 
YsVTr* . 

nROU = AMAX1 (DROIJ, ABSdJO ) , ARG(VO) 1 
|F(TyRF .FO. ?) GO TO 141 
YG=Xr*PnR 
YC = Y-~*PD7 
xn=Xn*PnP 
YD=YD*PD7 
pn x=yTe-#.6 
SR FONT TNI IF 

I F < X # GF • XC « AND • X#LF. • XD • AND • Y • GF • YC • 

* AND.Y.LF.YP ) GO TO 60 
GO TO 70 
60 I = Y+t • 

J = Y4- P # 

TF(F( T # J) #F0.06)G0 TO 70 
X° ( X > = y 
YP(k>n =v 

V =F + 7 

I F ( K .GT. LPB) GO TO 76 
NP=NP+ 1 
F ( t ♦ J)=FUL 

IF( F ( 1 + 1 « J) .FO.BND I U ( I # J > =U( I # J ) *ROS 1 J ) +U0* ( 1 «-WOS ( J) ) 
IF < F ( t + 1 « J ) ,N C #BND « AMD# F ( 1 + 1 # J ) #NF • OB ) • ) ( I » J > =U° 

|F (F( 1-1 « J) #NF*BNP# AND#F ( I — 1 # J ) • NF • OB ) U< 1-1 ♦ J)=UO 
I F ( F ( 1 # J-l 1 .NF.BND. AND.FI I *J-1 > .NF.OH) V ( I » J- 1 ) =V0 
IF(F( I , J+1 ) .NF.BND) V < t » J ) »V/0 
70 V=V+VTF 

!F(X,t.T# TBAP1GO TO 
y=y+vtf 

!F(Y,l.T# JBAP) GO TO 
GO To 40 
76 PP TNT 1 1 
RFTl ION 

H n TFIDPOU #<;F. 0.000001) DROU= V6CALF*DR/DROl ) 

\ A-5 


u c 


LMSOHREC D 225632 


PPTMT 1<~i» MP 
A r >«:!rM pon to kpft 
r,n Tr> M r ' 

TO COMPUTE THF RASTf: R POINT COUNTS ASSOCIATED WITH THF GFOMFTRY OF AN 
I NFLOw-OUTFLOW AND/OR OH PPOB« AND TO DEFINE OH CELLS 

jon i p ( pr «fo» n #and# ijl «gf« o»ooii so to i ?o 

PRINT P3« TYPE i I 1* LP* L3» LA « LB « L6« L7« NXB« NYB« UL ♦ UR 
IF (L7iTO«H) GO TO tOG 
I L st. n+ ^ 

!P=|..A +1 
JTF =| 7+1 
no 1<"iA JrF,sJTF 
no 1 nA I = I L ♦ I R 

lOA F ( I «J)=OR 


TL 1 


I 1# I pppc 

4* 

IYB 

IL? 

« 

LP^lPPPr 

4- 

iyb 

IL3 

~ 

L3*I PPPC 

+ 

IYB 

I LA 

ss 

L 4* I PPPC 

+ 

IYB 

tL* 

- 

L cs*TR ppr 

+ 

I XL 

IL6 


l**i pppc 

4- 

I XL 

IL7 

- 

L 7* I PPPC 

4- 

IYB 


XOTS s 0 # B/AM AXO ( NXB « 1 ) 

vnis = i . o/amaxo(nyb» i ) 

YFIP=< • S* YO I S ) +L 1 
WIN = MYP* (LF-L1 ) 
roLS-n,o 
no i io j=i » jp? 

IF ( J.GF» <L 1 +? ) • AND. J.LF. (LF+1 ) ) 00 TO 107 

106 IF < J«CF« ( L3+? 1 • AND* J«LF • ( L4+ 1 ) ) CO TO 109 
GO Tfl 1 1 O 

107 U< 1 « J)=UL 
BC.LTI J)=-l .o 
GO TO 106 

1 no u( IPi , J) =UR 
RCPT( J)=-1 . 

IF (UP .GF. o.nni ) GO TO 110 
ROG( J)=0. 

«COT ( J1 = 1 • 

1 1 o CONTI NUF 
GO TO AO 
l?n PRINT PI 
PFTI IOM 

C COMPUTATION OF RASTER PT CNTS« INITIALIZATION OF VEL COMP AND ASGMT OF 

c marker particles for flds in an axisymm tank with ellipsoidal bkhds 

1 A 1 CONTI NUF 

HO 161 N= 1 ,N c OMTS 
I CEGDP (N ) - RPPCl L*RCOOPn(M) 

I SFGP7 ( M ) = RPPOt.l *7COPRO(N> 

I SFGPR ( N 1 a ISFGPR(N) + 1 XL 

161 IGF GP7 (Ml s TGFGPZ(N) + I VO 
I GF GoP ( NGGMTG + 1 > = I SFGRO < 1 ) 

I Sf 'OP7 (NGGMT C + 1 ) = tSFGR7(11 

Ml = o 
N? = 1 

HTI 0| m = 1.0' PHI 
NG( iMt 1 = NFGMT C - 1 

171 Ml = Ml + 1 
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A c c , I r.N 177 TO KRTT 

7 - nrnnpniM?) - Prnnpn(MI) 

ill - oronon ( mo+ | ) « pronpD(MO) 

1 F ( Adp ( 7 ) # OT • 0.0001) OO TO 170 

OM 1 7 A TO YPFT 
vn r prooon (Ml) 

OO To vpft 

170 yr n (?roopn(M0) - 7oooRn(Mi ) ) /? 

y r r 7 ooopn(Mi) - yr#prooRn ( ni ) 

!F<Anr (III ) »l.,T • 0.000 1) OO TO 1 77 

ynio a < 70 oopn < no+ i ) - yronon (M ^ ) ) /ul 

yoio r 7 ooopr » ( mo ) - yni c «proopp ( mo ) 

1 77 vn ® ( Y - Y O ) /xr 

7 = (Y - YniM/ynlo 

TF(Y .OT. 7O0OPO(M0) ) YP= AVAX 1 ( YP » 7 ) 

17A OONTTMMF 

y = 0 . 6 *XTF 

I F ( JOUR .OT. o .AND. Y .GT. STH) CALL FCTN? < X * XTE . Y ) 
y a v + YH 

t F ( X .OT. yn) OO TO 177 
J s o.o + Y 
176 | = 0.0 4- X 

yp ( v ) = y 
yp ( <+ n = y 

y = y + f 

!F(Po .FO. O) OO TO 170 
yp(y> a vw - yo(y-o) 

XP(F+n = Y 
y = y + o 
NP = mo + 1 

170 I F ( K .OT. LPR) OO TO 76 
MP = KIP + 1 

IF(F(I«J) .FO. HOUR .OR. F(l»J> .FO. FUL ) GO TO 177 
V ( t«J) = VO 

IF(F(! + 1«J) .MF. W) U(t»J)=UO 
!F(F(T«J) .FO. HYB) GO TO 176 
F( I , j) = FUL 
OO To 1 77 

176 P( T ,J) = HOUR 

177 V = v + VTF 

tr<y ,|.t, yn- on to 176 
y = v + YTF 

1F(Y .or. HTI.OUP) OO TO 170 
1 F ( y ,|.T. 7OOOPn(N0)) GO TO KR^T 
OO TO 171 

1 70 IF(PO .FO. O) GO TO tBl 

moomti = m c omth 

Ml = IPO + 1 

OO loo J=1»JOO 
F( tm, j) = r ( i h+ i . j ) 

IF( CT(JH,J) .FO. HFUR .OR. F<1H.J) .FO. FUL) U(IH.J)=UO 
OO 1 OO 1 a J HI « T P 1 
IK TP7-1 , J) r lllliJ) 
v( fPR-1 « J) = V ( T 1 J ) 

IPO r ( IP-l-l , J) r F ( 1 , J) 

1 p 1 Ml = o 

MO - 1 

Ml = 1 
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o f on = 1,0 

MO n n 

1PP MO = MO + 1 

too c|ON n "1,n#r|GM 

MRHA| P !• N1 f 1 
|F (MO , OF , ^>1 OO TO 1 BO 

vnnPv(MOi r proonn(Mp) 
ynnPY(N7) n 7rnopn<NP) 

MO r Nn + 1 

NO o MPOMTH + 1 

1 RO Ml n Ml + 1 

NP a N° + 1 

TF(M'» • OT« MPOMT c ) OO TO 107 
I F ( Pr ,FO, 1 .AMD, NP ,OT» NO) OO TO 107 
7 s ornc'Pn(Mp) - RononniMi) 

T F ( An e ( 7 ) ,LT , n.onni) r,o TO 1P7 

T F ( e T ON , L T » 0,0 , AN'H , 7 .IT. ) OO TO 1 HP 

ipo vr s ( 700000 (N?) - ?ooopn( ni ) ) /*' 
vr s 7FOORO (Ml ) - YC*P000RD< N I ) 

yn s proopn(MP) + o,onni*PiON 

N4 = PrOODD(Nl ) + 0,0001 

1 F ( e. I ON ,0T, 0.0) N4=RCOORO< N1 ) +0.9009 

7 = N4 

1 P7 7 « 7 - PI ON 

IF (P TOM) JOO,1PO,1P<5 
1RO IF (7 . LT. VO) 00 TO 107 
OO TO 1°1 

ion IF (7 ,OT. YO) OO TO 197 
19i yPnPv(MP) = 7 

YRDRv (N.7) b 7*XC + YC 

M7 s K|7 + 1 

OO Tn 1 P7 

ion yppPy(M7) = POOOPO ( NP ) 

YROPy(N7) s 7FOOPO ( NP ) 

N7 s N|7 + 1 

OO Tn 1 OF 

107 NRORvP = NO - P 

1F(NA . LT • P) NnnRYP=NRORYP— 1 
I F ( no . FO ♦ O) on TO ?0 7 
MRnPvP = p*NRDOYP 
M s Kl7 - P 
J = N7 -1 

no poi 1 = 1 « n 

J = J - 1 

yRDRv i N7 ) a xui - ypnRY(J) 

YRDRv ( M 7 ) a Y«DRY(J) 

701 NO = M7 + 1 

P07 CONTINUE 
Nl = 0 

HO P11 MbUNPDPYP 

Ml = Ml + 1 

M? e Ml H 

y o = ynnPY(MP) - xpopyimi > 

1 F ( Aop ( yo ) ,OT • o.nooo) 0,0 TO P07 
Ml b Ml + 1 

MP b Ml 4 1 

yr = ynnPY(NP) - xrory(M1 ) 

P0 7 YC = YPOPY(NP) - YRORY(Ml) 
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alpha CM) = apscyc/xc ) 

ALPHA (M) - ATAM ( A| PH A (N ) ) 

ni c m ) n sqptc < xc*op>**o + <yc*d7)**0) 

|F<vr ,ot. n,m r.o to r<m 

yr - AP r .<yrvYO) 

fAMMACN) " AT AM ( VC ) + 7» 141ft 

IFCYr .IT. O.P) GAMMA <N> nft.Omp-ATANCXC ) 

rn to o i n 

000 yr r a p f ' < yr/yr ) 

r.AMMA(M) rr ATAM(YO) + 1.0700 

1F(Y0 .LT. 0.0) OAMMACN) '\ ,R70R~AT AN < XC 1 

01 n yROPY ( N ) ff 0,<T#( ynoRY (Ml) + XRORYCNP)) 

Oil YRDPY(N) = 0.0* ( YDOPYCNI) + YRDRYCNP)) 

I F ( I H|P .LT, 1 ) 0,0 TO ?1« 

CALL AeOOMP(VH) 

IF (Pr ,i-0. O) 0,0 TO 01 R 
< *« 1 

|. s 9*N‘-P + 1 
m = r**MOP - 1 
y s rnAP 
no Pin 1 = 1 , mop 
opt ( y ) = y - op(Y) 

c--.PT (1^+1) = OP(M+1 ) 

OPT ( | ) S OP CIO 

OPT Cl. + 1 ) = SPCK+1 ) 

K = y + 0 
L = L + 0 
o 1 7 M B M - 0 

NFP . t3*M.e;p 
* = 1 

no oi A 1 = 1 « NOP 
SP < K 1 = CPT<Y) 

CP( K+1 ) b OPT (K-+1 ) 

O 1 A y = y + 0 

01* CONIT T NI IF 

IFCMPP-MPC .OT. PO) PP I NT Plft 

Plft FORMAT ( /// • *•#•«•* ERPOR — NPEMPC EXCEEDS PO IN OO-LOOP 35« «) 

PRINT ?]o, N'PFMPC , ( I FMP C I ) , 1 = 1 * NPFMPC ) 

PIP FORMATC /// • VALUES OF NPFMPC., IFMPCI) AND JFMP(I) • 10X. * NPFMPC 

I =. \p// ( 1610) ) 

PRINT 17, C JFMPC I ) , 1=1 , NPFMPC) 

NRDCl I = NRORYP + JRAP - NJCFLL 

IFCPO ,FO, 1) NR0CLL=NP00LL+JRAR-NJCFLL 

PPIN- Pin, (ALPHA! I >, 1=1 , NRORYP) 

PIP FORMATC// • VAHJFS OF ALPHA C N ) AND DLCN) • // C16FR.3)) 

PR t NT PA, (Dl.C I ) , 1 = 1 .NRORYP) 

PRINT u?n, NRORYP, CXRORYC I ) , 1= 1 ,NROPYP) 

PP O FORMATC/// • NBDRYP =» 14// • VALUES OF XBDRY ( 1 ) , YLLDRY C I ) AND GAM 

1MA(1) • // (16FR.P)) 

PRINT oa , ( VRDRYC I ) , I =1 , NRORYP) 

PRINT PA , C 0,AMMA C ! ) , ! = 1 « NROPYP) 

PFPIYT - FLOAT CJPAR - NJCFLL) *07 
IF ( Ro ,FO, 1 ) PFR!NT=P*PpR!MT 
no pr>p 1 = 1, NRORYP 
00 0 PFPJwT = PFO»mt + OI.Cl) 

NI = 1 

NO = R 
NIP = 1 
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ma - Ninnri I 
no <- m jr 1 + ? 

M7 r JP1 - NJPP 

irtPr .ro. n npmhh-i 
nn< n n o.r*dl < l ) 
rM<(Mnn<"U-) >- PPPfMT - nnm 
no pof. JrP.JPI 
no r*O0 IfMP,IP1 

V r t 

!F(J .OT. M7) KMP1+NP-I 

mo = r<i* «j> 

00 TO < «?P4 ♦ r’?A * PP6« PV3* ) « NO 

,y?7 mi - mi + i 

NT e NT +1 

K|4 e M4 “ 1 

nn(N*»> = drint-i ) + o.r*<ol<N1-i ) + 0L<Ni)> 

|F(Pr .ro. 1) no <N4 ) BPFRIMT-nn<N3) 


r 

c 

r 


on to ^of, 

IF(K .Mr. ! P 1 .OP. J .LT. NO .OP. J . OT . N7) 00 TO ??0 
NT = NT + 1 

N4 = N4 - 1 

DH(MO) - OP(NT-l) + 07 
!F(Pr *F0. 1) OB ( N4 ) =PPR I MT-DR ( NT ) 

ojaft FONT TNI IF 

TF(Po ,FO. O) 00 TO 40 
N 1 = Ml + 1 

N = NRDCLL/P 

PRINT = DR(N-I) + O .<=;*( DL ( N 1 - 1 ) + DL(Nll) 

DP < N+ 1) - OR(N) + 0.R*<DL(N1 ) 4- DL < N 1 + 1 ) ) 

FALL OLFOUT 
OO TO 40 

TO T NO J CATC Fill, AND FI IP OFI..LO 

COMPUTATION OF FORCER AND MOMENT MOD MAy BE NEEDED IF DR NE DZ» 


MOMFNT I f COMPUTED AROUT TANK GO. 

PRO FONT I Ml IF 

IF(T .OT. OLVLTT ( INDO+1 ) ) INDP= INDO+1 
Ml = 9*1 NDO 

TEMPI s 0|.V/LTT( INDO+1 ) - OLVLTT ( I NDO ) 

C,P = ORT(Nl-l) + (T “ OLVLTT ( 1NDG ))*< CRT ( N 1 ) 
07 = fiZT(Nl-l) + (T - OLVLTT< INDG) )*<GZT<N1 I 


GRT(N1-1 ) ) /TEMPI 
OZT < N 1 — 1 ) ) /TEMPI 


no Pkr j=t«JP? 


no ?c«=: 1 = 1 , TP? 


P( I « JT = O.o 
OI I . J) = 0.0 
?RR KF ( t . J > = O 

I ND = 1 

INDFV = 1 

I TPP = O 
ITFP\/N = O 
NPT = 0 


K - 1 

pR-7 I=yP(K)+?. 

j=yP(K+l ) + ?• 
KF( t . J ) =K 
K = K + o 


npt=mpt+1 

IF(NPT.LT.MP)00 TO ?R7 
I F ( TvP c .NF. ?) OO TO ?FO 
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nn prp j~ i , jpp 

|FIF(|P?(J) .ro, BMP) CO TO P r >0 
F ( I Po » J ) n FMP 
F F ( l O p » J ) n O 
!F(Pr ,rn, n) on TO pnn 
F( 1 »,|) ~ r MP 

FFI ) , J) n o 
pr-P CONTINUE 

l mo nvn • i mo cmp + i 
po pro k h 1 ,MPrMPr 
V - , JFMP(J) 

Ml n I pp - |PMP(.n 
nn pro I rN|] ' ! P 1 
n pi/) n run 
FFI!,^) r o 

IF (Pr , F 0 . O) GO TO PRO 
Ft IP-7-1 .F ) s FMP 
FF< I PP— I ,F ) B n 
PGn rPMTTMHO 
PAG FONT I Ml IF 

C TO modify thf vfl comp of a newly created FMP OR HYB cell 
no paf j=p«jpi 
no pa<= 1=P«I°1 
MO s F t I « J ) 

GO TO t 2AB.P6B.2A1 .2AB.PAG.PAR. PA3 ) « NO 
PA 1 l F t F F ( I ♦ J ) • NF . 0 ) GO TO PAG 
F ( I , J)=FVIP 
GO TO OAA 

PAO t F ( Ff < I « J) .NF. O) GO TO PA«= 

F( ! . J) = HYP 

P64 TFtF(l+l«J) .FO. FMP .OR* F(I+1«J) *FQ. 

J F t F ( I— 1 . J ) .FO. FMP .OR. F(l-l.J) *FQ. 

I F ( F ( t « J+ 1 ) .FO. FMP .OR* F(t«J+l) •FQ» 

J F t F ( I ♦ 1 ) .FQ. FMP .OR. F(l«J“l> »FQ. 

?6 C CONTTMIJP 

C TO OFF IF A SUR-CFLL SHOULD RFCOME FUL.OR A 
DO J-^P.JPl 


HYP) Ut l . J) =0.0 
HYB) U< 1-1 , J)=0.o 
HYB) V( I . J)=0.0 
HYB) V( 1 iJ-1 ) = 0 • 0 

FUL-CELL SHOULD BECOME SUR 


no p*?G I=P, IP1 

IF(KF( I i J) .FO. O .OR. F ( I ♦ J ) • GE • OB) GO TO 270 

Ml = n 

|P(fr ( i + i,J),FO.FMP.OR.F( I-l^J) .EQ.EMP.OR.F ( I ♦ J+l ) .EQ.EMP.OR. 
1 Ft T » J-1) .FO, FMP) N1 = 1 

IF(F< I+t .J) .FO.HYB .OP. F( 1-1 . J) .EQ.HYB .OP. F < 1 , J+ 1 ) . FO . HYB 
1 F( I . J-1 ) .FO.HYP) N t = 1 
MO = F ( I ♦ J ) 

GO TO ( P70.PA7.PAP) « NO 
?A7 l F ( M 1 ,FO. 1) F(t«J)=cuR 


. OR . 


00 TO 7 70 

PAP IF (Ml .FO. O) Ft I « J ) = FUL 
?7 o CONTINUF 

ASSIGN PRO TO KRFT 
t F ( TvPF .NF. P ) GO TO ABO 

I F ( I r UR »LT. P .OR. LFRROR .NF. 0) GO TO A^'O 
CALL. CURVE 

tF(l e-RBOP .NF. O) GO TO A c 0 
MPT . ?*NFP - 1 


k- 1 - n 


WO - n 


A-ll 
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K"t n 1 

no P'7 f- Kn| , MPT « P 
| - rpiK) + P,0 

J n rP(K+1 1 + P»0 

jr < | ,ro«K j .amp. j »ro. Kp) on to P 7 G 
mo - r < 1 , J ) 

OO TO <P7A,P7A,P7A,P7r.,P76,P70,P7A ) , NO 


07« K - 1 r. J 


►•0 r ,) 

r ai.| roper < v ,r i ,fp , r 7 »r a , ka , oti • , r '.TV ) 
ir<K/i ,ot, n) n rju*rnrAT 

1 r ( ht ,ot, o) 0(1, rn> r r lTW »rnrrT 
r»7o comtimhf 
M - rs'ri r 

1 1- ( N , r 0 » o .OP. N.KO, 0 ,00, N,h0,40 »0R. N.FG.GO .OR, N,EQ»00> 
1 r.o to 077 
OO TO oon 


0 7*» PPtNT 070, GYCl.F _ , _ , 

07P. FORMAT! // • **#* CHFCK P(!,J) AND 0<I,J) ****UOX, • CYCLE 8 

no Jcl, JPO 


U > 


M = J 

070 PRINT 1AB4, (P(I,N),!=1,|PP)» ( 0 < I ,N ) , I = 1 , IP? > 

O o To A c ' n 

PP.o I F ( ry(~LF ,NF .0 ) GO TO 700 
INITIAL otiPrACF PPFFCURF 
nn ooo J=0, JP1 
no poo I=IPMM,IP1 

IF ( F ( I , J ) ,F0, SUR 1 THFTA ( I , J I 8 ( A+P*C0S (C*( J-I«5)*DZ) ) /OT 


pnr FONT I Mt IF 

COMPUTATION OF PPFUDOPRESSURF OF SUR, HSUR, 
inr I CHF ok' s 7P0 

no iio j=?,jpi 
nn onn I =P , I p 1 

I F ( I NJDFX ,FO, P ,AN0, F( I , J) • NF • HOUR) 
IF(lMOP-y , FO • 0 ) OO TO 0 1 R 

l F ( CYCLF • NF * 0 .AND. F(I,J) .NF. HSUR ) 


FUL AND HFUL CELL? 


OO TO 7?0 

THFTA( I , J)=0*O 


->r> 1 


tF (P( I , J) .FO. RUP )0.o TO nni 
IF(F(I,JI .NF. HFtlP) GO TO 71 P 
M = 0 

JF(F( I + 1 , J ) .FO. FMP .OR. F(I + 1,J) *FQ. HYD ) 

I F ( F ( I , J+ 1 ) .FO. FMP .OR. F ( I , J+ 1 ) .FQ. HYP) 

jF(F(!-1,J) .FQ. TMP .OR. F ( I -* 1 « «J ) .FO. HYP) 

tF ( F( I • J“ 1 ) .FO. FMP .OR. F ( I , J- 1 ) *F0. HYP) 

IF ( NPRT7.EO. 1 .AND. N.EQ.O .AND. CYCLE, LT»3) 


N=N+ 1 
N=N+? 

N = N+4 
N = N+R 

PRINT 656, 


I , J, ICHECK 


!t=-(M ,rn, o, on to 017 
GO TO <.7n^,310,?pO,.70P,7?0, 


7?0 , 720 ,310, 3P0 , 720 , 320 , 320 , 720 , 320 , 320 


1 ) » M 

7nc; THF TA ( I , J ) r THFTA ( 1 , J ) +0 OF 1 * ( U ( I , J) “U ( 


1 , J) ) 


GO TP 7RO 

3 1 o THFTA ( I , J)=THf-TA( I , J ) +C0F?* ( V ( I , > “V ( 1 , J- 1 ) ) 

GO To 7 . on 

7 ic. 1 F ( F ( I , J ) .M*-*, FUL) OO TO 7?0 

tFITvPP" ,FO, 1 ,AND, UR ,L P • 0,001) GO TO 720 

717 THETA( I , J) - G7* ( J—l 1 #DZ •+• GR*( I — 1 , F > *DR 


GO TO 7Pp 

710 no 7,o k=i ,Kpnn>yo 

Ml 7 yonPY(F) + .7,0 
M" 5 - yooioYtY) + ?«0 
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r 


IF (I .ro. N1 .and. J .ID. N3 •AND* I DUMP < K ) »FQ. ?) THLTA(I*J)= 
1 Turr a ( I « 1 1 + r I r-Nv/N*nnrT a*wm ( v ) 

T|n FONT TNI If 
7pn FONT I Ml II 


irUnpry ,ro. 

Abb frN 770 TO 


p ) an to i )On 
K-Prr 

thf , ji i r .T"Oi it' i nr 

r 74 4 
.1 ’P..IP! 

I r - P # I P 1 

ji ,ro»‘ i ip )Fo to 


rr.Tl \/n comp 


,>r> 1 


, I) *NF, HMIR) DO TO 744 


TO FOMPl itF 

mo | o nr ok 
no o/i/i 
no 'M A 
irifii. 

I r ( r < j , 

T‘l 1 N - O 

I F ( F ( 1 + 1 * J ) || <)• IMP •OP • f ( 1 + 1 . J ) il u. HYP) 

I f < f- ( i , j + j > .ro, r mp •op, i (i.j+n .mi. myb) 

!i < r ( i-i , j > .i -<>, i mp •op, i (l-i,.i) oFO, iiyp) 

IF < F ( I , J- I ) .ro, FOP .OP. F ( I » J“ 1 ) »FO, HYP) 

IF ( MPPT3 .FQ. 1 .AMO, N.fO.O .AND, C YCL.F .LT • 3 ) 

IF(N ,r-0. O) r.O TO 044 

4.0 TO ( 737 , 737 * 344 « 3.13, "HOD * 33.7, 344 « 333, 333 • 333*344 ♦ 333 ♦333*333* 


N = N+ 1 
NuNT, 1 
f l -M+4 
N-N+O 
PRINT 


640* I ♦ J* I CHECK 


1 344 )* N 

333 I F < F ( I « J+ 1 ) ,FQ • FuL .OR. F ( I « J+ 1 ) .FCI. bllR .OR. F(I*J+1) *FQ» 

1 HbtlR) OO TO 77* 

40 TO *» 3 7 

334 GO TO I 33b * 337 * 344 * 33t> » 33b « 337 » 344 « 3 37 » 33B « 337 * 344 « 336 « 33b « 337 ) ♦ N 

33b IF < F ( I + 1 * J+ I ) ,FQ, FMP .OR, F ( I + 1 * J+ 1 ) ,FO, HYP) 

* V < I + 1 * J ) sV ( I « J ) -OPOD7* ( I ) ( I ♦ J+ 1 ) “H < I ♦ J ) ) 

J F ( N ,FO, « .OR, N ,F0, 13) 40 TO 336 

GO To 44 4 

336 I F ( F ( I - I * J+ t ) • FO • FMP • OR • F ( I - 1 « J + 1 ) ,FO, HYB) 

* v( 1-1 « J)=V< I » J)+PR007*(P( 1-1 . j+t )**u< l-l « j) ) 

TF(M ,NF, 1?) 40 TO 344 

337 IF(F(I+1,J) »FO« FuL .OR* F(I+1»J) ,FQ, bUR .OR. F(I+1*J) .EQ, 

1 H<MiP) 40 TO 330 

40 To "*44 

339 GO TO < 344 * .34 1 ♦ 344 ♦ 344 * 344 ♦ 34 1 * 344 * 343* 344 ♦ 34 1 » 344 * 343 * 344 ♦ 34 1 ) ♦ N 

34 1 1 F I F ( I + 1 * J+l ) .FO. FMP .OR. F ( 1 + 1 « J+ 1 ) .EQ. HYB) 

* U( I * J+l )=U( I * J)-P70DR*( V( 1 + 1 ♦ J)“V< I ♦ J) ) 

1 F ( N .FO. lO .OR. N ,F0* 14) GO TO 343 

OO To 744 

34.7 1 F (F ( I + 1 « J-l ) ,FO. FMP .OR, F < I + 1 ♦ J- 1 ) ,F0. HYB) 

* U< 1 , J-l ) =Ut I * J)+0700R*( VI 1 + 1 ♦ J-l )-V( I ♦ J-l ) ) 

74 4 FONT INI IP 

1 F ( TvPF ,NF, ?) CD TO 749 
N 1 = NJF 1 + 1 
VI 1P3*N1 ) = V( TP) « N 1 ) 

1FIPO ,ro, )) V( 1 *M1 ) =V< 3*N1 ) 

N 1 = JP1 - NJC? 

VI I Po *N 1 ) = V I I P 1 *M1 ) 
t F ( Pr ,FO, 1) V ( 1 ♦ N 1 ) -V ( 7 ♦ N 1 ) 

C REFINFMFNT 1 6 NF.CF6SARY IF PROB ON PROP RF-OR lb TO BE SML TO* 

GO TO ( 74b«74ft,740) « IFAbF 
74 c VI 14*17) = VI 17, 17) 

U ( 1 7 , 1 7 ) c -V I 1 7* I 6)*COb< ALPHA! 1 « ) )/S INI ALPHA U ^ ) ) 

IFIF( 13*17) «FO. C MP) P( )?* 17)=UI 13* 17) 

U(17«1B) = ~V I 1 3 « 1 7 ) *COS ( ALPHA I 1 6 ) ) /b I N I ALPHA ( 16 ) ) 

!F(F( 1 ?« 1 ft) ,ro, FMP i U< 1 •’« 1 B) =111 1 7 * 1 P) 

11(13,10) = — V/ ( 1 7 , i p ) *FOb ( ALPHA (17)) /b I N I ALPHA (17)) 
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inrniii'i) .ro. hyp) un Moimi ipuo) 
r,n Tn i/io 
->4 0 mNTTMiir 
T/|0 rOMTtMlir 

irtryor »r r> • . r * > no in in;* 
no vo i -p , tr ap 

t M 1 • i ) -I ii I . •'* ) *nm 
'•o-.o i it i , iP '» ) -i ii | . mi ) nor t 
•u ,'i fOMT i Mi 11 

nn M ,,ipi 

ir< Tvnr .10. p . ami > . r*( .io. i . At ii » ■ i (i.i) *• ii. imp) on to mom 

V< 1 t.DnVP, • ) *1 *n T( .1) 


1 r |T. 

II ITYPI ' 

• i'Ui , 

* 1 • 

Midi r ( 

i p;s 

, 1 ) .ro, IMP) no TO MOO 


V< IP r '«.l) 

v< 

MM 

< , 1 ) M \C\J 

i < )) 


inn 

FOMT I Ml M 
F.O Tn t 1 ‘ 

IT 





07 '• 

FONT I Ml II" 






Aon 

IF ( n vF| I 

• ro 0 

n ) 

C\n Tn 

n 1 n 



H'tT'.'PI T 

• 1 r. 

T ) 

r.n *rn 

Anc\ 


An \ 

IF ( TWPPT 

• i * 

T ) 

on tn 

A07 



IT ( MOM r 

• t o 0 

1 ) 

PP1MT 

: >r i 

UIYP(PP). ^FOOPP(PP) 

An? 

JFJ Tli'F I M 

• t r t 

T 1 

PFTI IPM 


AHA 

FVFI r ■> 

ryrt p 

4* 

\ 




T s T + nr 
r,n Tn ‘.nn 

4 no tmiPLt S TWPLT + TPI.T 
iOA a^cion on 1 tn K'Pf t 
no Tn 4 pr 

40*7 twPpt = twppt + Tprrr 

onq AFnirM 4 n 0 TO KRft 
no Tn 4 nn 

41*= A'C|rM 4 in TO k'PrT 
no Tn 4 po 

41 F A^SIr.M 404 TO KRPT 
nn jr\ 4 *»n 

4?o Ann inN 4P4 TO kpi 
r Tn Pi OT THP APOMPTPV OF A PPOP 
4?1 CALL HRAMrvt?) 

tF(Tvor ,fo. n nn to 4 pp 
I F ( TvPF . fq, P) no TO 400 
n ALL l I NF V < I XL . I YP ♦ ! XP « I YP ) 
CALL l t NFV ( !XP. I YP. ! VP. I VT ) 
FALL l. T NFV ( ?XO » t YT « I XL . I VT ) 
r At. I ( TNFV ( t yi. . t VT . I VL . I VR ) 
4 0o nn to K'P i 
4^p FONT f Ml IF 

(“ALL I !NFV( IVL . 1 YP. tl-P. t YR) 
FALL I TMFV1 TI.R. f VR, !LR» It 71 
F ALl l 1 NFV ( 1 1. c , U. 7 . I L A, , 1 1. 7 ) 
FALL 1. I N'F\/ ( 1L n« II .7. IL6. I YR) 
FALL I TNFVt II P * t YP. IXP. ! YR) 
FAI.I | I MFV< IXP. I VR. I VP. It 7) 
FALL. I IMFVt IXP. II 4. IXP. IYT) 
FALL l INF\/( IXP, f VT. IVL. IYT ) 
FALt t INFV/I |VL. IYT. tXL« TL?i 
FALI 1. TNFVt I VI .11 t.IXL.lYR) 

FO Tn 400 
4 J?4 NPT-n 
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»* r- or 

k - , 

I - o 
MP! r: MP 

Apr I V 1 r Vf * ( K ) #PI ‘f T| I 

I v i - yr*(K+n hoppo. i 

|V| n I VI 4 \ Y\ 

! V 1 H J V t + I VP 

iriTwpr ,rn, o ) no to a \ o 

I - o 

l r ( vp (•-')•! t • o * n ,oi>, vi mi- ) .im .xr »om* vimh i )«i u')»n .op. 
i vp(i> i i ) .pi .vr > i.- i 

r.n to aim 

•TUi II I f Vf I I .<YI. I') (All MAUKTR ! I Y( I I *K»I, .M'TGMl 1 »M*XW) 

•T 1 • > i PUT l Ml ll 

Jill ,|’U, O) (All I >1 PPM I y 1 . I Y I * 7' >* U ) 

IP (| ,!-(), 1) MPU MI M-I 

y ' •* t- o 

MP 1 Ml»N I 

|l (MHT.I T.MP) r.n to Apn 
IP (I op ,fo, n ) 0.0 TP PJt T 

A Or | r.M a OO m FP 1 
IP ( TvPr ,Fn,p ) no TP A7P 
no to API 

r f or? plotting yn. vfctOP 
d?e L, ^ o 

po AOP JsP.JPl 
PO A OP l - P « ! P t 
M n = I- ( I . J ) 

ro TO < ApB , APT, APT, APB, APB « APR. A 07 ) . NO 
A ?7 tvl f (FLOAT! I ) - 1.n)*PPPCLL 

IVl - (FLOAT! J) - 1.0) 4RPPCLL 
IYP - O.F* (ii( 1-1 , J) 4- II! I . J) )*PROU*PPPCLL 

lY? r P.F*(y ( 1 , J-l ) + V! I , J)) *OROU*RPPCLL 

ivi = ivi + ivl 

I V 1 = IVI 4- IYP 

IXP = IYP 4 IY1 
IYP = IYP 4- 1 Y 1 

I F ( TvPF .NF. P .OR. F(I.J) .NF. HSUR ) GO TO A29 

V) s FLOAT ( 1 ) - 1 .F 1 

Y 1 s r|. OAT ( J) - 1 .f ), 

CALL MARKFRICYOLF.n.L.NFGMTl .M.VW) 

IF(L # MF. O) r.o TO APB 
APO rALL I tNry ( IX 1 . I Y1 « IXP. I YP) 
x = tyo - ivi 

Y •- IYP - 1 Y 1 

Ml = CORTiVHl? 4- Y**P) 

1 F ( N 1 . GT . 7 ) C Al L ARROUI ( I X 1 . I Y 1 » I XP ♦ 1 YP . 6 . P ) 

APB FONT TNI IF 

r-O TP FPf-T 

APO lFll.OP . FO . P) GO TO KRFT 

print ip* mamf.t.cyclf. gr* mpi 

PR l Mr IP 

0 TO CHOCK THF I NO0MPRFGIH1L I TY OF A FLUID BYSTFM 
A PA no aap JJ=1*J°b 
PO AAP I - 1 * IP P 
J» JPP4- 1 - JJ 
P( I *u)=n. 
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MO r F( I * J) _ 

GO TO <43B*436*436t43R«438t438i436) • 

A36 D( I « J ) n PRUI>*RDR*<P!P< I >*U( I .J> ~ RIP< 1-1)*U< l-l « J> ) 

^VrInt^IV.’j! Full). V(hJ) » PR1(1,JW THETA < I ♦ J ) ♦ 

1 n < t « J ) 

440 CONTTNllF 

|F(TvPr «NF • P) FO TO 44 8 
no ?*»n Id ,MnnrLi 
FN ( 1 > e 0.0 
FV ( M = n « r 
32 O FY ( I | e 0,0 
DMOM a 0,0 


PTV 

« 0.0 


RTV 

c 

• 

c 

II 


Ml 

s o 


M? 

= 3 


M3 

- n 


N4 

= MBnC.tL 

+ 1 

MR 

a NBOPYP 

+ 1 

N6 

a MJC1 + 

3 

NT 

a JP1 - NJC? 

1 F ( PO .FO. 

1 ) N2= IH1 

no 

343 J = 3 < 

JP1 

no 

34? I -N? 

♦ tPl 

K = 

* t 


1 F ( J .GT. N7) K=|P1+NF- 


= r ( K , J ) 


00 

TO (242. 

3 36.236.236. 

?33 Ml 

a N1 + 1 


M3 

a N3 +• 1 


N4 

a M4 - 1 


NR 

a NP - 1 


V 

= n7*<YH 

- YB0PY(N1 ) ) 

y 

= nP*( YBDRY(N1 ) - VH) 

N- 

o 



IF(F/Ktl • J; • ir i w • n n»r~ • 

j F ( F ( < « J+ 1 ) .FQ. EMP .OR. F(K«J+1) »FQ. HYB) N=N+E 

I F ( F ( K— 1 . J ) .FQ* FYP .OP. F(K-l.J) .FQ. HYB) N-N+4 

IFIFIK.J-I) .FQ. FMP .OR. F(K.J-l) .FQ. HYP) N=N+8 

FN(Nd = RHO«THFTA(K. J) 

I F ( N .FQ. O) FN(NB)=FN(N3) + RHO*PS I < K » J > /0T 

FY ( No ) = -C.0F3*FN <N3)*DL(Nl)*C0S( ALPHA < N1 ) ) 

IF(Pr ,F0, O) FY <N3 ) =XO0RY (N1 )*FY<N?) 

I F ( Po .FQ. o ) r-0 TO PP>4 
PM(Nj/>> = RHOfcTHFT* ( IPP-K. J) 

FV(NA) = -COF3*FN(N4)*PL<N8)*«lN< ALPHA(N*) ) 

FYtNVj) = -COF WN ( N4 ) *DL < N8 ) *COF ( ALPHA < NB ) ) 

FY (Mi) = cOFO*FN (N^ ) *DL < Ml )*B INI ALPHA IN 1 ) ) 

tFfJ .OT, JH) FY < N4 ) “ —F Y < N4 ) 

034 | F ( J .0-T. JH) FY(N3)*-FY<N3) 

?3B ?F(l< n .NF? IP! .OP. J .LT. MB .OP. J »GT • N7) GO TO 24? 

KH = M3 + 1 
M/l a Mrt - 1 

Y a rV'* ( VH - FLOAT (J) + |. c ) 

FN (Mi) a PHO* ( THFTA ( X . J ) + P c - 1 ( K » J ) /PT ) 

|F<Po .PO. O) OO TO ? 4? 
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F V | Mi ) - rnF0*O7*FN!N0) „ , _ 

FN(N4) r RHO K (THTTA ( IP3-K, J) + PG 1 ( ! P3-K , J ) /DT ) 

rv(M/i) --rrriKn7-KFH(NA ) 

PAP r Tv r FT V + FY(NT) 

IF <PG .ro, P) GO TO PAP 
PTY r. FT Y ^ F Y ( NA ) 

FTY e FT Y + F X ( N ^ I + FX<N4) 

RMOM n PMOM + Y * ( FY ( NO ) ~ FY!N4)I + Y* ( FX ( NO ) + FX < N4 ) ) 

PA'- CONTINUE 

PPIMT 44 1, PMOM* FTX » FTY 

AA\ FORMAT!// * FOPCFG AND MOMENT FXFRTFD ON TK WALL DUF TO SL LIQUID* 

1 inv, • PMOM -» 1PF1P.4, * F PX = • 1PF1P.4* OX* 'FTY =♦ 1PF1 P»A 

-a /// FX, • DO ( N ) • 1 1 X ♦ • FN(N ) * 11X* »FX(N) • 11X, *FY(N)« // ) 

00 AAP M = 1 * NBDGl. I 

1 = M 

AAP PPIf'T AAA* not II* FM (II* FX ( I I * FY ( I 1 
AAA FORMAT ( OF 1 6.0 ) 

PRINT 446* ( LGFG <11*1 = 1*41* < LSFGS < I ) » I = 1 » O ) * 

1 ( ( cmf ( I *K ) *l<= 1 * P ) * I = 1 * 4 I 

AAA FORMAT( /// • **** VALUFS OF LSFG(I)« LSFGS(Q) AND SMC ( I * < I ' 10X* 

t 4 J f , 10X* OJF // (OF14.4I) 

A4P FONT t Ml IF 

tF(CvCLF «GT* O #0R* TYPF *NF* P) GO TO KRFT 
IF (MOOT? * NF • 1 I GO TO KPFT 

TCHFoW = 4 04 

PRINT 400, CVFL.F* ICHFCK 

4 c;n F ORMAT(/// • **** CHECK CPU TIME ^OR PLOTS AND COMPUTATIONS OF DO- 

1LOOPF 600 AND 044 • OX* • CYCLE =» I?» » ICHFCK =• 14* • ****») 

GALL CLKOIJT 

4 OP FORMAT! // ( 1 4FQ #4)1 

404 FORMAT! 14FO. 4) 

HO 4 «f J= 1 , JO? 

Ml = J 

PRINT 40 ?, (U( I »N1 I * I = 1 * IP?) 

4 OF PRINT 404, (V! I *N1 ) , I =1 , IP?) 

GO TO MPFT 

ago no 4*t/i i = i , mggmt g 
N = I + 1 

CALL LINFV! ISFGRR! I ) , ISFGRZ! I ) ♦ I SFGRR ( N ) , I SFCRZ ( N ) ) 

464 CONTI NUF 
GO To 4 ?P 

470 IF! I p| OT *NF • 1 *0R* I SUR »LT» ?) GO TO 421 

*' = 1 
L = o 
MPT = o 

474 TVl = p P ( M ) «pnpGl l. 

I V 1 r GP(h"+1 ) *PPOGl..L 
I y 1 e I vi + I YL 
ivl = I VI + I vp 
GALL PL.OTV! 1X1 « I Yl *4-7*0) 
ur s ig + ? 

MPT = MPT + 1 

IF (r ,i P T .IT. MOP) GO TO 474 
GO To 4 "> 1 

r~ T n GOMPt 'TF THF Tit of \/FL GOMP 
f.po 1 F ( TvPF #FO« ?) GO TO 14 10 
i'll GOMT I N't ir 

no t J rr * t JP 1 
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no noo T r.o , jp i 
MO rr r ( l , J ) 

r,o to ( f-n/, , r-,rv 4 . npo, non , r ,°n . 504 ) « NR 

nn/, if I I ,ro, I PI) no to 5 in 
m i r ( I + 1 « J ) 

ro to ( 5 i o « f .nri « 506 . 5 l o « 5 j o « 5 i o • 5 O 6 ) ♦ N l 
eo/\ 7 | Pni i ( I « J ) *RRP I I ) *(R< I ) '*» K 1-1 <t j ) ~R ( ! + 1 > * , ) < t + 1 « J n 
I )VF) D-Jp . pf-tt ( II ( ! < J ) +i M I « J” 1 ) '*<V< 1 « J- 1 ) +v < ! + 1 « 1 > ) 

IJVTPcO.pn* (III T . J ) +U I I » J+ l ) I *< V ( I ' J) +V * J + 1 * 1 1 

l IT II ( I , J) rill ! , J ) +DT* ( PORK ( 7 1 P+THPT A ( I . J ) -THFT A It + 1 4 J )) 

* +RH7-M ( l IVPR-I IVTP ) +CR+NI I* ( R D73* ( U< I ♦ - , + 1 ) + l 1 ( I . J" 1 ) “?• *U ( I ♦ J ) > 

* -pnnD7 « (VI 1 + 1 • J ) -V ( 1 + 1 « J- 1 ) “V C I « J ) +V ( I 4 J~ 1 ) ) ) +P I I » J ) ) 

rio mi s FI I 4 J+1 ) 

ro To (Fpo,<^l p,5i P,‘ r 'Po,«-,po«^?0«Flp) , N 1 

rr^o 7 1P=\/( T *.J)*(V( T 4 J-1 )-VC I 1 J+1 ) ) 

t jv/TL srR I P ( I - 1 )+n. P'S* ( U ( I “1 4 J ) +u ( I - 1 « J+ 1))*(V( 1-1 *J> +v * ^ ' 

UVTR=PIP(t) * 0 .?^tUII.J) +U(l«J+n) *('./( 1 + t « J)+V( I « J) > 

VT IL ( 1 « J ) =V ( I « J)+DT*< POZ*(ZtP+THETA< I » J) -THETA! I » J+1 > ) +GZ 

* +Rnp*pp !<!)*(( U VTL-l IVTR ) “NU* < R I P ( I > * C PDZ* I U < I 4 J+ 1 > “U < I 4 J ) ) 

* -RDP*<V<!+1.J>-V<UJ>>> -RIPH-U *(PDZ*<U<I-l.J+l>-U<I-i«J>> 

* -R no* ( \/ ( ! i J ) — \/ CT — lfJ))))) + ^ ( I ^ J ) ) 
fOMT T Ml JF* 

r to rowpt it^ d— of ful and hful ofilf 

n O *=^OC? J=r p , Jp 1 

UT I L ( I P 1 4 J ) =UTIL ( IP1 » J)*ROS< J)+UTIL< lPl~l « J)*< 1 .-ROSI J) ) 

| CHF r/ “k' = K «° 

OO 5„o J=P,JP1 
OO p-/)0 1=0, IP 1 

NO = F ( I , J ) 

CO TO ( 535 , 534 , 535 4 535 , 535 , 535 4 53 1)4 NO 

. ... .,4 41-4 4 /'/O <fr\ croc 


!F(F( I+l 4 J) « 
I F ( F ( 1 4 J + 1 ) « 
IF ( F ( I— 1 ,J) « 
I F ( F ( I , J-1 ) . 

IF (NPRT3 • FO 1 


, OR # F( I + l »J) .FQ. HYB) GO TO 535 
. OR 4 Fit 4 J+1) • FQ. HYB) CO TO 53<c 
.OR. FI 1-1 4 J) .FO. HYR) GO TO 535 
.OP. FI 1 4 J-1) .50. HYB) CO TO 535 


IFIMPPT3 .FO. 1 .AND. CYCLF .LT. 3) PRINT 655 4 I 4 J 4 ICHFCK 
504 4 ( ] ,J) = RR III) *PDR* I P T P I I ) *UT I L I I 4 J ) - R I P I I ~ 1 ) *UT I L ( I - 1 4 J ) ) 

| + POT* (VT It. I ! 4 J) - VTILI t 4 J-1) ) 

ro to 5^0 

50*= DC| ( t « J) = 0 . 

*=4 5 CONTI Ml IF 
CO^PUTATTOM OF PC I ( I ♦ J ) 

IFIlMPCy .CT. 1) CO TO 547 
ICHBric = PCI 

jr (TyPF « FO • P) CO TO B 10 
04 c; 1FI tMPFY . OT • P) CO TO 643 

IFIImOFY .FO. P .AMD. t TFR » CT . 0) CO TO 350 

K4 7 TOHrok- = 570 

ee^o 1 F ( TvB5 .BO. P) CO TO r,r,7 
OO 1=3* TP1 

PCI I T 4 1 ) =POT I I 4 0 ) 

DC J ( I 4 JP3 ) =P C 1 I t 4 JP 1 ) 
ere;-? rOMTtMIlF 

no j= o 4 JP 1 

l F ( F ( 1 4 J ) .FO# BMP ) P5 I I 1 4 J ) =P5 t ( 3 4 J ) 

550 lFIFlIPF.J) .FQ. BND ) P5 I ( I P3 ♦ J ) =P5 1 ( IP 1 » J ) *P05 I J ) 

T F ( 1 MO.ro. 0 ) 00 TO 600 

IMHrrn 

ITro-iTrn+i 
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HO 570 jr-p,jpt 

no | ~p , T F> 1 

MO r- 

r.O TO ( f, 70 , , r,7n , r,7n , r^n , n70 , nfi? ) , NO 

667 1 F ( F ( ! + 1 , J ) .FO. FMP .OP. Ft 1 + 1, J) .FQ. HYP) GO TO 070 

l F ( F ( 1 , J+ 1 ) .F : Q. FMP .OR. F ( I , J+ 1) .FO. HYP) GO TO 070 

! F t F ( 1-1 , .J ) .FO. FMP .OR. Ft 1-1, J) .FO. HYP) C.O TO 070 

IF(F( I , J-l) «F0. FMP .OR. F< ! , J-!> . FO • HYP) GO TO 070 

IFCCyCM .LT. 6 .And. ITFR .FO. 1) PRINT 656, 1, J, ichfck 

=67 PGIRrP'- I ( 1 + 1 , J) 

p c io=p r ' 1 1 1 , j-n 

poll -PT! ( T-1 , J1 

I r t TvP r ^ ir • 1 ) GO TO =66 

T f ( F t T + 1 , J) .ro.Pp) P c > l P = p G I ( I , J) 

IF(F( T , J-l) .FO.OP) p O i n = p r, T ( I , J) 

IF ( F ( T — 1 , J 1 .FO. OP 1 PS l L = PS I ( I , J ) 

=6= X=W#(-D( I , J1+RR1 ( I )*RDP?*(RIP< I )*PSIP+R1P( I “1 ) »FS I L ) +RDZ7* < PS I B+ 

1 PS! ( 1 ♦ J+l 1 ) )-ALP*PS! U . J) 

Y = AOc ( X ) —APS ( PS 1 ( I , J ) > 

7=Ape(y)+APS(PSI t 1 , J) ) 

PSI ( T « J) =X 

1 = (7 .l c ’. o.nnoooi ) GO TO 670 
!F( Ap6 ( y/ 7) .GT.FPS) 1ND=1 
=70 FONT t Ml IF 
GO TO =60 

r COMPUTATION OF U(1,J) AND V(I,J) 

600 CONTINUE 

DO 6on J=7,JP1 
DO 6oo j -p , IP 1 
MO = F ( I , J ) 

GO TP (670,604,604,670,670,670,604) , NO 
604 N 1 = F( 1+1 , J) 

GO TO (610,606,606,610,610,610,606) , N1 
606 U(I,J) = UT IL ( I , J ) - PDR-MPSI ( 1 + 1 , J) - PSI(I,J)) 

610 N 1 = F( 1 , J+1 ) 

GO TO ( 670,61 p , 61 7, 670,670 ,670 , 61 7 ) , N1 
61 7 V(1,J) = VT I L ( I , J ) - PD7*(PM ( I , J+l ) - PfU(NJ)) 

670 rONTTMUr 

ASSIGN 600 TO KRFT 
N 1 = APS( VN( 1 1 ) 

no 601 M=7,NPORYP 
N7 = APF(V/M(N)1 
IF (Mo .GT. Nil Ml =N? 

671 FONT INUF 

I F ( I tFR .LT. 7500 .AND. Nl .LT. 40) GO TO 675 
PRINT 677 

677 FORMAT ( /// • ***# RUN IS STOPPED BECAUSE OF ITER OR VEL EXCEED TH 

IF G I WFM LIMITS (STMT 677). *) 

PFTl |r>M 

67F ITPPx/M = I Tf*P\/M + 1 

GO TO 6=o 

600 ASSIGN 640 TO K'RFT 
GO TO 070 

640 I F ( TvP F .TO. 7) GO TO Pin 
64 7 .ASSlr.N TOO TO KRFT 
645 IF(TYP5.NE. 1 )GO TO KRFT 
DO 647 Jr 7, JP 1 
no 647 I- 7, TPJ 
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I 


l.J( I » J) =1 l( I , J+1 )*BCB 

v ( 1 1 J)=v< I - 1 . J)*bc.r 
V( I « J)=V< 1 + 1 , J)*BCR 


A SUP OP HSUR Crui 


• FO* 

FMP 

.OR. 

F( 1+1 ♦ J) 

• FO. 

Hvn) 

N=N+ 1 

• FO # 

FMP 

.OR. 

F( I . J+1 ) 

• FO • 

HY P ) 

M=N+? 

• FO# 

FMP 

.OP. 

F ( I - 1 . J ) 

• FQ* 

HYH) 

N=N+4 

• FQ • 

FMP 

.OR. 

F< I , J-l ) 

• FQ# 

HYF) 

N=N+B 

1 • A ND ♦ N 

• FQ. 

o .and. cycle. 

LT • 3 ) 

PRINT 


**** CELL F(I.J) IS TREATED AS A HFUL 
• . IN DO-LOOP* |F, *.***#*) 

no TO 6RO 


6S6. I i J » I CHECK 
CELL « I =• 1 3 » 


IF(F( I , J) , Mr. OB) on TO BA 7 
f F( F ( 1 . J+1 ) .Nf .OB.AND.FI 1 + 1 , J) .rO.OU) 
l F ( F ( 1-1 , J) • NF • OR • AND • F ( l « J+ 1 ) .FO.OH) 

IF<F( J+1 . J) .NF.OB.AND.F ( I . J + 1 ) .FO.OB) 

64 t FONT 1 Ml IF 

on x n fc'PP'T 

TO OOMPUTr THE VFl COM- OF 
6B0 t CHF rtf r RRO 

no Boo Jr?,J«1 
HO BpO Jr?, IPl 
IF(F< I * J) .FO.BUPJGO TO BBS 

jf<f(!«j) .mf. h*sup) no to bro 

BBT Mr. n 

I F < F ( I + 1 • J > 
tF(F( I . J+n 
t F ? F ( 1-1 , J) 

I F ( F ( I . J— 1 ) 

IF (NPRT3.FO. 1 
6*56 FORMAT < / * 

1 * J a* n« 

|F(M .FO. O) 

GO TO (661 . 66? . 663 ♦ 664 . 66 1 ♦ 665 ♦ 662 . 667 » 668 . 659 » 6E8 * 669 « 66 1 .662. 
1 661 ) * N 

BBR V(I.J-I) = 0,B*V( 1-1 « J-l ) 

V( t . J) = O . B#V ( I - 1 » J ) 

661 U<I*J)*(UCI-1*J)*RIP<I-1)-W<! )*ORODZ*<V< I . J)-V( I . J-l ) ) ) *RRP ( 1 ) 
OO TO BBO 

6RO V(I.J-I) = o,B*(v/( 1-1 . J-l) + V ( I + 1 « J— 1 ) ) 

66 P V ( I . J 1 = V ( I » J— l )-DZODR*RRl ( I ) * ( U ( I . J ) *R I P < I ) - J < I - 1 « J ) *R I P < I “ 1 > > 
OO TO 6RO 

BBS U( 1 . J)=U{ 1-1 . J)*PMORP< T ) 

OO TO 666 

664 U ( I - J * J ) a ( U ( I ♦ J ) *R I P < I ) + R ( I ) *DRODZ* < V < I » J ) -V < T « J- 1 ) ) ) *RRP ( I - 1 ) 
OO TO 6P0 

660 U( 1-1 ♦ J)=UC I « J)*RPORM( J) 

666 V< I ♦ J)=V( I . J-l )-.?R*D7*<U< I « J)+U( 1-1 « J) )*RRI( I )*< 1 -PC.) 

OO I O BMP 

667 V ( I . J- 1 ) = V ( I i J ) +D70DR*RR I ( I ) * ( U ( I ♦ J ) *R I P ( I ) -U ( I - 1 « J ) *R 1 P t I - 1 > > 

OO TO BBO 

66 R U ( I . J ) =U ( I - 1 « J ) *PMOPP ( I ) 

OO TO 670 

660 U < I - 1 . J ) =U ( I . J ) *RPORM ( I ) 

670 V( I . J-l ) =V( I . J ) + .? c 5*DZ*RRI ( I )*< U< I .J)+U( I“1 . J) )*< 1 -PC) 

6R0 CONTI MUF 

OO TO K'PFT 

700 MP T = p 

I F ( TVPP .MF. P) OO TO 700 
! F ( 1 <M IP .I T. P) OO TO 700 

r - 1 

KM = 1 

NPN s o 

701 I F ( NOT .OF. MSP) OO TO 70R 
in = SP(K) + 7.0 

HP v = in - i.o - SP(K) 

HMy r 1 . n - HPY 

JP = SP(K+1 ) + 1 . F 

HPY = JP - O.s - sp ( K+ 1 ) 

HMy = 1.0- HPY 


»* 
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( IK* - HPy*HMY*U < I 0“ t » JP+ 1 ) + HMY*HMY*U( I D , JR+ 1 > 

) + hPY*HPY*U ( I D— 1 , JR ) + HMY+HPY+l I ( f D » JR ) 

|P a cp<<) + l ,F 

HPy rr TR - 0,<n - C P(K) 

hm y „ t , n - hpy 

jn = «p( y+n + p » n 

HP v - jn - i#p - gp(K+i) 

hmy = i.n - hpy 

V K - HPy*HMY*V< !R» JD) + HMy*HMY*V( IR+1 «J0) + HPX*HPY*V< 

1 + HMY*HPY*V< IR+1 « JO-1 ) 

c;P(<m) = spc*) + uk+otodr 

«;P(KM+1) = SP(K+1 > + VK*r>T007 

YC S FP( K+ 1 ) 
no 700 M=l«MC,r,MTl 


1R, JD-1 > 


MO a M 

IF(vr • Gt • 7C00RD ( N ) • AMO • YC • LE • 7C0CRD ( N+ 1 ) ) GO TO 704 

-O? rONTTNHF 
r.n to 707 

704 7 = OPOORD(NO+t) - RCOORn(NR) 

IF(AnC<7) »LT • O.ooi) oo TO 70G 
y C. * (7C00RD(NR+1 ) - 7C0ORO ( NO ) ) /7 

yr s 700 ORD ( N Q ) - XC*PC00RD( N9 > 
yn = <«.p(y+n - YC)/y c 
go TO 70* 

70 ^ yn s POOOPn(MP) 

70* TF(«.P(t*') , GT « yni go TO 707 

MPM - KlPM + 1 

KN = KN + P 
7 A 7 K = fe r + ? 

N p T _ MPT + 1 

GO TO 701 
700 M«;P s NPN 

IFUcuR ,FQ« O .AND. INDSMP *GE * IC.YCLE) CALL CHKSMP 
700 MPT a O 
NPM = O 
K= 1 


YM= 1 

TO COMPlJTF THF movements OF A MARKER PARTICLE 
7ir> 1 P ( MPT * OF • HP ) GO TO 71« 
in = YD(K-)+?, 

HPY: TP-1 .-VP(K) 

HMy= 1 ,-HPY 
JB=YP(Y+1 ) + 1 * G 
HPVs jP-0,G-YP< Y+1 ) 

HM Y= 1 .-HPY 

I lK = HPV*HMY*t l( in-i , JR+1 ) +HMY+HMY 

* *U< ID. JR+1 ) +HPX*HPY*U ( ID-1 « JR)+HMX*HPY*U< ID. JR) 
lR:Vn < K ) + 1 

HPVs I P— o , c-yp ( K ) 

HMvr 1 ,.-HPy 

jn=yn(H4-i )+?, 

HPv-jn-1 ,-yP(^+1 ) 
mmy= t ,-HPY 

V K':Hrty*H«v«u( IP, jn )+HMy*HMY*v< I R+l , JD)+Hpy#HPY 

* *v< ip, Jn - 1 > +H'/y«HPY*v( I p+1 ♦ JD~ 1 > 
yP ( km ) = yP ( i< ) +t |y*n-ronR 

yp ( K m+ 1 )=yP(K+l 1 +VK+ 0 T 0 D 7 
I =yp ( ym ) + ? , 
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c 


jByPfK-M+t )+P. 

IF ( TYPE .NF. 1)00 TO 71 « 

| p ( vn ( k*M ) • OF . ! p AP ) OO TO 7PO 
71 r yNnKM+P 

NPMpMPM+ 1 

7?n IF(F( 1 , J) ,ro,FMP)OP TO 7ft0 

TF(F(f«J) * FO • HYP) OO TO 7PA 
770 Y = K+7 

MPT=MPT+ 1 
r,o TO 710 
77 c: MPwMOM 

IF(Tvpr.Njr. 1 1 00 to pfo 
740 yaULWRDR* ( T +DT ) — YD I S* < 2.#C0LS + 1 • ) 

!F(Y,l T.O. )G0 TO PFO 

TO INTRODUCE NEW MARKCR PARTICLES IN CASE OF AN 
rOLF = OOl F+1 • 

Y=VFTP 
MPNaKlP+M I N 

7 FO yP ( YM ) s Y 

yp ( k*m+ 1 )=Y 

KNssKM+p 

tsy+p. 

J=Y US 
NPsMP+ 1 


INFLOW-OUTFLOW PROS 


—ne; 


7*0 

7*4 

770 


PI o 
PI F 


Y=Y+YO I F 

TF(F( t 4 J) .NF.PMPIOO TO 7FF 
F ( 1 . J) rFUR 
lit T . J)atiL 

IF (Km • OT • L°<7> OO TO 7F 
T F ( MP « L T » NPN ) OO TO 7F0 
no TO 740 
F ( t . J ) « SUR 
OO TO 7*70 
F ( I . J ) = HSUR 
I F ( F ( 1 + 1 * J) 

IF(F( 1-1 . J) 

I F ( F ( I .J+1 ) 

IF < F ( I . J-l ) 
no TO 770 

no TO (PI^.PPl). IMDFY 
Ml =0 

HO Pol N=1.MFDRYP 
I a yppPY < M ) + P.O 

J a yODRY(N) + 2.0 
IDI IMP (Ml = O 
VN < N 1 a 0.0 
IF(F( t « J ) .NF. HFUR) 

IF(F( l+l « J) .FO. FMP 
I F < F ( I ♦ J+1 1 *FQ. EMP 
1 F ( F ( 1-1 . J) •FO. EMP 
I F ( F ( I .J-l 1 .FO. C MP 


EQ. 

EMP 

.OR. 

F( I + 1 . J ) 

• EQ. 

HYB) 

FO. 

FMP 

.OR. 

F ( I — 1 ♦ J ) 

• FQ. 

HYB) 

EQ. 

EMP 

. OR * 

FC I . J+1 ) 

• EQ. 

HYB ) 

FO. 

EMP 

.OR . 

F< I ♦ J-l ) 

. FQ . 

HYB) 


P7 1 


T F ( C VOLF • 
no To n-»i 
I F ( N 1 .FO. 
IDl)Yp(N) = 
M? a M 

FONT 1 Ml IF 

k- = 1 


o) 

1 


Ml =N 


U( I . J)*U< ID. JD> 

U( I-i» J)=U( ID-1» JD) 
V( ! * J ) *V ( I D ♦ JD ) 

V( I ♦ j-l )=V( ID. JD-1 1 


GO TO 831 

•OR. F ( 1+1 , J) .FO. HYP) 
.OR. F< ! « J+1 ) »FQ. HYP ) 
.OR. F(I-l.J) .FO. HTP) 
.OR. F< I . J-l) «F0. HYP.) 
LT. A) PPINT 1» J* t CHECK 


GO 

GO 

GO 

GO 


TO 

TO 

TO 

TO 


8?1 
821 
epi 
8? 1 
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[ 

] 

I 

! 

I 

} 

I 


MPT n 1 
N4 « 1 

MO r O 

TNnrv b ? 
no opt ! e i , o 
nn o*»7 jci,mbhpyp 

077 DUMP( | . J) n r.n 
Pin IP (MOT . r,T. MP> CO TO 007 
I B V0<K) 

j = ypt^+n 

no Pol N#N 1 .NO 

IF(tnilMP(N) .PO. O • AND* intlMP(N-n .RO. P) N4=N 
IF(lOUMP(N) .PO. O .AND. IDlJMP(N+l) .TO. P ) N4=N 
!F(|oi|MP(N) .MP. 1) CO TO R<M 
N6 b yPDPY(N) 

M7 = YPnPY(M) 

I P ( I .PO. N6 .AND. J .FQ. M7) CO TO 841 

PO TO 

P41 7 = <yP<K) - XRDRY(N) >*COS(GAMMA(N> ) + <XP(K+1) - YRDRY ( N ) ) *S I N ( 

1 GAMMA (Nil 

I p ( apc ( 7 ) * OT • DPPG) CO TO 8*1 


TOUMP(N) = 

7 


N6 a XRDRY ( N ) + 

0.0 

7 = M6 
DUMP( 1 «N) 

= 7 - 

l.o - yPORY < N ) 

OUMP( p.N) 

= 1.0 

- dump< 1 .N) 

N7 = YPOPY(N) + 

1 . e 

7 = M7 
nt IMP < 7 . N 1 

= 7 - 

n.o - vonRY(N) 

nilMP(4.N) 

= 1 .O 

- DUMP(O.N) 

Mfi s XPnrjV(M) + 

1 .* 

7 = M6 
DUMP < P « N ) 

= 7 - 

0.* - XRDRY(N) 

OUMP(6.N) 

c 

• 

u 

- nUMP(O.N) 

NJ7 a VRDPV ( N ) + 

0.0 

7 = M7 
nUMP(7«N) 

= 7 - 

1.0 - YBPRY < N ) 

D< JMP < P . M ) 

* 1 .O 

- DUMP<7.N) 


NO = NO + 1 

MA a K + 1 

7 = (XP<K) - XRORY(N) )*COS(CAMMA(N) ) + <XP<K+1) - YBDRY <N ) ) *S I N < 

1 CAmma (N)) 

IF(NPRT3.F0. 1 .And. N0.LT.R1) PRINT 844. NI,N2,N.XP<K>.XPCN4».Z 
04 A FORMAT ( / * CHFCK DO-LOOP 8*1 * 10X. 31*. 7F1P.*) 

M4 S N 
GO TO P** 

081 CONTI NUF 

ppe 1P(M4 .NP. Ml) CO TO 8*0 
N1 = Ml + 1 

POO |P(NA .NP. N? ) GO TO 003 
MO = MO - 1 

POO t P (N9 .LT. Ml ) CO TO 071 
* 8 in + P 
NPT = MPT + 1 
GO TO OOP 

067 PRINT 060. CYCLP. NPT. ( I DUMP ( I ) « I * l .NPDRYP ) 

060 FORMAT (/// ***** NO TP — THFRE IS A HOUR RUT 7 I S LT DPPS. ******* 
1* ♦// P I P/ / <1610)1 
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071 ir<cycLF .ct. 3 .op. nprtp .nf. n go to ear 

dp I kit ft7P. OVOLF « NPT* < IDUMP! 1 ) * 1 « 1 *NRDPYP) 
n 7 ? FOPMAT (/// * **** VALUES OF CYCLF * NPT* I DUMP (!) AND COMPUTED AREA 

1 f, #**# • // ?1R// (lfi!P)) 
no i t=i *o 

Ml = 1 

o 7 r, PP 1 MT ( PUMP ( M 1 . J ) • J n 1 * NRDRYP ) 

00 TO OAF 
MR 1 CONTI MtlF 

no not N=1.NRPPYP 

IF( inUMP(N) .NF. 2) OO TO R c > 1 

in = ynnRY(N) + ?.o 

jp s YnnPY(M) + 1 

IP s xnnRY(M) + I.f 

UK * DUMpTi*N)*0UMP(4*N)#U<!D-1.JW+1) + 

1 JR+n + DUMP(1.N)*DUMP<3*N)*U<!D-1.JR> + DUMP ( 2 « N ) *DUMP < 3 ♦ N ) * 

^y/K* ' = ^ DUMP ( R «N ) *DUMP ( 8 *N ) *V ( 1R* JD> + DUMP < 6 * N ) *DUMP < 8*N)*V< IR+1 • JD) 

1 + nUMP(5*N)*DUMP(7*N)*V< IR* JD-1 > + DUMP ( 6* N ) *OUMP ( 7* N ) *V ( I R+ 1 « 

P jn-i ) 

\/N(l.) s uK*COS ( GAMMA ( N) 1 + VK*S I N < GAMMA t N ) ) 

P91 CONTINUE 
N) 1 = O 

no O0 1 N= 1 * NRDRYP 

OO} IF ( A«S< VN ( N ) ) »LT . VFPS) N1=N1 + 1 

1ND = 1 

t F ( N 1 .FO. NRDRYP) t NDFX=3 
tF(MnPT2 .MR. 1) OO TO 545 
CALL Ct.KOUT 

PRINT O10. CYCLF. ITFP. ( VN ( I) . I = 1 ♦ NRDRYP > **** , ✓/ oio// 

910 FORMAT!/// « **** VALUES OF CYCLE* ITER AND VN( 1) **** • S' 218// 

1 ( 1 KFR.4 ) ) 

00 TO 545 

14in OO TO ( 1420, 1 440 ) « I NDFX 
1420 DO 1470 J=1*JP2 
PO 1 4 I = l*t PP 

G ( I « J ) = U ( I * J ) 

1470 H( I « J) = V( t » J) 

1440 IFICYCLF^OT. 2 .OR. NPRT2 »NE» 1 .OR. ITERVN .GT. 3) GO TO 1445 
CALL CLKOUT 

,44? FOr’^T ! 4 // • CHECK U(l.J)* VU*J) AND PS 1 < ! * «J ) PE FORE ADJUSTING TH 

1 FT A ( l . J) • ) 

DO 1447 J=t«J°2 

1447 PW!NT J 14R4* ( U ( I * N 1 ) * 1 = 1 * 1 4 ) * < Vt I * N» 1 ) * 1 = 1 * 14) * ( PS 1(1 *N 1 ) • I - 1 * I 4 ) 

144F CONTINUE 

DO 14^0 J=1 «JP2 
OO 14^7 1=1, tP? 

IK 1 . J) = C-( 1 * J) 

\IT, %\cyrxi .OP. NPPT? .NE. 1 .OR. .TtRVN .GF. 3, GO TO 502 

nnfKl . 1 /) -?n . rvrt !Trp» ! NHFy $ SIGNVN 
1470 FORMAt! // • CHECK CYCLE* ITFR* 1 N D FX ♦ S 1 GNVN AND THETA ( I « J ) • 10X* 

1 71F, l-A,?) 

1 a R o j=i,jp? 


ITERVN 


3) GO TO 1445 


ITERVN 


3) GO TO 502 


I NDFX ♦ S I GNVN AND THET AC^J) 1 1 OX « 
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Ml r3 

\ PRINT 1 AHA * ( THFTA < ! *Nt ) • I = 1 ♦ 1 A > » <U< I ,N I ) * ! « 1 ♦ 1 4 ) ♦ < V< ! «N 1 > « l » 1 » 1 A ) 

i/in/i format! / <mfo.4)> 
ro Tn fo? 

FND 

♦ FOP* TO eO'j.nn'i 

FI IRPOI IT T M F r|. KOI it 

CALL FCLOCKTDATf , T I MF *FGFC * F60SFC ) 
tifPITr (F, ]nnn) T T MF 
i non FORMAT ( AHOTIMFb at?) 

O**#** moot FIFO FOP FXF-C O VFPFIOM 
CALL CPI IT!M( ITIM) 

FSFC a FLOAT ( IT! M)/l .FF. 
mP I Tf <*, pnon ) FSFO 

?OOn FORMAT! 1 3HOFSFC (CPU) = F14.A) 

PFTUPM 

FMD 

♦ ftcM * I L CPt )T I M * CPUT I M • USF AS CALL CPUTIM(ITIM) 

Tl( l) aXP'T’ • WHERE ITIM IS FLAPSFD CPU TIME 

rm itim* i A AC. f23. ARRAY) • IN MICROSECONDS 


♦ asm, IL 

CPUT I M 

« CPUT 

IM • 

T. ( 1 ) 

AXPS 


• 

CPUT I M* 

LA 

AC 

♦ (23* ARRAY) 


FD 

PCT* 


1 A 

AO 

* ARPAY+?? 


MF I « 

yu A n 

,000 


SA 

An 

♦*P.X1 1 


J 


XI 1 

AP n AY 

oc-S 




FMD 




♦FOR* IS SO4*S04 

SUBROUT I NE MARKFR ( CYCLF « K * L « NSGMT 1 « M * XW ) 

COMMON/LF/ RCOORD(SF)* 7C00RD!3F)* XPUOOO' - ')* XI* Y1 
L = c 

YC = XP(F+i ) 
tF(^ .FO. P) YC=Y I 
DO 3A?n N=l*NSGMTl 
MO s m 

IF(YC • GT • 7C00RD ( N ) .AND. YC »LE. ZCOORD(N+l)) GO TO 3030 
■JAPO COMTINIIF 
L = 1 

POTI IOM 

TO TO 7 a nOOOPD ( NO+ 1 ) — RCOORD ( NO ) 

lF(AnC(7) .LT • 0.001) OO TO 3O0F 
y r a (7COOPn(N°+1 ) “ 7COORDINO) )/7 

VC a 7 FOOPD(NO) - xC*PCOORD! NO) 

YO a (VP (!<•+!) - YO/VC 


I r t* 

.FO. 

0 ) yn=r 

ro rr 

-JO/lC 


VP) s 

prOOPn(MO) 

I r 

,NF . 

C .AMD 

!F<K 

• FO. 

0 .AND 

IF( V 

.ro. 

O .OR. 

yn s 

X'M ~ 

XD 

! r CK 

.NF. 

o , AMD 

TF (K 

.FO. 

O . AND 

Dfri inM 


FMn 



f r rnr 



K) .GT. XD) L* 1 
•GT. XO) L=1 


XO) L=! 


iPPol IT I N|F A^OCMPIYH) 

C0MMON/L 3/ SP(4fl3)* DS» STC* STH* ^TR* ST7* STR?« ST7Z* T ZJ?« 
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LFFG(B), U r Fr r , (H ) , SMC<4«0) 


V* n 

n 





M rr 

1 





M r >P 

- 1 





VI 

P. t *0' 

r 

VH 


VI r 

cth 





VTP 

F* O ^ 





yn n 

err 

- 




Vr m 

rTH 



nr, 

T7 

^P( 1 

) r? 

VI 




^P( P 

) o 

Yl 




K| 

p M r 

P 

4> 1 



I F ( 

irp • 

CnT 

. ?41 ) 

no 


l?n c-o TO (3100,3140) . N 
n \ TO V p e Vt *f VTF 

y? _ ST r - FORT(FTZ? - FTZR?*<X?“XH)**?) 

7 s cQRTMX? - X|)*»? + C Y? “ Yl)**?) 

IF (7 .LT. •"'«) CO TO 31«C* 

XTF - O t Ofl*VTF 

r-o to -»i? n 
^140 y? = Y1 + VTF 

yr> a CORTCFTR® - FTR7?*<Y? - FTC.)**?) + *H 
7 = eORT ( ( v? - X|)**? + ( Y? - VI)**?) 

IF (7 « OT • oc, ) r,0 TO FIFO 
XTF s 1.0?*XTP- 
no TO -»1?0 
omo fp(wt) s y? 

FP ( K* 1 ) * v? 

IF(YO .OT. YD) RETURN 
K = *" + ? 
vi = v? 

Y 1 = Y? 

I F { N ,rn, ? • OR • Y? .LT • YO) GO TO 3110 

M SS o 

VTF s Y? - FP(K-3) 
no to Olio 

pl^4 FORMATlV// ' **** ERROR — ASSIGNMENT OF SURFACE MARKER PARTICLES 
1 Fyrr*FO« ?4 1 LIMIT. **** •) 

c **** A T SuRFACE marker particle will be ADDED OR ELIMINATED between *** 

C. TWO MFinH^OPtNO PAPTtOLPF. 

FNTRy OHKSMP 

COMMON/L4/ SPT<?)« DUMP (P. 60). NRWT3, CYC.LF ♦ IND^MP. LERROR » l SUR 
IMO c -MP = 0 
y? = o # 4 *DF 
y? = 1 ,**DF 

J « *> 

K = 7 

FPT ( 1 ) = FP( 1 ) 

FPT(O) = FP(?) 
no ">opo n = ?.NFP 
I F ( J .OT. 401) no TO 3 ? 0 O 
y 1 = CP < K ) - CPT(J-?) 

y 1 s CP( K4-1) - CPT(J-I) 

-r s enOT(V)*->? + Yl**?) 

!F(T #LT. V?) no TO 3 ?P 0 
IF (7 .OT, y? ) no TO 3??P 
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rpT(j) - f'P(k') 

r nT( J+l > " r P(K'+l ) 

J " J + P 
on To ippr. 

ippo J r < An r ( y t ) ,rr, AP^YII) 

M| r 1 

Y 1 r rPTOI I 4 P. n *Y1 
r,n to opof 

•ippr* Ml n O 

VI n rnT(JnP) 4 P.nflVl 

r **## in oA r ,r or MnNr>p« curve 
0 M-l AMP N instead or 

TpTO |.n K 

ir<N ,rn, mop) Lbl-P 
LI c. L - A 
LP a I - f* 

I 1 = L + P 
r i b cp< L 1 ) + oPtl o) + OP(L) + c.P(L’ , l 

rp = <5P<1 i)*#r» + OP ( I o)**P + on ( L. > * * P + e P<L7)**P 

CT = c .P(l1)**0 4- PP(L?)**7 + ^ILIW + e P(LT)**^ 

<■*4 b ep(|.| )*«A + SP(LP>**4 + PP(l.)**4 4 SP(l.*M**4 

r K = SP(L-O) + PP(L-l) 4 SP(L4l) + FP(L+?) 

ro B AP(L||I^IL-1| 4 «.P(LP)*SP(L-1 ) 4 OP t L ) *SP < L+ 1 ) 

1 + ^PdTK'APII. + T) 

C7 « pP< l-IHSPCL 1 )**? 4 SP< L- 1 > *SP < L2) **? 4 SP ( L+ 1) *SP ( L ) ** 2 
I 4 ^PIL + ^Ht^PILT)*^? 

or = o?**? - oi*o? 

on u ri^r^ - 4.o*CA 
rio - r 1 *o? — a • o» 0 T 
r 1 1 = r i #*o - a . o#o° 
n? B r?«rA - 01*07 
n 3 a. O?*03 - 01*04 
N? = *» 

JF<Ane.(OP) .LT • o.PPPOl) N?= 1 
1 F < Asc (Oil) .LT. 0.0000 1 ) N?=N?-1 
OO TO ( PPAO.TPAP^PPRO ) « N? 

COF-> e O 1 P/C. \ 7 
OOF 4 s <0« - 0 t P*C0F? )/C 1 1 
OO to 1 PFK 
-*?4P rnpi - ro/rio 

oof a = (01? - n won) /cp 
on to •a?c^ 

OPFO OOF 3 B (0PF0O - 011*01?)/<PP#C10 - 011*010) 

OOF4 s < OO - 0lO#0OFP)/01 1 

OOFo s f>i?*i*(OR — 01*OOF4 — Op*O.OFP) 

OO TO ( PPftO,TPPF) « N 1 

COFfr = OORT ( 1 .0 - 4.0*C0F3*(C0F5 - Y 1 ) / ( C0F4**? > > 

Yl a -O.O*0.OF4*( 1 «o - 0OF61/COFP 

00F7 = OOF'H + OOF4*\ 1 4 OOF7*Y1**P 

)F< Anc( Y1-COF7 ) .LT. 0.0001) OO TO PP 70 

VI a -O.P*OOF4* ( t .0 4 OOFF>/COFO 

OOF^ « O.OFF + '*Of'4*Xl 4 OOFP*X1**P 

I F ( Aoo ( y | -OOF 7 ) .LT. 0,0001) OO TO 3P70 

OPIKIT TPAP. OYOLF. K 

FORMAT < /// • **** FRROR - COMPUTATION OF XI FROM Yl tS INCORRECT. 

| • IPX. • OYOLF s* t^. PX« • K" =' 14) 

Pt‘TI inM 

OPRP Yl = OOF^ 4 OOF4*X 1 4 0OF7*Xl**P 


r,n to pppR 


rirriNG 10 marl my uisinc, n-t. m-p. 

N«P« M~l» N AMP N4I. 


*### 




••to 




A-27 


LMSC-HREC D225632 


O rpy(.j1 = VI 

r PT( .1+1 > « VI 

J - J + P 

I r < J .FT • Wl) cn TO ^OO 
no Tn TP in 
TpnP V n yr + P 
HP n J + 1 

iriMoPTT .Mr, 11 no TO r, r*on 

Ml n OKMPP 

pp | mt TPH/i ryri r.M r P» .1 

TpfVt FORMAT < ///•***■* VALUI r > OP ryci f « N r iF'» NP ( t r iP(ll 
PP I MT 7?W>. ( Ff'( ! 1 * I r \ .Ml 1 

tppp r opm at < // norn.pi) 

Pp I MT TP no . ( PPT ( 1 1 « 1 p 1 * NP 1 

ipon t i r n n (,!•♦! ) / f 1 
no Vfioq |n1«MP 
( M n ( ! ) 
nr t( ipm 


AHH OPTMP U » 0 ) 


OPIMT 

wri inM 


r 

r 

r 

r 


»«M« THF FF OP A FLUID IF PRESENTLY RESTRICTED TO A f ^^PLF ™NTCUW 
DFFINt-D NY NOT MOPE THAN a FCTNS WHICH AWE: MONOTON I C ALLY IN R JR Z. 

NO aiWVF CAN HAVE' 61 OW MORE' PTE, AND ASSIGNMENT OF PR MARK'FRF EOLLOWL 

Tor row patp of an ompty domain. 


DIMFMCION FMXI661* FMY<66)« C<°1 
vr s op(r>i - AO( | ) 
yr - J01- F»< t) 

jp(y^ # r-T. O .AND. ABF ( Yr ) »LT» YC) r,0 TO Pa 1 0 
|F(yr ,LT, n .AND. abc, ( yr ) .LT. Apc(vci) 00 TO 
JFIVr ,OT. O) no TO ta^O 


L CFO f 11 = 4 

no Tn 1/15F 


1 A 1 n 1 CFP< 11=1 
no To T/|p<=; 
-»alt- LFFO( 11=0 

ro TO - 5 A on 


-»4?o L.FF0(11 = P 

-»4P* iF0F(1) = 1 

M 1 - '»•* ( MFC) — A 1 — I 

NP = I FF0< 1 ) 

NP = 1 

KM= c 1 

no taoo k'='’.Ml»T 

M" 1 = N*» 4- 1 

yr u n P ( K ) - fp(f-t) 
vr = cD|y+n - np('. -i ) 

OO TO (TAAO.TAB0.TA60.TA70), NP 

7 = vr/yr 

|F(M 1 . OT . nil no TO TAA7 

IF (7 .IT. 1 .P 1 no TO TAOO 
rAU FOTM1 (W. 1 .P.L » 1 ) 

|F(L »ro. o) no TO Taor 
*M 4 t t cFn ( mc+ 1 ) P 
no To -lann 
•tflFO 7 a vr/vr 

tF(M-) .OT, nil no TO T4B7 
1P(7 .OT. -1.71 no TO 7400 


7 4 70 


TO 7400 

7400 


r , o TO -JAflr, 
GO TO 7404 


-j#KtCP - i 
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CALL rCTNI (K,"! .O.L.M 

f r ( u. ,ro, O) no to 7400 

740-7 |.<Tr,|M'-+| ) 
on to 7400 
nftf.n 7 n yr/yr 

IDMn , r,T » on on to 7407 
I r ( 7 ,it, um on to 'moo 
rAi i mini (K « i .P *1 • n 
1 1 (| ,ro, o) no to 7 /ion 
iao 7 i orn fM r '< I) n a 
rn to pa no 
7470 ? p vn/vn 

tr <m .oTt 01 ) on to 74 77 
tr<? ,ot. -!•;*) on to -vino 

TALI F0 7M1 <F,-1 .P7L.P1 

If (l ,ro, o) no to 'a oo 

7/, 77 I Of r ( MP4 1)^1 
76 fo HP ” NIC’ 4 1 

1 4 * ( Hr ,OT» r ) OO TO 7600 
tr(ir ,fOi Mil 00 TO 7606 
M f i = i p.rr, ( mo ) 

| err e- (MR 1 v w» 

M"* u 1 

7/,0n CONTINUE 

( c.fAP (MM+I 1 = 7*MCP - J 

rn to 74 ; | o 
7 /jO/i i orr(MO) - o 

l r-FCC(Nn) a 7*MSP - 1 
no to 7 t - 1 o 

7600 N J K rVCLT - 1 

74 c,* FORMAT^// / ERROR — PRESENTLY. FREF SURFACE IF LIMITED TO 

14 SFGMFNTS. SURFACE TENSION EFFECT IS NEGLECTED AFTFR CYCLE l » 
n • ( CTMT 740 c; > • « > 

I rponp e 1 

prTl |PN 

7c^in po 7ecn L=l»4 

t F ( L*FC ( L 1 .PO. O) GO VO 0*^0 
Ml = t°FGG(L ) 

IF (L .GT. 1) N1 =M 1 -1 n 
M? a I SFGS ( L+ 1 ) “ 0 

M a (M? - Nil /7 

no Iron 
j = (t + 11/7 

FMV( J> ts SP ( I 1 

eMv(j 1 - «P( ! + n 

N7 a l nr c ( L ) 

GO TO ( 74,70 , 7 c 70, reaper, 7C 70), N7 
7rapr, CALL ( FOPF 1 ( SYy , <-7MY , O « M » 7 ♦ C. « l FPR ) 
rn to 7 t 60 

iri n no 7n7n I sMt ,M7« 7 
ikik c,My( 1 ) - “ SMy< l ) 

C All L cTJ^F 1 ( SR"Y * GVy .O.M.7.C. IFPR ) 

7040 IFtlrPP .i'lr. n ) no TO 7 o<f 0 
no lead f=l.o 

nvr < t , 1 ) = c ( ! ) 

74,40 ini .no. P .OR. L *F0. 4) SMC ( L « I 1 = “C ( 1 ) 


7040 


4) SMC<L« I »*-C( I ) 


7 «no rONTT^ilin 
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PrTUPM 

N1 n OYCLF - 1 

PRINT 3S04«N1 . 

*,^64 FORMAT < /// • **** ERROR IN SUBROUTINE LSQPF1 (STMT 3340) • SURF AC 

IF TENSION EFFECT IS NFGLFCTED AFTFR CYCLE* 14) 

l.n’PPnP n l 

return 

rMn 

• FOR* IS cn^,,cnf 

SI (BPnt IT I NF ARROW ( I X 1 * I Y 1 « l Y? ♦ I Y? » I H I TE * I B ASE ) 

HP I O.hT b IH1TF 
BASF “ I BASF 

y? = ix? 

Y? = IY? 

nY = IY? - IYI 

OX s IX? - IX 1 

SO s 1 ,/SORT(DX*OX 4 DY*DY) 

FACTY = RASF * I4WDX) 

FACTy = RASF * <SO#OY) 

ys s x-> - height # <so*nx) 

YP = Y ? - HFIGHT * <SO*OY) 

1X4 = <Y3 4 FACTY) 4 .G 

TY4 s <Y3 - FACTY) 4 *G 

TXS s <XS - FACTY) + •« 

IYS 8 (Y3 4 FACTX ) + «S 

CALL LINFV (1X4* l Y4« IX?* IY?) 

CALL LINFV ( IXR* I YG* IX?« IY?) 

RETURN 

END 

• FOR* IS SO"? « S07 

SURPOIITTNF FCTNl (K*X*L*M) 

C0MMON/L3/ SP ( 48? ) * DS* STC* STH* STR* STZ ♦ STR?« STZ2* STRZ2 * 

1 STTRF* NSP* LSFG(S)« LSEGS(B) ♦ SMC<4*9) 


<SO*OY) 


( 1X4* IY4* IX?* IYP) 
( IXS* I YG* IX?« IY?) 


OA 1 = 1*4 
N4 4 ? 

( 3BI O * 3R 1 G ) « M 

r,P. j / tSP(N4 + l) - SP(N4-1 ) )/<SP*N4) - SP(N4“2)) 

I F ( 7 *LT. X) RETURN 
SO TA PR30 

TRIG 7 = (SP(N4) - SP(N4~?) )/(SP(N4+l )~SP(N4-1 ) ) 

IF (7 *OT. X) RETURN 
PP3C CONTINUE 
t. * 1 
PFTI ION 

FNTRV FCTN?(X*XTF»Y) 
y = eOPT(STR?-GTP7?*( Y-STC)**?) 

XI = P.g#XTF 
Y 1 = n,GOP1*XTF 

PP4 0 I F ( ApS ( X— V 1 ) *LT • Yl) GO TO 38G0 
Y 1 = Y1 4 XTE 

no TO PR40 

-»P^n y = Yl 
PETI IDM 
FNO 

•EOP*tS cno*aoo 

SUBROUTINE FORCF(K*Nl *N? *N3« N4 ♦ NS ♦ STU. STV ) 
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c ***# FORrF IS COMPUTFO ONLY WHEN SP(K) CROSSFS XC OR YC ONCF. 
njMFMSfON SMX (3), SMY (3) 

C0MM0N/L3X SP<482>. DS. STC, STH. STR. 5T7. 5TR2. STZ2. STRZ2* 
\ f'-jyPp « NSP* LGFG(*3) ♦ LSFGS(G) « SMC(4.R> 

NA n 0 


MO <7 o 

X - rl OAT(Nl ) - 1 
y s rl.f'AT(M'») - 1 
HMy w V - 0,F 
MPV If Y + n.K 
HVIy s Y — o • s 
HPY b Y + n «F 

l F ( ► , GT • LSFOS(N7+1 > ) N7=N3+1 

NR s t. SFG(N7) 

GO TO ( 371 0 . 3770 . 3710 . 3730 ) « NR 

** 7 io smx ( 1 ) = y - o.s 

SMy ( p ) = X 

s y + o.s 
no 371 ^ MbI ,3 

SMy(m) s SMC(NR.l) 
no 3**1 F N=1 .7 

37 IS SMY(M) = SMY(M) + SMC (NR . N+ 1 ) *SMX ( M ) **N 

! F < SMY ( 2 ) .LT. HPY .AND. SMY(2) *GT. Y) NS=N? 
j f ( SMY ( P ) »GT. HMY • AND* SMY(2) *LF. Y) N5*N2-1 
!F(SMY(1> .LT. Y .AND. SMY ( 2 ) .GT. Y .AND. SMY(3> 

I F ( SMY ( 1 ) .GT. Y .AND. SMY ( 2 ) .LT. Y .AND. SMY( 3 ) 

I F < SMY( 1 ) .GT. Y .AND. SMY( 2 ) * GT • Y .AND. SMY( 3 > 

I F ( SMY ( 1 ) .LT. Y .AND. SMY ( 2 ) *LT. Y .AND. SMY( 3 ) 

TF(Na .FO. 0 .AND. NS *F 0 . 0 ) RFTURN 
XC = (CMY( 3 ) - SMY ( 1 ) ) / < SMX ( 3 1 “ SMX( 1)1 

VC = SMY ( 1 ) - XC*SMX( 1 ) 

YH s Y + X/XC 
GO TO 7740 

7730 SMY ( 1 ) = Y - O.S 

SMy ( <? ) s Y 

smy< 7 > = y + o,s 

no 7-?7F Ms 1 4 7 
SMy (Ml B SMC ( NR . 1 ) 

nn 777 F n=i .7 

773 R SMy(M) a SMX(M) + SMC (NR »N+ 1 ) *SMY( M ) #*N 

I F ( SMX ( P ) .LT. HPX .AND. SMX< 2 ) *GT • X) N 4 =N 1 
IF(SMX< 2 ) .GT. HMX .AND. SMX( 2 ) .LF. X) N 4 »N 1-1 
IF(SMX(n .LT. X .AND. SMX(?) .GT. X .AND. SMX( 3 ) 

| F ( SMX < 1 ) .GT. X .AND. SMX ( ? ) .LT. X .AND. SMX ( 3 ) 

| F ( SMV ( 1 ) .GT. X .AND. SMX ( P ) .GT. X .AND. SMX( 3 ) 

t F ( SMX ( 1 ) .LT. X .AND. SMX ( ? ) • LT. X .AND. SMX( 3 ) 

!F(N 4 .FO. O .ANO. NS .TO. O) PFTURN 

ye s (evyn) - smx ( t))/( smy ( 3 ) - smy<iii 

yc B SMXM) “ XC*SM.Y( 1 ) 
yn = X + V/XC 

7740 y a yr*(vn - vo/u.o + xc**?) 

y s ewr (MO,?) 

7 s 0,0 

no j *>fo n= 1 .s 

y = v + FL 0 AT(N+ 1 >*SMC(NR.N+?>*X**N 
07 SO 7 = 7 + FLOAT (N* ( N+ 1 ) )*SMC < NR «N + ?) *X** ( N—l 1 
yr-MPo - (1.0 + y**7)**2 

GO TO ( 37*0,7770, 77F.0 . 7770) . NR 


• GT. Y) N4 ®N 1 — 1 
•LT. Y) N4»N1-1 
•LT. Y) N4*N 1 
•GT. Y) N4=N1 


•GT. X) N5*N2-1 
.LT. X) N5=N2-1 
•LT. X) N5*N? 
•GT. X) NS»N? 
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77 ftn t^mpi 3 y*n,n + y*#P> 

‘’TV r 7/TFMPP + Y/TFMPJ 

«*,TU n -Y#7/TFMPP - Y**P/TFMP1 

PFTtiPN 

7770 TFMPl a ( YO - X/XC ) * < 1 . 0 + Y #*P) 
<=iTV p -y#7/TFMPP + Y/TFMPt 
FTIJ is 7/TFMPP - 1 • 0/TFMP J 
PFTtJPM 
FNO 

•MAP«!L MA|N«MAIN 
L | R c,ye;<r,*MCFr:*. 
tm <;oi 

• rOPOtJT TPF«* sLHMACP • 

• PFW | NO * I | HMACP • 

• FRFF LHMArP. 

• XOT MAIN 

**** DATA RFCX **** 

C 

• FIN 
♦FIN 
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Appendix B 


LISTING OF LHMAC 3 PROGRAM 


(: 


, ARRAYS! 5000 > « 

( ARRAY? *U), ( ARRAY 3 1 V ) ♦ 


SO 




I 

I 

I 

I 

I 

( 

I 

1 
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Listing of LHMAC3 Program 

M sG N^LOC k-HFEO-HUNTSV I LLF 3D MAC PROGRAM (TAPF GFN 0Q0272) 

»ASG,T LHMAC7*T*GAVF06 • LOCKHFFDMAC3PROGRAM 
• nrw | NO LHMAC3, 

• POP * I « ,en) 

C LOCKHEED/HUNTSVILLE 3D MAC PROGRAM ( LHMAC3 * G4K CORE SPACE ♦ 4 DRUMS) 

r COMMON/L l / t T YPF * I RAP* JBAR * KRAR* IP?’ JP?’ ^P? * l ARP AY ! 5000 ) * 

1 AfrpAYl (5000) * ARRAvp(DOOO) * ARRAY3 C 5000 ) 

0 ARRAYS ( coon ) 

FOUlvALENCF ( I ARRAY * I FLAG ) » < ARRAY 1 «PHI ) » 

1 ( ARPAY4*1»') » ! ARRAYS, 0) 

1 MTFAFR TYPF ( ?? ) 

OATA tyPF/1?h P TN PAPPR * ?0#6H / 

FORMAT! 1 61 c ) 

CALL ! DFNT ! O * T YPF ) 

pp Ap SO* I TYPF * | BAR * JRAR* KBAP 

J P? s T BAR + ? 

JPP s JBAP + P 
FP? s YRAR + P 
! F ( T TYPB ) 7S,6S*6S 

CALL MAIN! T^LAG^PHI *U*V*H*n) 

GO TP ss 
CALL FMOJOB 

C.JOO 

FNO 

• pop « is c A9,<n? 

SUBROUT INF MAIN! I FLAG* PH I »U*V*W*D) 

C0MM0M/L1/ ITYPE* IBAR* JBAR* KBAK* IP2’ 

ARPAYl ISOOO) * ARRAY? ( 5000 ) * ARRA Y3 ( 5000 ) ♦ ARRAY4< *000 ) ♦ 

ARRAY >1 BO 0^ 100)< J7(l00)» XI?COO), Y(?OCO)* 7C?000)* PLTX« 

MARGNX* MARGN7* NSEGPT* NIMRKR* COSPSI * 

C0NTR3* CONTR4* XV ( 3 ♦ 3 1) ♦ ZV!3*31) 

ISH* IFUL* IFH* IB* ICST1* ICST3* 

MLMT 

DYODX* DYOD? 

I FLAG ( KP? » JP? ♦ IP?) * 


65 


7R 


1 


COMMON/L?/ 

1 PLTY* rppul* 

? CONTR 1 « CONTR2 * 


IS 


JP?* KP2’ I ARRAY ( 5000 ) * 


SINPS1 « 


IP1 


COMMON/L 3/ IF* IFH* 

1 jni* Mot * ICYCI.f* ICfll!6), 

COMMON/L 4/ DXODY* 0X0D7* D70DX* D70DY* 

C "tmn/’l. ' PHI CKP?* JP?* IP?) * U(KP2* JP?« IP?’ 


DIMFmsION 


** 


0(KP?* JP?* IP2) 


1 V( K'P?* JP2* IP?») * W ( KP? * JP? « I P2 ) * 

DIMENSION TITLFC1?)* I0PTC16). IPLTI16)* ^ M6) » GRTO 6 ) * 

1 GRy ( 70 ) ♦ GRY ( 30 ) ♦ GP7<30>« BDRY( 6) » NSFGV ( 3 ) * JPLANF ( .. ) 

DATA U* IFH* IS* ISH, IFUL* IFH, I B/ 1 000 * 2000 , 3P00 * 4000 *5000 ,6000 * 

1 7OP0/, ICST1 * ICST2* ICST3, ICST4* ICST5* I CST6/ 1 0 * I CO * 1 000 ♦ 2 ..00 * 

? 1 noon * GPOn/, CST1 *CST?,C«T3/0. 0,0. 0,0.0/, 

** P I .HP 1/3.141 5P?7* 1 • 047?/* ISIGN/-1/ 


)00 
1 1 o 

1 1 4 
1 70 
1 40 
1 CO 
1 6* 
1 60 
170 
1 74 
) an 
1 00 


popmat( PAA) 

FORMAT ( 161«) 
FOR»^ ATI 1 FF«, I ) 
POPMATI p R 
POPVATI 1 H 
FORMAT! // 
FOP M AT< / 
FORMAT! // 
FOPMAT! // 
FORMAT!// 

1 

FORMAT! /) 


pp 1 

n 

.4) 


1 HI 

% 

1 ?A6) 


// 

( 

1610) ) 


/ ( 

t 

6FB • ? ) 

) 

// 

( 

?6! p ) ) 


// 

( 

1 op 13, 

6) ) 

// 

( 

RF16. 

R) ) 

-M* 

« 

no. 

CP?0 


.P ) 
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Ml H r 
Ml IP - ° 1 

Ml 17 - 7 ? 

Ml I/I r> nn 
pritijMr i 
prAn ,pp. 


Tjxi F 


PFAd 1 1 0 « LNTH1* LNTH7, LNTH3, LNTH4 » LNTHB, LNTH6, LNTM7 
RF An 110, NMPPUX, NMPPUY, NMPPU* 

9FA0 110 , ( 10PTU) , 1*1 » lft) ♦ (IPLT! I) , 1*1 ,16) , ( IPRT < I ) , I ® 1 , 16 ) 
RFAD 110, NGRT, LHT « NVPLT, (NSFGV ( I ) , 1 = 1 *3 1 » ( JPLANF ( I ) , I = 1 * 3 > 

no pcs] J =1 ,MVPLT 
v- s MFPOVtJ) 


p=An 114, (XV(L» T ) » I = 1 ♦ <) 

PFAn 114, (7V(L, 1) , 1=1 «K) 

XV ( L , K + 1 ) = XV (l ■ , 1 ) 

pcs i 7VR,K+1) a 7 V ( L , 1 ) 

MVpDr a NMPPUX*MMPPUV*NMPPU7 
Jl a MGRT 4- 1 
JP a ?*NGRT 

RFAn l?0, ( RORY tl)»T**l,6) 

PFAP 1^0, ( OPT (t)»!*l,Jl) 

PfrAn 170, ( GRX ( I ) * I = I , JP ) 

PFAP 170, (GRY! 1 ) , t=l , J?> 

prAn 1^0, (FP7! I ) , I =1 » J3) __ nuA _ K11 . 

PFAP 170, Dr, DBFTA, DX, 0Y, 07, FPPA, FPGD, FPSP, FPSV, PHO, RNU, 

1 \/ e <“ALF, V'ALL 

RFAn 170, TtN, TPLT, TPRT , TCOMP, TFIN, TCPU 
I F ( I PRT ( 1 ) ,FO. O) GO TO 701 
PRINT 140, T I T|_F 

CALL FI_KOl.IT ( TCPU ) „ . . . 

PRINT 710, ( I OPT ( I) , I *1 , 16) » < IPLT( I ) , 1 = 1 , 16) » ( I PRT (1), 1*1, 16) 
FORMAT!// • VALUFB OF IOPT(I), IPLT(l) AND 1PRTCI)'// (1618)) 

PRINT 714, I TYPE , IBAP, JBAR, KRAR, LNTH1, LNTH? , LNTH7, LNTH4, 

1 LNTHF, LNTH6, LNTH7, NMPPUX, NMPPUY, NMPPUZ , NMPPC 
714 FORMAT!// * ITYPF at 17, 4X, • l PAR = • 17, AX, • JBAR = » 17, 4X, 

1 , KpAR = • 17, 4X , *LNTH1 = . 17, • LNTH? = • I?, • LNTH7 = • 17, 

i LNTH7 a * 


LNTHA a» I**, * 

t NMPPl )X a t 17 , 


LNTHF at 13 , t 

, NMPPUY a , 


LNTH 6 a* 17 // 


NMPPU7 = * 


t i , MMPPF = t IT) 

PRINT 31B, NGRT, LHT, NVPLT , ( NSEGV (t),I“l,7)» ( JPLANfT (I), 1*1, 3) 

718 FORMAT!// • NGRT =• 13, 4X, 'LHT =• 17, f X , .NVPLT =• 17, 4X, 

t iM^rOV(t) « 1 I 3 * 7X« • JPLANF ( 1 ) 313) 

7?P FORMAT?//’ • ( VALUFs' OF ’ GRT ( 1) « GRX(J), GRY(J), GRZ ( J > AND BPRY(I). 
1 // ( 1 ff 1 7 # A ) ) 

PPJM-r 170, ( GPX ( T ) , t = 1 , > 

pp I N'T 170, ( OPV ( 1 ) , t = 1 , J? ) 

PRINT 170, (GR?( I ) , I * 1 , JP) 

DR | NT 170, (PORY! 1 ) , 1 = 1 ,6) 

PRINT 7?6, OT, DRFTA, DX, DY, D7, FPSA, FPGD, FPGP, FPGV, RHO, 


DX a. F7.4, BX« 


1 PMm, V'-FALF, 'HALL . . _ v 

-?F FORMAT!// • DT = • F10.7, • DBFTA =. F7.3, ♦ DX =• , 

, ,nv =. F 7 .A, Fy, toz a. F 7 . 4 , «X, ' FPS A a. * 6 . 4 , AX, T.PBD =‘ 

? F 6 , 7 , 4X, • FP8P =' F 6 .A// » FP5V *• F K ,?, OX, » RHO =. FB,.-« 3x , 

ir 10.-7. • l/^PALF a, F 6 , 3 , • 'HALL at F4*l) 


, toy a, F7.4, ‘X , »n/ =’ ^ - 

? F 6 . 7 , 4 X» t FP 8 P =• F 6 .A// » FP 5 V =* F K ,?, OX, » RHO =• F 8 ..-, 3 .,, 

O ,r MU =, F 1 o , 7 , • VFCALF =• F 6 . 3 , • 'HALL =» F 4 . 1 ) 

PR I NT Tin, TIM, TPLT. T PPT, TC.OMP, TF IN, TCPU 

FORMAT!// • TIN =* F 6 . 3 , BX, 'TPLT a, F 6 . 3 , 4 X, 'TPRT =» FB. i, 2 X , 
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j * TrOMP »t FF.F, ' TFfM "• FH.A, • TCPU ~ • F«»l/// 

o « v/AlllFC OF y\/(J*|) AMO 7V(U»t>*/> 
no Onn ,j K 1 , MV/PI T 
V r MfPGVCJ) + 1 
I " J 

npjMT 1^4* (VV(|_« ! ) « Ini *F> 

Tin PRINT 1 r >4 , ( 7V < L « I ) « I " 1 ^ > 


lO| rOMTfMIIP 

J1 a IPP«JPF#FP* 

|P(J1 .ip. irsTM oo to 4 on 


PQ 1 NT 400 , 
400 FORMAT!// 

1 AT I ON <400 J 


jpp, IP?, KHp, .11 

#* rRROR "* EXECUTION IS TERM I NATFD DUE TO SIZE 
VALUES OF IP?« JPF , KPF AND THF I R PRODUCT ARF • 


LIMIT 
3 l 4 « 


o l a, • • • ) 


prri tr>M 


AOe: PONT I M if 

rorpp T s roc (PS I ) 

S IMPS I = St ►KPS I) 
tpmd i s | NTH 1 
TPMP? = L.NTHF 
TFMPi s I.NTHF 

OO TO F!l ««5»1 ♦ «*! ,npi ) ; ITYPP 
Nl J1 = TFMP1 + TEMP?*COSPSt + 5*0 
jp - ypMpo + t^mpf^SINPSI + 1*0 


on to R71 

o? 1 CONTTMUr 

nEJ NIX. .F.6.4 

«=fS4 FORMAT (// ' ** ERROR #* 


THE MARGINS OF A PLOT HA? NOT P.EEN DEFINE 


I n vf t ■ 071 ) • * ) 

PFTI ION 

C * CO WP RctR PTS FOR PLTS* * ASO MPKR PARTICLES* 
«=,7I NRPPllL = 1 0F3/MAV0 ( J1 « JF) 

RPPUI = NRPPU- 

marcmv * < i of? ji*nqppul>/p 

MARON'7 * ( 1 OFF — JP*NOPPUL>/F 

OO TO (nRl ,701 ,701 ,7 n l ) , ITYPF 

*FRl J1 = tpmpf*cospsi*pppul 
j o _ TFMr?* e i mps I #o°pul 

JF a LNTH1 4MPPHUI. 

JA a l KITHF*NRPP<IL 
VPMOV a FLOAT( J1 )/F»0 
V/PMC.t a FLOAT ( JF I/F.O 
PO Foe I a 1 ,NVPLT 
Y a K|CFOV<!) + 1 


PO epe J 

= 1 ,K 


y v < t , J i 

- pppi l| X y\/ ( I , J ) 

+ vpmoy 

7V < I * J ) 
M = 1 

a RPPl IL*7V ( I « J ) 

+ VPMG7 

no sol I 

= 1 ,? 


! V ( N > = 

MAPOMY 


1 Y ( N+ 1 ) 

a MARONX + JF 


lY (N+F) 

a 1 V ( N+ 1 ) + J1 


1 Y ( N+O ) 

a I V ( N + F ) - J*» 


1 Y( M+4 ) 

a l V ( N ) 


T7( l ) = 

V A ROM? 


17( !+*) 

- MARON7 + JF 


17(1 +n ) 

a MAP0M7 + J4 


t 7 ( 1+7) 

a MARGM7 + JA + 

JF 
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nOj M r? M 4- ^ 

f<* r: n 

M n r 

nn a^i I nn , i ^ 

T V ( 1 1 n fV(K + ?l 
1V( n 1 V</+J ) 

1 V( 1 +A 1 r I V ( / ) 
l?(M> * T7<M~/|> 

!?( I 1 « I7(Z+? ) 

17C I+O ) « 17(K+1 ) 

I7II+/I) o I 7 < z ) 

/ r> ^ - R 
r,n i m n m + o 

NOFOPT a 1 ft 
70 1 FONT t Ml IF 

TF(lPPT<?) ,E0, n ) GO TO 711 
PRINT 70/1, NGFGPT 

7 O 4 FORMAT!// • NSEGPT = ' I3» • AND VALUFS OF I X < I > AND IZ(I) ARE GIV 

1 FM nrt.Oui.M 

PRINT 1ftO« ( !X( I ) « I»1 «NB C ‘GPT) 

PRINT lftO« ( ! 7 ( I ) f 1 = 1 ♦NFFGPT) 

711 CONTINUE 

IP1 = TRAP + 1 
JP1 = JRAP + 1 
/PI = /(TAP + 1 
CONTP 1 a I BAR 
CON TP 7 = JO AP 
C0NTPO = /BAP 
CONTRA * 0,0 
PO 7A 1 /si »ZP? 

OO 7 A 1 J =1 % JP? 

no 7/> 1 I = 1 ♦ IP? 

I FL Ao ( K , J , I ) a IF 
PH I < K « J « t > = 0 • O 

U(K«J«|) = 0,0 

V(K«J«U a 0,0 
'«!(/, J« T ) a r '*° 

o</« J« n = °. n 

TFMPl s 1 . P/FLOAT (NMPPU* > 

TFMR-> = 1 . P/FLOAT (MMPPUV) 

TrMPs *t 1 ,o/n.OAT(NMPPU7) 

yr a tOAR 

YC = JBAP 

7C a LHT 

yy = o,r,# T FMPt 

YY = ",K#TFWP? 

77 e o.r;#TFWPo 
PLTY a TF^P? 

MMAP/D = 0 
M s A 

GO TO <701 , 1 IOI , 1 IOI , 1 101 ) « t TYPF 
70 1 TFMP/i = 7C - TE^PO 
PL.Ty a yc - TEMPI 
7qc. k' s • o 4* 7*7 

7A1 J a + YY 

7*e; t - + yv 

MMAPkP a NMAPKP 4 i 

M - M + 1 
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X( M 1 n yy 
Y(M) n YY 
7(M) n 77 

IF(M ,Nr, J F F T A 1 GO TO 771 

WRI Tr INI I) ) ( X ( L) *UitJ « | CSTA 1 # (Y(U) « l- n 1 ♦ I OFT A 1 * (?<U) iU n t * I CRT A ) 

M r fS 

77| rONTIMUG 

| FLAT (**« J« ! ) n in It f I 

yy n yy * tgmpi 

lF<yy • I - t • y o go to 7ft n 

yy <7 n,G*TPMP1 

YY n YY + TFMPP 

IF(Yv «LT. YG ) GO TO 7ft 1 

yy o,q# T pvipj 

VY = o,«#tfmp? 

77 e 77 + TFWP7 

IF(77 .LT. TFVIP4 .OR* 77 .GT * 70 GO TO 770 

M | MPkO b NM/VPIFO 

yy r n,q*yy 
YY a '',A#VY 
TFMPl a n„F*TFMPl 
TFMPo a n,c*TPMP? 

77*=: IF (77 « L T » 70) GO TO 700 
|F(M .LT. 1 ) GO TO 701 

WttlTF(NUl) (X(L).Ls|«M)« ( Y ( L ) ♦ La 1 . M ) . (Z(L).LaliM) 

701 CONTI MUG 

PO 700 JaP>,JP1 

no 7 pf i=p, i p i 
TFLAr(t.J.T) a to 

700 I FLAP ( yp? « J . I ) = !« 

PO 701 k'aF.KPI 

no 701 JaF.lPI 
IFLAP(y.l.T) = in 

701 I FLAP (K» JP?t 1 ) = in 

no 7oo K = p,K'P1 

00 7on JaF.JPJ 
IFLAo(y.j.l) a tn 
7 HP tFLAO(y . J« TP^) a in 
11 ° 1 FONT I Ml IP 

nxonv = nx/py 
nxon7 a nx/n7 
n7onv a n7/ny 
n7onv a 07/PY 
OYonx = nY/nx 
nvon7 a nY/n7 
nxnx a nx*nx 
OYny a oy*oy 

0707 a 07*07 
nynv a ny*ny 
0YO7 = 0Y*07 
O70X = n7*ny 
cnri e o,o*DMi i/ny 
FOFP a P . n*PMI i/nv 
COFO a P,n*PNll/07 

COF/i a (1.0 4- FPFA )/( F.O/OXOX + ?• O/OYOY + 7.C/O70Z) 

ITPMP a |PPT(C) 

IGtMT a t HT + P 
KLIMT a KI.MT 
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prw|Mn mma 

U/R I Tr ( NIJA ) ( ( ( p ( K . J * I ) « I k 1 « I PP ) » J* 1 # JPP > • K~ | * l<|_ I M I T ) 

CALI. rri_|_ JP( | F|_ Aft « PH I oil. V.ul.D. I TFMR) 
f rvr | r = r 
T n n ,n 
JCMP n n 
JPI.T = n 
,JPPT n n 
J Mm i|i - 1 

1 OH i |-P( T n T/TPI.” 

FORT n T/TPPT 
k-rMP r T/TCOMP 

IPit'ni t ,mf, JPI.T) on to 1RP1 
I F < 1 PI T < 1 ) .FO. O) go TO iaoi 
17 n ITVRf 
1 ■> - 1 OPT ( 1 ) 

f A = 1 

prmiMO Mdt 
M n MMAPk'P 
CALI FPA»1FV/(T) 

1A01 1 F(M ,r-T t IG^TA) MsIGGTA 

pFAn(MM > (VtL ) »L*1 «V> ♦ (V(L) »L=1 «M) . ( ? (L ) «Le 1 .1* ) 

FALL PLOT < 1 . tTYPF. 10. I4.M) 

M s mMAPFP - 1 A*IFGTA 
1 F(M .IT. 1 > OO TO 1 A 0 1 
TA s 1 A + 1 

OO TO 1 A?t 
) A o ) FONT 1 Ml !«=■ 

IFItPLTtF) • FO . P .OR. KPLT .FO. 0) GO TO 1FR1 

k|Ct m t s l^?! 

OO le-’l L=1.MV/PLT 
FALI FPA**F\/(0) 

J = ,|0| ANF(L.) 

1A = Mcrrvd.l 

Man 

no IRIK f = P.F|.mt 

•77 S OPPllL* (FLOAT <F ) - l.*) + VPMG7 

no |e|| I =?. IP) 

!F( fF| AO(K, J, t ) .LT • JO) GO TO 1 R 11 
M B M + 1 

V ( M ; - ppm l| * { ri OAT ( | ) - 1.0) + VPM GX 

7 ( M ) s 77 

Y(ICCT0+M) = O . f=»*RPPUL*VSCALF# (U<K»J.I—1 ) + U(K.J.t)) 
7 ( 1 COTO+M ) = o . e^pppuL* VFC ALF* < V < K- 1 , J * I) + W(K.J.I)) 
i *■ \ 1 OOMT |MI ir 
1 k 1 f FONT TNI l 17 " 

jF(f/ .OT. tocto) GO TO inn 
FALL PLOT (0 . |fct7.L« IA»Y) 
inr>i GOMTtMllo 

|nq| JPLT a JPI.T + 1 
)ftn| rONTfMItr 

I F ( l OPT { 1 ) .FO. O) GO TO lOOl 
IFft'oriT .MR, JPPT) GO TO tool 
PFliltMO LM I A 

PF AO (^ I UA) ( ( < O ( F ♦ J« I ) ♦ I = 1 ♦ t PF ) ♦ J= 1 « JP 9 ) «^~ 1 « FL 1 ^ I T ) 
OO |7HK FsP.Ft.MT 

no t 7nn j=->, jot 
no 1701 l=p»!oi 
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f rr in. /\r,(K« J. n .1 t. IFUi ) r.n to I7r>i 

n<K, j, f > n orr. J* f ) 

i " <nni < k*j* i4 n + ph i ( k , j * i - n - p#o*phi cr« j* i n/nxnx 

p " ( PH |(K«J+|«|) + PH! <K. J~| * I > " P.O«PH! ) /OYDY 

7 - ( PH | ( K+ | « J « f ) + PHI (K- J * J* I ) - f'«n*Pll! |K» J» I I I /D7D7 

17ni font r Mir 

1 70 n rONTtMlir 

print ipio, t * cv * ov. r. 7 * irv<"ur’« »tpp 

jM|n P0RM/\T(IH1» • V/r t_OF I T I F A AMI t PUrTiAl iPf AT T I Ml ’ • rfUA* • PP-P* OX ni 
l po, 4 * * r.v n 1 ro,4« • 07 n • r 0 . 4 * * icyri.r rl !A* • I TRR r • ! r // 

P * u- j j in AO hi ov, *yi oy, him ov, »n* py, *PH|» oy, 

7 «| in AO III oy, »\m oy, » um ov , m» OV, »PH!*) 

00 1070 kp 1 ,y| MT 
pn 1070 jri,jn? 
on io7i id , |pp, 0 

1 c 1 + l 

IFd. .rr, mm no to ioai 
1 r (1 ,rn, 1 ) PPIMT 100 

PRINT 1070* K , J, I , !FLAG(K« J« I ) * U(K*J«1>* V(K*J*I>* 

1 DIk,Ji|)i PH I ( m * J * l ) » L » I FLAP! K * J*L ) * 

,7 W<y,j,L)« 0 ( K « J « L ) * PH 1 ( K * J *L ) 

1070 FORMAT! 7 17* I A* K F10.4* 71 A* pniO.A) 

no to 1071 

1061 PRINT 1070* K * J* I* 1FLAG ( K « J * 1 ) * U(K.J*I). V(K*J«1>. W ( K ♦ J * l ) * 

1 O ( v * J* 1 ) , PH 1 ( K « J« 1 ' 

1071 CONTI MUF 
1 P7 e FONT J Nt JF 

JPPT a JPDT 4- 1 
1 PO 1 FONT 1 Ml IP 

IFtKrVD .MO, JFMP) OO TO 1^)0] 

F7 a 0,0 
FV a 0,0 
F7 a 0,0 
RMy a 0,0 
RMV a 0,0 
RMy = 0,0 

tP = ? 

no to ( lPl 1 « 1R11 ♦ tool ♦ 1901 ) « ITYPF 
loll no looi yap ,kt mt 

TFMDi a K 

ToMpn a 07MTFMP0 - 1 ,d 

no 1o71 Ja?«JPl 
TFMF'7 e - JOAR/P 4 J 
TFMp* a nV*(TFMpo _ 1,et) 

TFMP 4 e RHO*DYD7#(PHl (X* J* IP1 ) - PHl(K*J*tP)) 
rv s Fy + TRMP4 
PMy a RMy + TF M P7*TFMP<a 
1071 PM7 - RM7 - TFMP7#TFMP4 

no io?R 1=P«1PI 

TFMp 1 a - 1 r AP/P + I 

TFMP 1 s OV# ( TFMP 1 - 1 .<i) 

TFMPr a RH0*07DX* (PH! (K - * JP 1 * I) - PH I ( K * 1 ? * I) ) 

FY = FY + TPMPA 
P My - PMy - TFMp7*TFMP e ^ 
in?q PM? rn PM7 + TFMPl •*T pr MPR 
1071 CONTI MUF 

no 1041 l a O « 1 P 1 
TRMPl = - IOAP/7 4- 1 
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tfmpi b nX*<TFMP 1 - i ) 
no idai 

TPMPo " " JPAR/f* + ' J 
tpmpo - nv#<TPMPp - i«ni 

TfMP/. r. PHO-cpyPV*(PHt (KP 1 *,.J. I ) PHI ( I f* # I > ) 

f? „ r "> + tfmpo 

PMy n PMV + TfMn f *HfWPA 

) Ort j PMy n PMV " TFMP | #TPM n A, 

nn i op \ ynf*«y| mt 
no inn i ,*F> r >«jP1 
no i no i l •-> P • I P I 

jo<~, 1 0 <K«,lt I ) o PHO#PH I 1 F « J * 1 ) 


I ^ ? 

ir> 7 r FORMAT!//// • PWFftWJMr 0 1 STM I BUT I ON < N/M* *P ) ON TH« ‘■EKTJK ♦ J* *>• W 
|GHT<y,jC|IMI. FRONT 1 K « i? ♦ I > « HACK < K « JP | » I ) ♦ HOT TOM ( g * J « I ) ANO TOP!*. 
!/ • WALLS OP THP CONTAINER, HeW-KCTIVELY * > 

no ! mt 1 op 

no lo7i Jo?* JP 1 


N ® J 

1071 PRINT \91P % (0(K « J«L) «K=?«KLMT) 

107*5 rORMftT( 10F1 7, fi) 

PRINT 1O0 
00 1077 Je P % JP I 
M = J 

1 P77 PRINT 197?, (D(K« Ji |P 1 ) «Ka?»KLMT) 

PRINT 1O0 
‘ io7A PsP^LMT 


9 V 

1« ^ -PINT 197?, (0!N,L , I 1 IP1 1 

print l on 

no 1070 KsP^LMT 

N a V 

JQ7K PRINT IOTP, !0<N, JP1 , I > ♦ T=7>« 1P1 ) 
PR! NT IPO 

no 1070 j=?»jpi 

N = J 

IP 7 F PRINT 1 07 ? , !D!L,N,1>«I=?*IP1> 

print ioo 

no 1077 J*?»JP1 

N 9 J 


J 077 PRINT 147?| 
inpi JCMP r JCMP 
PRINT 1994» 
1 ORA FORMAT!// • 
1 * ♦ F7 = • 

? r i o«4 , • 


iniFPi « n» n » i «?♦ ipi ) 

+ i 

7, fry, FY, F7, RMX, RMY , RM? 

T =• FlO* K » • ! SFC ) ♦ FX =• F10.4. 

F10.4, • !N>» RMX s» F10.4* •» RMY 

< N-M ) • • ) 


1 oo j rONT 1 Nl IF 

r #C.k IF nnPY prfsfurf nffoo to of aojustfo. 

OO ? 1 «« F=?«KLMT 
no P i J=?*JP1 
no p i«i i=P*l ra l 

iFMrl A 4 CK.J.I) .OT. IOH 1 OO TO ? 1 1 1 
IF( IrL AO!F« I 1 .LT. JO) OO TO ?141 
M = ( IFLAO(F» J« 1 ) “ l S ) / 1 CST 1 

1 P ( M ,OT« 10 .OT 7 ) W = V-I 00 T 7 
1 F(M *FO, 1 »OR • M «FO» A > OO TO Pill 
IF (VI # rCj* P #OR« M ,FO, fl ) OO TO ? 1 P 1 


FY = • F l 0»4 * 

FI 0 »4 ♦ H RMZ =* 
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IP CM ,FQ. IP) • OR • M • FQ « 3D GO TO PI 31 
GO TO PI A 1 

Pltl OH I < k* , J , | ) « rori *(UC !'♦ -j* I > " IJ < K , J , I'M 1 I 

GO TO MM 

PI PI PHI (K*J» I ) ~ rOFIIIHIKi U I I " WCK-I «J» I ) ) 

GO TO P1M 

p| 3 1 PHICK».J«M f rO(-p*CV(K« I I V(K*J"1«II) 

GO To DM 

P 1 A 1 PH I i r « . I , I 1 r n • n 
MM rOMTTMHf 
M r *n. PONT I Ml IF 

a p g f gm Pinn to mf 

P 1 M 1,1 T 1 

no Pi 71 *lo P , JP | 

no pi 71 I i'P « t f 1 1 

ni-t| <i i , J, I ) a phi <11 + 1 «J« 1 1 

PI v 1 PH! (FPP, J* I 1 PM! CKP1 • J« I 1 
no pi7G M=p»rLMT 
no p i 7 r loPtim 
PHI ( K ,L 1 , 1 ) = PHI PML1 + 1 * 1 1 
017c, PH t ( y* « JPP, I ) o PHI CM JP1 , I 1 
no pipi K 2 ?,^lmt 

no P1«1 JeP, jot 
PHI <ir, J.L1 1 e PH 1(K"»J«L1 + 1) 
p 1 O 1 PHI dM J, 1PM = PHI PM J« IP1 1 

GO To Vf 
Ojnc; CONTTMIIP 

ICVCI. F “ 1CVCLF 4- 1 
T e T 4- OT 

1F(T .OT, TMM) Oc-TUON 
nr i)i i mo mi ip 

WRl TF (NUP ) ( ( ( PH I ( K , j » 1 1 , I S 1 « IP?) * J= 1 * 1 « K::: 1 ,KLMT ) 

no p i o i 1=1, MOOT 
| F* ( T , r T , OPT ( 14-1)1 ro TO MM 
J = - 1 

TFMOl = <T - CRT < I ) 1/CGRTC 14-1 ) - GRT < I) ) 

OV s ORV(J) 4- TFMP1 # ( GPY< J4-1 1 - GRV (J1) 

GY = ORYtJ) + TFMP1*( GRY( J4-1) - GRYCJ11 

G7 = GP7CJ) 4- TFMP*. *CG07C J4-11 - GR7 ( J M 

ro To Moo 
MM CONTI Nl)F 
r> 1 on FONT 1 Mt in 

pn ^i-iT) Ha.MM MT 
no pooo jsMJPl 
no t =o , ipi 

1FC IPI . AOCK, J, I ) .LT. !M GO TO ??31 

■pr VIP \ e 1 1 ( < , J « I 1 « ( U ( K , J , t 4- 1 ) —U (K,J,I**l)1 /OX 

TPMPp = V(K, J« I ) * ( VCK» J+l , t )-V( K', J-l , 1) 1/DV 

TOMP-j e \»i ( K , J , 1 1 * < W ( K f 1 , J , I 1 ~VI ( k— 1 , J , l 1 ) /07 
T rMP/| = ( ( j ( K « J4 1 , 1 ) 4-1 M K , J , I ) 1 * ( V ( K , J , ! 4- 1 ) 4-V ( K , J , t ) ) /A • 0 

TOMPo = < V ( K 4- 1 , J ♦ 1 ) ♦V (K«J> 1))*<UI<K,J4>1,1 ) +[•! (K,J, 1 1 1/4,0 

TFMPr = (V'(K, J, 14-1 l+P'CK* J, I 1 )*<1'<K+1 1 )4-U<<»J» 1 ) 1/4,0 

f =Mp-f - (U(F , J, 1 1+lllK , J- 1 , I 1 1 * ( V< K » J— 1 , 1 4- 1 1 4-V < F , J— 1 ,111/4*0 
TFlPPn = ( VCK, J, 1 l+VCF. «J, 1-1 1 1*CU<K, J4-1 , 1-1 14-UCK, J, 1-1 1 1/4,0 

TPMPo = ( V C k; ♦ J , 1 ) 4-V < K- 1 , J , I 1 1 ( VM K- 1 , J4- 1 , I 1 4-W < F - 1 , J , 1 1 1 /A , 0 

TFMP t P= ( n»< K , J , t 1 C K , J- 1 , 1 1 1 * < V < K4- 1 , J- 1 , t 1 4-V ( K , J- 1 , II 1 /4 , 0 
Tf*M( ' t 1 = (WCK, J, t 1+WCK.J, 1-11 1*(UCK+1 , J, 1-1 1^U<F,J, I-I 1 1/4.0 
T rMP|*s (UCK, J, I l+UCK-l , J, 11 )*< WCK-1 ,J, 14-1 14-WCK-l , J, 11 1/4.0 


\ * 
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TrMPl T n RNU*( IIKK.J. l+n+UIK^M"! * * > /DXDX 

1 + ju(K . J+l « l > +1 1 f K « J-1 i 1 )~7.0#U(K» J » I ) l/BYOY 

0 + iijiKH.J.IItlllK-I.J.II-^IlK.J.mAVn?) + C, x 
irMi>i4 r rnh# < <\/(K« J« h i ) +\/< i* « j . i-i i — p • o-w v < k . J« I > )/nxr>x 

1 f cvtr, J+ 1 , d+'/ip.i-i . n-r.n«\/(K, j, n i/pvnv 

•> f ( w(|« ' + 1 . J. I HVIK-I •-!« 1 )-7.0* Y< K »J. l ) )/(V07 ) + GY 

‘Trip n UNl.«M W (^J.lH)^(K,J,l-n-P.O«W»K,J,!))/nxnx 

1 f (W(l- »J+l ♦ I 1+WIK « J “1 * 1 )-7.0*\i'<K. J » I > l/HYOV 

r» + ,viini , j. |)+W(X -1 ,j. i >-o.o#w(K, j, n )/pwi + C7 

fl'l AO, ( k'.f 1 . J« ! ) . r Q. IP) GO TO 2r 1 1 
+ OT-MTFMPIO + ( PH l(K.J.l) - PH I ( K + 1 , J ♦ I ) ) /OZ 

- T f-NIPB)/t)Y + (TFMP11 - Tt-MPft ) /OX ) 

1 FLAG < K « J+ 1 ♦ 1 ) .F O. IP) GO TO 2221 
+ DT#(TFMPM + ( PH I ( K » J ♦ I ) " PH I ( K « J + 1 ♦ 1 > ) /BY 
_ j|-Mpc;)/D7 + ( TFMPfl - TFYP4)/OX) 

I FL AG < K « J ♦ I + 1 ) .FO. 1°) r ’0 T0 2231 
+ DT*(TEM p l^ + ( PH I I K « J » l ) “ PH I ( X . J « 1 + 1 > ) /DX 

- T cr^p5)/-)7 + ( TFMP7 - TFMP4 ) /DY ) 


II* (K .ro. 
W(KiJ*I ) ! 

j _ T rMP?. 
7211 IFIJ .FO. 
wiKtJtl ) 

J - T FMD? 

7221 IF (I .FO. 


KP 1 .AMD. 
; u,l ( K . J . I ) 
+ (TF^PIO 
JP1 .AND. 
s V(K. J* I ) 
+ ( tempo 
I PI .AND. 


U < K « J ♦ I > = U ( K » J . I ) 
j - t^MPI + (TFM017 

277 1 CONT TNI IF 
7777 CONTI NUF 

A^ e TON 77R1 TO Mr 
*> 071 = CALL VFLFTYI IP’LAO.PHI 
1.1 = 1 
M = 1 

no 7041 J=7»J p 1 

?F( lFLAG(Lll!u l ) .FQ. IP) U<L1.J.I)=B0RY<N)*U<L1+1.J.I) 
IFUFl AG(Ll.J.I) .FQ. IP) V(Ll.J.n=BDRY(N)*V(Ll + l.J.I) 

IF ( IFL AG I XP7 . J » I ) .FO. IP) \/(XP 7 , J. 1 )=nORY(N+1 )*V<KP1 «J. I ) 

.FO. IB) UCXP2.J. I )=BDRY(N+1 >*U(KP1. J. I ) 


774 1 


77*1 


.FO. 

• FO. 
.FO. 
• FO. 


IB) VI X. J.L1 ) *FORY( P) *V IK. J .L 1 +1 ) 

IB) V‘M X » J ♦ L 1 ) nBDRY I N ) *W I X . J . L 1 + 1 ) 

IB) v ( K « J . I P2 ) =BDRY ( N+ 1 ) #V ( K « J . I P 1 ) 
IB) W < K . J « I P2 ) =BDRY ( N+ 1 ) *W I K « J . I P 1 ) 


77B I 


770 1 


• FQ. 

• FQ. 

• FO • 

• ^O# 


IB) U(K«Ll«!) =BDRYIN)*U<K.Ll+l . I ) 

IB) UMK.Ll ♦ 1 )=HDRYIN)*WIK.L 1+1 . I ) 

IB) U ( K « JP? . I ) =HPRY ( N+ 1 ) *U I X ♦ JP 1 « I ) 
IB) w ( K « JP2 . I ) =BDRY ( N+ 1 ) *W C K . JP 1 » I > 


IF < IFLAGIXP7. J, I ) 
m = m + 7 

00 77«1 Ks 7 .Kl.MT 
no 77*1 J= 7 .J n l 
IF ( I*LAGtK. J.L 1 ) 

IF ( I PL AG (X . J.L 1 ) 

I F ( IFLAGIK.J. I P 2 ) 

I F ( IFLAGIK. J. I P 7 ) 

K 1 x m + 7 
no 7 *>^i k= 7 «klmt 
no 1 7 *i i = ?« ip 1 
IF ( in. AT (X .LI » I ) 

IF< IrLAGIK.Ll . I ) 

I F < I Fl. AG< K . JP 7 « 1 ) 

IF ( IFLAGIK. JP 7 . I ) 
no TO Mr 

no 7 *>o=; k = 7 ,kL m T 
on ?->nc J= 7..)01 

7 r t t'Zr AC ! K « j!*i ) .LT. IFUL .OR. IFLAGIK.J. I) *GT. IFH) GO TO 72 p l 

nixlj'n = culx.^n 1 uik. j. ,- i ) )/ox ♦ (vik.j.d - vix.j-i.m/DY 
^ 1 -pMP j ^'l^u'l^K I J . I — 7 )^U(X«J«t"*l ) +*J 1 X ♦ J « I )*( U< K.J. 1 + 1 ’ “J* 5 K ^ J- 1 T I )) 

T r,v,pp = V(X.J- 7 .I)*V(K.J- 1 .I)+V(X,J.I )*<V(K.J+ 1 . - 2 . 0 *V X.J 1 . 

TrMP-, . W< K-?.J. I )*WIX-t.J .1 I+WIK.J. I )*( WIX+l. J. !)“?• 

TFMP 4 - (UIK.J.I ) + UIX.J+ 1 .I))*(V(K.J.I ) + VIX.J.I + 1 ) 

1 TtulK.J-l.l-l) + U(K.J.I-l) )»(V(X.J- 1 . 1 - 1 )+V<X,J-l . 1 ) ) 

-> - M ) ( K « J— 1 . 1 ) + UIX.J. I ) )* (VIX. J-l « I ) + V(X, J-t « 1 + 1 ) ) 

; - ;; )( K.j,t-n t uiK.jM.i-m^vix.j.i-n + V(X.J.D) 
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+ U(K+|«J*ll)*IW(K»J*n + W < K , J « I + m 
+ 1 1 < K * J, |"1 ) )*(W{K~1 . j» t~l )+W(K-l « J« ! ) ) 

U(K , J. ! ) )* (UK t'-l « I ) v VCK-1 • J» l + l ) > 
tMK + 1 « J« |-t ) ) *(W(K, J. t-1) + W ( K « J » I ) ) 

t v(K+i » j« n ) + < w<k,j, i i f w(K*J+t»in 

+ V(K« J-l ♦ I ) HM W(K-| « J-1 « I ) + '*MK~1 « J« | ) ) 

\H Y «J« t ) ) + <U/( K~1 • J * I ) + Ui(K" J « J + t « I > 1 
\/(K-+1 « J-l ♦!))*((•'( K'* J-t ♦ H + UKK.J.!)) 
j, n/HT - TFMPl/nyDV - T -IPP/DYD Y - TFMP3/DZD7 
TrMpp</n7Py + tfmpa/pyoz ) 


jpMPr r (U(K » J« I ) 

1 + ( ll( K- l • J« t ~ l > 

n - (ii(K-i « j« i ) + 

3 - ( 1 1 1 k « j , i - n + 

TPMP* n (\/(P » J« 1 ) 

1 + ••-J’ - '1 •! ) 

T - (\/(K-l , J, ! ) + 

•» " ( \/ ( K* ♦ J" 1 • t > + 

ph I ( w « j « n r n ( i< # 

1 - n,n#{ TFMPA/OyOY 4 

Tin) rOMT \ t- Jl If” 
noon OOHTfMUF 

no 

no J=P,JP1 

?FC IFLAgUCjIh .LT. IFUL .oh. |FLAG(K,J*n «GT. IFH) GO TO 2301 

PH I ( * t J • I) = PHI (K» J. 1 ) + RNU*«(OCK,J,! + n + n <^ 1 " 1 'i^vl 
i k.j.i n/nxox + <n<K t j+i«n + n<K*j~i,n - p.o*o(k, j. n >/dyoy 

p (p(K+l*J«l) + n(K-l « J. I ) ~ 2* 0*D < K . J » 1 ) ) /DZPZ ) 

^5 rOMTINUF 
?nnr rONTTMtjF 

no V - P « Yt. MT 

no " Je’.JPI 

no poll l «? « T n 1 

ncF.j*n = . n .° 

1F( lFLAG(K« J« I ) .LT. IFUL .OR. 
n ( K , J . !) = PH ! ( X % J * I ) 
o'* 11 fONTTMUF 
PP in CONTTMUF 

K*L 1 Y T T = K'L’VT 

pnniwn NU4 

WMITF (NU4) ( ( <D(K, J. I ) « 1=1 ♦ IP?) « J =1 »JP?> ' K=1 .KLIMIT) 

RpAp(MUP) ( < ( p HI (K, J. I ) . 1 = 1 « IP?) « J =1 • JP2) *K=1 « KLMT ) 

TTFO s 0 

AFFtON PA 1 1 TO MF 
r>A 1 1 TTFP = I TFP 4- 1 
L = o 

no PAP* Ks9,k-LMT 
. HO PUP* J~P * JP l 

IF( IFLAG(K, J, I ) .LT. IFUL .OR. IFLAG<K,J«I> .GT. IFH) GO TO 2421 

PHII^M) - ((PHI(K,J,I + 1) + PHI (K’.J. l-n )/nxny + PHKK.J+UI) 
I + PHI (K. J-l , I) )/PYPY + (PHKK+UJU) + PHI (K— l.J« I ) ) /O^nZ 
o - n( u'. J« 1 ) ) #ror4 - fpc,a*PHI ( J* I ) 
jru .r-Q. 1) on TO P 4?1 

jrMp>» - An c <(PHI (Y ♦ J« I ) ) 

TFMPo = APMTFMP? - TF*AP1 )/( TFMPP + TFMP1) 

[riTFwm ,r.T, FPFP) L = 1 
fONTTMl IP 


I FLAG (K« Jt I ) .GT. IFH) GO TO 2311 


-Vir* 1 


1 

• #MV 


CONTI M IF 

IF(L .no. O) GO TO ?431 

GO TO P 1 F 1 

PALI ri.FOUT ( TOPU ) 

Mftpk'PD PART TOl. FF, 

PF'.tl f MO NUl 
PFwiKtn Min 
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M - mmAPKP 
I f\ ~ 1 

nr. Pi |p ( M .r-r. TPSTA) MbIOFTA 

prAn(NHI) ( X < L ) iLb 1 «M ) , ( Y <L ) ♦Lb 1 *M ) « ( 7 (L ) .Ln 1 ,M) 

no prp 1 Mr- 1 «M 

I R y(M) + p,0 

J n V(M) + P.n 

x n 7(H) + p.n 

Tr*MP i s j 

Tr*»P? = j 

TrMPl a x 


vr 

n 

rrMn j 

- 

1 , F 

vr 

Cl 

TrMpp 

- 

1 . F 

vr 

8 

T^MP? 

- 

1 .F 

no 

PP71 1 a 1 

« 

3 

V! 

8 

TFMP1 

- 

O.F - Y(N) 

Y 1 

8 

xrMPp 

- 

o.F - Y(M) 

71 

8 

tfypt 

- 

O.F - 7<N) 


n = t + i 

ji = j + i 

xi = y + 1 

no To ( Pf-ii ,7^41 .p^fi ) « L 
? c 01 T 1 = ! - 1 

Y 1 = v(M) + p.n - tfmpi 

VO r 1,p — XI 

tr<Y(N) .OT. YO) CO TO P*3*. 

J 1 = J - 1 

Y 1 = Y<N) + P.« - T p MP? 

PF^r YP = 1.0 - Y 1 

1F(7(\|) . GT • 70) OO TO PF.FO 

xi s y - 1 

71 = 7<N> + P.F - TFMPO 

PFpO 7P a 1,0 - 71 

XX a DT# ( Y I #7 1 * ( X 1 *U < K » J , I ) + XP*U ( K , J « I 1 ) ) + Y 1 *7P* ( X 1 *U ( K 1 ♦ J , I ) 

1 + X?*U(K1 . J. t 1) ) + YP*7l*(Xl*U(K. Jl . I ) + X?*UCK. Jl . I 1 ) > 

p + yP* 7P* (XI *U (Kl«Jl«I) + X?*U(K1 . Jl « ! 1 ) ) ) 
on To pr-71 
'•*4 1 Jl = J ~ 1 

Y 1 = Y(M) + P.O - TFMPP 

YP = 1.0 - Y 1 

ir(X(N) .OT. XO) GO TO PF 4F 
11 = 1-1 

XI = X(N) + P.f - tfmPI 

PF^r yp - 1,0 - XI 

I p ( 7( M ) . O.T . 70) GO TO PF 40 

k- 1 a X - 1 

71 = 7 ( M ) 4- O.r - TFMpO 

Or 4 o 7P a 1,0-71 

YY a OT*(71*Xl*(Yl*V(K, J.I > + YP*V ( K , J 1 « I ) ) + 71 *XP* ( Y 1 *V ( K , J * 1 1 ) 

j + y p*V(K,Jl « 1 1 ) ) + 7P*X 1 *( Y 1 *V< K J I ) + YP* V C K 1 « J 1 » I ) , 

p 4 7P*XP# ( Y 1 *V (Xl«J«ll) ♦ YP*V< K1 , Jl , I 1 ) ) ) 

0.0 TO or 71 
pr.r 1 K 1 a y - 1 

71 a 7 ( M ) 4- P.o - TFMPP 

7P = 1.0-71 

!F(V(M) ,0-T. XO) GO TO PFFF 

II = 1-1 

Y 1 a Y(M) + P.F - TF^PI 


B- 12 



LMSC-HREC D225632 


y pa 1 , 0 - y 1 

fFfVfN) • OT • YD CD TO PFFQ 

Jl n J - 1 

VI n Y(M) + P. r > ** TFMPP 

77 n 0T*<X1*Y1*<71*W<K» J. I > + ZP*W< Xl . J. 1 ) ) + XI *YP* ( 7 1 *W < K« Jl « I > 

1 + 7p#W(K 1 « Jl » I ) ) + X?*Y1*(?1*W<K« J. I I) + ZP*W(Kl . J. ! 1 ) ) 

o + y?#YP*(71 #l»l<K« Jl « I 1 ) + 7?* W ( K 1 « J 1 « ! 1 ) ) ) 

*r:7i rONTTMUr 

y (N) s y(M) + xx 

Y(N> = YIN) + YY 

prni 7IN) s 7(M) + 77 

wRITf(NUP) ( X ( L ) «L = 1 « M) « < Yt L ) « L= l ♦ M ) » <7<L).L=1» ) 
m = mMAPKP - 14*TCFT4 
IF ( M .LT. 1 ) on TO proi 
14 = T4 + 1 

no to PRPl 

pro 1 I NOl 110= I N0U 17*1 FI ON 
NIJ1 = Ml IP - INOU1P 
Ml 10 s Ml IP + tNOUlO 

no ?-?i i xsp.ypi 
no P7l i J=P. jpi 
no p -»1 1 ! = P. IP1 

P71 1 n ( K ♦ J . I ) = 0,0 

PFUl T mo Nlll 
M = MMAPyO 
= 1 

TFMPl = I PAP 
TFMPP 8 JFAP 
TFMPp s XOAP 

P7Pt IF(M .OT. ICFT4) M=ICST4 

PFAniNUI ) ( XtL ) . L=1 «M ) . ( Y ( L ) »Ls h w ) * (7(L)«Lsl«. ) 

nO P7P1 Nrl,M 
I s v ( M ) + P*0 

J = v ( M ) + P.0 



y a < m ) + 

P 








p 

!F( 1 

.LT. 

1 

♦ OP# 

I 

.CT. 

TPM 

GO 

TO 

?7 


T F ( J .LT. 

1 

• OP# 

J 

• OT# 

JPP) 

GO 

TO 

27?* 

1 

T F ( y .LT. 

1 

• OP# 

K 

• OT • 

KPp) 

GO 

TO 

27? * 

i* 

!F( I 

,=o. 

1 ) 

XCN) 

s 

0# 1 






IF(J ,ro. 

1 ) 

y <n‘) 

3 

0# 1 





r 

1F(K- .FO. 

1 ) 

7(N> 

e 

O# 1 






I F < I 
IF ( J 
IF(K 
OO Tn 


.FO. 

,=o, 

.FO. 

0771 


IPP) X ( N ) = TFMPl “ 0 . 1 
JPP) Y (N ) sTFMPP-O. 1 
ypp) 7 (N) =TFMPP“ 0. 1 


P7PF PRINT P7 PM IC.YCLF. ITFR. 1* 
07 PA rOPW/lT ( /// • *#** FPPOR 

prri irM 

->771 •’MK.J.T) = 1 .° 

v> 8 MMAPyp - 14 *?cfta 

IF(M .l,T. 1) OO TO P 741 

T4 8 fA + 1 

OO TO F7P1 

774 1 OO Pprr XsP.K’l.MT 
HO Ppre- J=P.JP1 

nO ?*?<= 1 I ap « 1 P 1 
m 8 n ( y « J « T ) 


, K, I 

MR. 


J, X ( N ) « Y(N)« Z < N ) 

7F 1 4 . M 
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TF(N ,ro, o) an to Pin \ 

1F| Ao(K , J# I ) n TF1_AG(K« J# I ) + N 

J1 r r 

07^1 FONT T Ml > r 
nynr FONT ( Ml IF 

KLMT « J1 + P 

l F < ! PPT ( 5 ) ,FO. 1) PRINT P 770# ICYCLF# I TFR 
, n 77 n FORMAT ( /// • I CYCLE e» 17# BY# • J TFR 13 • ! r ' ) 

1 F f MT • 5f . KT* 1 ) MTbK'HJ 
ITMin b 1 PPT ( n ) 

CALI cri.L ID< I FLAT- « PH I # U # V # '" # D # ITFMP) 

A l i |? | fm 1 op 1 TO Mr 
on tp 7?3 R 

•50) t nRIMT p014# NFT M T 

7014 FORMAT!// • FPPOP ** EXECUTION OF Cl IRRFNT CASE IS TFPMJNATFD D 

11 IF TO C TOR A OF LIMIT AT I ON ( • 14# •).*) 

RFTI |om 

fnd 

I r OR, TO 00 7 #00 7 

SUBROUT INF CELL I O ( I FLAG# PH I #U#V#W#D# I TEMP) 

COMMON/L!/ ITYPE# 1 BAP, JBAP# KBAR# IP?# JP? # KP? # I ARRAY ( 50 00 ) # 

1 ARPAYl (500P) ♦ ARRAY?(FPOO) # ARRAY3C5O00) # ARRAY4 ( 5000 ) « 

? ARRAYS! R PPP ) 

COMMON/L 3/ IF# I EH « IS# ISH# IFUL# IFH# IB# ICSTl# ICST3# IP1# 

1 JPl, KP1 # I CYCLE # ICFLL(S)# KLMT 
DIMENSION IFLAG(K-P?# JP?# IP?) # PHI (KP?# JP?# IP?) ♦ U ( KP? , JP? ♦ I P2 > ♦ 

1 V( KP?# JP?# IP?) # W< KP?# JP?# IP?) ♦ 0(KP?# JP?# IP?) 

OAT A tl # I?# 13# 14# IF# tft/1 #?#.7#4#F,6/ 
r. *CK IF Am FMP CLL SHOULD BFCOMF A SUR CLL AND vice VERSA. *ADJ THE 
r \/r-| COmPC OF A NF'i'LV CRFATrD FMP FLL • 

OO 4 11) Ks ? # " LMT 

00 4111 J=p# JPl 

no 4 101 I = ? # I P 1 
m = tfi_a«<k# j# n/tcsTs 
tF(M .FT# 4) 00 TO 4 IP) 

1 . = tflaoik, j# n/tcsTi 

1 3 T H. A C, ( K , J # T ) - I CST 1 *L + 1 
OO To ( 40 ? 1 ,4 r ' c U ) # L 

4 P? 1 F-0 TO (4 IP 1 #4 1 01 #4031 #4041 ) « M 

4031 1 FL AO ( K # J# I ) = 1^ 

OO TO 4PF1 

4041 IFLAr(K#J#1) = IFH 
OO TP 4PF1 

4051 GO TO (4001 #4071 #4061 #4071 ) # M 
4001 IFLAOCF# J# I ) = t 7 

00 TP 4 I 0 1 

4071 TFLAr(F#J.1) = T AH 
OP TP 4101 

4PR1 TCFLt (11)= |Fl..AO(K,J, t + 1 )/10ST7 
10FLI (1P)= t FL. AG (K+l «J# 1 1/JOST7 

1 CFl | (17)= 1 FL AO. ( K # J # 1-1 ) / 1 CST3 

1 CELL ( 14 )= ! = LA0(K-1 , J# t )/tCST7 

J CFt t (15)= 1FL. AG(K# J+1 * I 1/1CST3 
T CFL 1. ( IF)= t FI AO ( K , J— 1 # I ) / I OST3 

J1 = p 

np 4 pp 1 Ms 1 # A 

?F(tPFtl.(N) .O-T. ?) OO TO 4PR1 

OO TP ( 4PM7 , 4PR4 # 4 PR5 #4PRF # 4PR7 , 40PP ) # N 
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rtf fn U(K. J. 1 ) a n,n 

r,o to /jnno 

A PR A UK K • J. t ) n 0.0 

on to a ooo 

acrh i i c k ♦ j. t~n - b.o 

GO TO A OOP 

APR* *iKK-t « J« ! > n o.n 
r.n to aooo 
A 0B7 VCK. J. I ) a n .° 

CO TO A POO 

AORO V(K« J-1 ♦ I ) « n »P 
A OOP J1 a Jt + 1 
A OP 1 OONTtM/F 

t F ( | TFMP .FO. P) fin TO A 101 

I F ( J 1 .ro. O .OR. J1 .FO. 6 ) PRINT 4094, ICyCLF. J1 , 
AC'RA FORMAT <//• ** COMMFNT ** Jl=0. THERE IS A CAVITY. 

11.0 MOVFS TOO FASTI A 101) •• / ♦ VALUES OF I CYCLE. Jl« 

0 C I*) 

A 1 o 1 roMTtMUF 
A 1 1 1 CONTTNUF 

C, #C.K IF A SUR CLL SHOULD BECOME A FUL CLL AND VICE VERSA. 
r Cl L. 

no AOM K’sO.YLMT 
no A-JiM J=2,JP1 
np AoRl 1=2, IP] 
m = T F| AOIF. J. n/IOST** 

TF(M .LT. 7) OO TO ASFl 
1 CELL ( I 1 ) = I FLAG ( K". J ♦ 1 + 1 1/TCST3 
ICFLL (I?)® IFLAGfF+l . J. I )/TGST7 
l CEl-l ( IS) a I FL AO < K « J « 1-1 ) / 1 CST3 
1CFL.I ( I A ) = TFLAG(F-1 , J, I 1/T0ST7 
ICELl ( I 5 ) = I FLAG ( K, J+l » I ) / I CST3 
1 CELL I !6)= I FLAG (K. J-1. I )/ICST3 
no APR 1 L= 1 «F> 

I F ( T ort_L (L ) .LT. 3) GO TO 4261 
A??1 CONTINUE 

r-0 TO (43R1 .47S1 ,4271 ,4241 .A271 ,A?41 ,4751 > , M 
4271 IFLAGIk*. J. t ) = IFUL 

0,0 TO A 3«1 

424 1 IFLAOlto. J. t ) = I FH 

GO TO 4^1 

A 261 GO TO < 43^1 ,4351 ,4271 ,42«1 ,4271 ,4281 .43<=1 ) , M 
A 27 1 t FLAG ( K« J. I ) = IS 
ro to af°1 

428 1 IFLAG(F.J.I) = ISH 
A2°1 M = o 
*1 = 1 

no aoP) L=1 ,6 

IF( im L (t. ) .LT. 7) N=N+K1 

ARO, \S\ - o#l<| 

f FLAG ( is. J. t ) a IFLAG(K.J.I) + N*tCGTl 
APF1 GONTtMtlF 
A "'FI FONT TNI IF 
DFJI ION 
FMO 

•FOP. IS CPA. SPA 

SUBROUTINE VFLCTYI I FLAG .PH I .U.V.W.D) 

C0MMON/L 1 / I TYPE. I BAP. JBAR. KRAR, IP2. JP2 , KP2 « 


I ♦ J. K 

Jl=6. FLOW FIE 
I. J. K ARE • 


*ID OF A SUR 


! ARRAY (SOOO) ♦ 
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\ appavi i^nno) , arrayf( 5000 ) , APRAY3 (noon) , arraya< oooo ) « 

lrH. IS. I SH« 1FUL. IFH. IF. ICST.. ICST3. IPl. 

1 JP1. XP1 « I OYOLF, 10FLL.(ft)« Kl.MT 
C0MMOM/L.4/ DXOOY, DXOD7 « D70DX. D70PY , DY°OXi DVOD? , p , 

niMFMOlON 1FLA<3(KPP,JPP.!PP>. PHI (KP7, JP7, IPS) , U ( l*.PP « JP2 « t PP 
, V (XP7, JP?«!PP) , ««MKPp*jP?*|PP)* D(FPP* JP?« IP?) A r „ 

vcolp THr U’l COMPS OF A SUR CLL AND OF FMP CUUS JUST OUTDIOE OF A FS. 

no AnM KwOffe'LMT 
no 4n6 1 

no 4 of-i |BP«IP1 

V c JF| AO (Kf J< I )/!OfiT3 

| p ( M # | y # 7 • OP • M #OJ# 4) O0 TO 40^1 
l = f TFI-A«<K*J»! ) - M#|CST.T)/ICST1 

IF (I .OT. 47) 0,0 TO 4P71 

IF(1 • OT • 31) 00 TO 4771 

i r 1 1 .rt, 1 « ) 00 TO 4671 


- nZOOX*(U<K» J« I > 
DZODY* ( V ( K * J , I ) - 


- v<K« J-l « I ) ) 
W ( K- 1 » J » I ) ) 


- 1 )) 
V(K, J-l , I ) ) 


IF (I .OT. 47) OO TO 
1F(l .OT. 71) 00 TO 4771 

I (f ( (_ ,rr. 1 « ) 00 TO 4671 

00 = T0 (4531 ,4541 , 4535 ,4545 ,453 1 ,4551 ,4541 ,4561 ,4555,4541,4531, 

1 4565,4561 ,4545,4*31 ) , N . , ... 

4571 U(K, J, I) = U(K,J,I-1> - nX0DY*(V<K,J,I> - V<K«J- J »> 

U ( K , J , I ) = U(K,J,I) ” DXODZ* (W(K,J,I) - W(K-1,J,I>> 

00 TO 4«5l 

4575 U(K,J»t> = IKK, J, 1-1 ) x , ...» 

454 1 w ( k , j , 1 ) = W(K-1,J,1) - nzoox*(u(K, j, 1 ) - u<k,w, n 
W(K,J,T) * W(K,J,1) - DZ0DY*(V(K,J, I) - V(K,J-1,D) 

r.o to 4 P 5 i 

4545 U(K,J,!-1) = U ( K , J , I ) + DX00Y* ( V ( K , J , 1 ) - VtK.J- .1) 

U<K,J,J-1> = U<K,J,t-t) + nx0DZ*(W(K,J, 1 ) - tf(K-l,J,I)> 

OO TO 4051 

4 55 1 I I ( K , J , 1 - 1 ) = ll(K, J, I ) 

00 TO 4541 

4 CF 65 IKK.J.l) = U(K,J,l — l) 

4561 WtK-1 , J,l) = W(K,J,l) + DZ0DX*(U(K,J,l> - 

W(K-1,J,1) = W(K-l,J,n + DZOQY*(V(K,J, I ) - V<K,J-1,!)) 

o-o to 4051 

4565 IKK»J,1-1) = U ( K , J , I ) 

no To 4561 

ar ’ ? ' So* T r Is63S,s631 .S6SI .46SS.S6W.S63S.4SS1 .S6S1.S66S.S6SS.S63S. 

1 4671,4651,4665,46^5,467*), N 


4 55 1 

45 pp 

4561 


« ll(K , J, 1 ) + 
a U(K , J, I-t ) 

= ll(K, J, I ) 

U(K, J, I-l ) 

= W(K, J, 1 ) + 

= w<K-i,J,n 


DX00Y*( V(K, J, 1 ) - V(K,J-1»I)) 

+ DX0DZ*( W(K, J, 1 ) - W<K-1 ,J, I ) ) 


DZ0DX*(IJ(K, j, I ) - U(K,J,1 — 1)) 

+ DZODY* ( V(K, J, I ) “ V ( K , J— 1 , 1 ) ) 


4 67 1 


4671 IKK, J, ! ) = IKK, J, !-l> 

4675 V ( K , J , 1 ) * V ( K , J— 1,1) “ ( 

VCK, J, 1 ) = V(K, J, 1 ) - HY< 

(IO TO 4051 

464 1 V ( K , J » t ) = V(K,J-1,1) 

0,0 To 4 54 1 

4 64 f * tl(K, J, 1 ) - l K K , J , 1 — 1) 

4651 V ( K , J , 1 ) = V ( K , J— 1 , 1 ) 

1F(L .50. 75) OO TO 4661 
!5(t. .50. 75) OO TO 4655 
'•K K , J , 1 ) = W(K-l , J, 1 ) 

15(1 .50. IO) OO TO 4051 

4 65,- tl(K, J, 1-1) = IKK, J, 1 ) 

15(1 .50. 77) OO TO 4051 

!F(l. ,50. 70 .OR. I- .FO. 

A f,61 <*KK-1 , J, l ) = '»K K , J , 1 ) 

r.o to 4051 


464 1 

464 5 
4 65 1 


V(K, J-l , 1 ) 


- OYOOX*( U(K, J, 1 ) - U(K,J»I-1)> 
DY0D7* ( W< K, J, 1 ) - W ( K— 1 , J , I ) ) 


4 655 


4661 


TO 4635 


B- 16 


4 


LMSC-HREC D225632 


4060 V ( K * J ♦ I ) ~ V(K.J-l.I) 

r,o to 4061 

47 1 cloVo (473G.4771 .4741 .4746.4761 .4736.4701 *4741 .4766.4746.4735* 

1 4771 . 476 n, 47 n 6 , 4361 »47 , »6) . N 

/|771 1 K K * . ) « I ) n ll(K t Ji I *■ 1 1 

,. 7 7r, VIK.J-I.f) n VlK.J.I) + PYOnX*<r<K,J,n - U<K, J. l~l ) 1 

V<K.J-t*t> n V(K«J-Jil) + nY0P7#(W<K,J,!) - W<K~l«J*l)> 
oo tp 4 Roi 

A74 1 V1K.J-1.I) n VlK.J.I) 

60 TP 464 1 

4 74 6 U(K.J.l) « U(K.J.I-l) 

IF (l_ .TO. 41) 60 TO 4760 
4 7«il W(K.J.I) = WCK-l .J. 1 ) 

47 ®!r v(K«J- 1 »!) a VlK.J.I) 

TF(l #F0. 41) 60 TO 4661 

IF(L «F0. 40 .OR. L .FO. 46) 60 TO 4661 
IF(L .FO. *»«) 60 TO 4R^1 
4761 IKK. J. 1-1 ) » IKK. J. I ) 

I F ( L .FO. 36) 60 TO 40*1 
|F(L .FO. 44) 60 TO 4661 

GO TO 4 736 

?>1 <30 TO ( 4635.4531 .4541 ♦ 4? = >35.4545 .4635.4661 .4541 .4561 .4555«454 l ♦ 

j 4671 .46*6,4661 .4646.4641 )« N 


4061 CONTI NUF 
40*1 CONTI NUF 
K a 1 

ctim K a K + 1 

J = 1 

610* J = J + 1 

I = 1 

F 11 1 I a T + 1 

M = IFLACKK. J. I ) /I06T3 

IF(M .LT. 3 .OR. M • 6T . 4) 60 TO R 461 

N a 1FLAG1K. J. 1+1 1/IC6T3 

IFCN .LT. 3 .OR. N .6T. 4) 60 TO 5121 

LGC1 a < I FLAG ( K+ 1.J.I1/I CST3 ) * I CST 1 + I FLAG I K+ 1 » J . I + 1) / ! CST3 
IF1LGC1 .FQ. 11 .OR. LGC1.FQ.12 .OR* L6Cl.FQ.2t .OR. LGCl.EQ.22) 
t U<K+1 , J, 1 ) cUtK, J. I ) -DZODX* < »K K , J. 1 + 1 )-WCK. J. I ) ) 

LGC1 a ( IFLAG1K-1 . J. 1 ) / I CST3 ) * I CST 1 + I FLAG ( K- 1 . J . I + 1 ) / I CST3 
IF(L6C1 .FQ.l 1 .OR. LGC1.FQ.1P .OR. LGC1.EQ.21 .OR. LGCl.EQ.22) 
1 U(K-1 . J. I ) bU(K. J. I )+DZODX*< W(K-1 . J ♦ 1 + 1 ) ~W ( K- 1 ♦ J. I ) ) 

LGC1 = ( 1 FLAG (K.J+1 » I )/tCST3)*I CST 1 + 1 FLAG ( K . J+ 1 . 1 + 1 ) / I CST3 
IFILGCl.FQ.ll .OR. LGC1.FQ.12 .OR. LGC1.FQ.21 .OR. LGC1.EQ.22) 
1 U ( K « J+ 1.1) =U ( K « J . I ) -PYODX* 1VIK.J. 1+1) -VlK.J.I)) 

LGC 1 a ( 1 FLAG t K« J- 1 ♦ I )/ 1 CST3 ) * I CST 1 + I FLAG ( K * J— 1 . 1 + 1 ) / I CST 3 
IF(LC,C.l «FO. 1 1 .OR. LGC1.FG.12 .OR. LGCl.EO.21 .OR. LGC1.EQ.22) 
. tl(K. J-l . 1 )*UCK. J, 1 )+DYODX*(V(K. J-l « 1+1 )-V(K. J- l . I ) > 


'121 N = t FL AG (K.J+1* 1 ) / I C6T3 

tr <m «LT. "* .OR. N.GT.^lCOTO^PFl 

L6C1 . <!FLAGCK*J.!-1)/!CST3)*1CST1 + IFL AG(K . J+ 1 . I - 1 ) / I CST.7 

IF1LGC1 »FQ. 1 1 .OR. LGC1.FQ.12 .OR* LGCl.EC.21 *0R. LGC1.EQ.22 
i V(K»J« I“1 ) a V ( K » J » 1 )+DXOOY*<U( K.J+1 .1-1)“L(K.J«!-1)) 

LGC 1 = ( I FLAG CK.J.I + 1 )/!CST3)*t CST 1 + I FLAG! K . J+l . 1 + 1 ) / I CST7 
I F ( LGC 1 • EQ • 1 1 .OR. LGC1.EQ.12 .OR. LGC1.EQ.21 .OP. LGC1.EQ.22) 
1 v(K. J. 1 + 1 ) *V<K. J. 1 )-DXOOY* UK K.J+1 . I ) ~U< K , J . I ) ) 

LGC 1 a ( | FLAG ( K+ 1 . J « I ) / 1 CST3 ) * 1 CST 1 + I FL AG ( K+ 1 , J+ 1 . I ) / 1 CST3 
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i 


1 


IF (UGC I 

UGC J r 
I F ( UGC 1 


t 


nf>r>1 M n fF| 

I F ( N .1 

U.GC.1 b 
IF (UGC 1 

1 

UGC 1 a 
I F ( UGC 1 

1 

UGC 1 o 
I F ( UGC 1 

1 

UGC I e 
IF (UGC1 

1 


•FQ.1I • OR • UGC1.F0.13 »0R. UGCt.EQ.31 .OR. UGC1.EQ.33> 
V/(K + 1 * J * I ) «V ( K . J • I >~nZOnY*<W<K. J+l .11 "W< K . J ♦ I ) ) 

( I FUAG < K — 1 «J»l )/ICST 3 )*fCRTl + I FU AG ( K- 1 . J+l . I> / I CST 3 
.FQ.ll » 0 R. UGC 1 .FQ .13 .OR. UGC 1 .EQ .31 .OR. UGCJ.EQ. 33 > 
V(K"l . j. I> «VCK* J. I > 4 -D 70 DY>MW(K-l » «J<M « I ) — W C K— 1 . I ) > 

AG(X+1 . J. I >/|OGT3 

T. 11 .OR. N #GT. A) GO TO 5461 

( |FUAG<K,J,I-1 )/!CSTG>*lC6Tl + 1 FUAG ( K+ 1 . J » I “ t> / l CGT3^ 

.FQ.ll .OR* UGC 1 .FQ .13 .OR. UGC 1 .FQ .31 .OR. UGCl.EO.r-A> 
W(K, J, 1-1 > nW(K. J. I M-r>XOnZ*<U<K+l ♦ J * l”' )~U(K» J. I ~1 > ) 
( IFl AG( K . J. I + 1 )/ 1 CST 3 ) * I OST 1 + I FU AG< K+ 1 . J ♦ I + 1 > 3 I CST 3 
•FQ.II .OB. LOCI. EQ. IS • OR. LSCI.CO.SI .OH. LQC 1 .E 0 .SS) 
VI ( K . J . I + 1 > eW ( K « J ♦ I > -nxonz* ( U ( K + 1 . J . 1 > -u ( K « J . I > > 

( l FUAG < K » J~ 1 . I > / 1 CST3 > * I CRT 1 + !FLAG(K+ } , J-l . I > /I CRTS 
•EO.ll .OR. UGC1.E0.12 .OR. UCC1.FQ.21 .OR. UGCl.EQ.rA> 
W ( K , J- 1 , I > °W( K » J » l >+DYODZ* ( V(K+1 » J** 1 * I ) “V < K. « J- 1 » ! ) > 
( IFLAG(K.J+1 .1 >/ICSTR)*!CSTl + IFUAG(K+ 1 ♦ J+l . I ) /ICST3 
•FQ.ll .OR. UGC1.F0.1? .OR* UGC1.EQ.31 .OR. UGCl.EO.32> 
W(K, J + l . I )aW(K, J, I >-DY0r>7*<V<K-H , J ♦ I ) -V ( K , J » l> > 


545 1 CONTINUE 

I F ( I .UT. IPt ) GO TO Sill 
1 F ( J .UT. JPt ) GO TO S10G 
I F ( K .UT. KUMT > GO TO *101 
RFTURN 
PNG 


7(3000) 
CORPS 1 i 


• FOR , 1 r SOP, SOP 

GUBROUTINF PUOT (tl,t?»t3»14,M) 

C0MM0N/U3/ ! X ( 1 00 ) « 17(100), X(?000), Y(?OPO), 

1 PUTY, RPPUU, MARGNX, MARGN7 , NSEGPT, NIMRKR, 

? C0MTR1, C0NTR3, C0NTR3, C0NTR4. XV (3, 31), ZV(3,31) 
GO Tn ,ST31 ) , II 

pp?1 IF (I A .GT, 1) GO TO 654 1 
00 SP31 1 = 1,0 


, PUTX, 
SlNPSl 


J a T + 1 

PP31 CAUU L TNFV( IX( I ) , I7( t ) , 1X( J) . I7( J) > 
no BP3 r ' 1 = 1 1 ,NSFGPT,2 
J a T + 1 

p,c,1K CAUU L I NE V ( I X < I),I2( 1 ) , I X ( J ) » I Z < «J ) > 
MCOUKIT = O 


PP/| 1 CONTf Ml IF 

00 p-?P 1 I = 1 , Y 
MCOIImT = MCOUNT + 1 
VP 3 x ( I ) 

YP a Y ( I ) 

7P a 7(1) 

I F ( MrOI INT .GT. NIMRKR) GO TO 5601 
GO Tft ( PPfil , ppo1 ,gpoi ,cpo 1 ) , tn 
5561 IF(YP .UT. PUTY .OR* XP .GT* PUTX) GO TO 6601 
GO TO «7P1 
P'-RJ CONT tMUr 

GO TG 5701 

5^51 t F ( t p . FO • 1) GO TO r -6Pl 

GO TG ( R 61 1 , R 6A 1 ,564 1 , c 64 1 ) , t 0 
5611 1 F ( XP .UT. C0NTR4 .OR. XP .GT. CONTRl) GO 

IF(YP .UT. CONTRA .OR. YP .GT. C0NTR2) GC 

IF(7P .UT. CONTRA .OR. 7P .GT. C0NTR3) GO 


TO 6701 
TO 67C t 
TO 5701 


OG TG PGP1 
P64 1 CGMTlNItr 
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BOR) jy n PPPlJ|_*<VP + YP*COSPSl) 

J7 n PPPUL*(7P + YP* r I NPS I ) 

jv n jy + maponx 

J7 P J 7 + MAPGN7 
CALL PI_OTV(JX« J7*RH,r ) 
n7Ri roMTfMiir 
prn ipm 

n7?1 rOMTlMljr 

no n«f<l 1 I n 1 , J f\ 

.JV r V\M !*■»»! ) 
p v V ( T o ♦ ! ) 

KV c XV ( ! 7 « I + 1 ) 

Y7 n 7V ( l 7 « l + 1 ) 
r7T) CALI | IMrV( JX « J7«FV,y7) 
no o*?a 1 1 » 1 » m 
jy e y ( I ) 

J7 B 7(1) 

KX B y ( JP+I ) 

k'7 b 7(17+1) 

TALI L!NPV( Jy« J7«KX*K7> 
yp a k'X - jy 
7P S K7 - J7 
J a PODT(VP*»? + 7P##P) 

B741 I F ( J • GT • 7) CALL ARROW ( JX « J?. ♦ KX «K? « 6 ♦ 2 ) 

ROO| CONT T K H l r 
RFTUDN 
PMO 

iFOP, ! S SP) »rp) 

SlIRRnilTINF CLKOUT ( TCPU) 

CALL crLCCK (OATF « T l MF »FSFC ♦rftOScx ) 
wp i tp ( *>« 1 ono ) TIxf 
(OOn POPMaT( AHOTIWPb A) 7) 

CAL! CPHTIM(tTlV) 

F RP C a FLOAT ( T T 1 M)/t *F6 
WR!TP(ft«7000) FRPC 

; pn OO FORMAT( IB+rFRBC (CPU) = F 1 4 • 4 ) 

1F(Fcpc .LTi TCPU) RETURN ! 

PRIMr mrn 

ii 7000 FORM AT ( * ** MFRSAGF ** FXFCUTION IS TFRMINATFD TO PROTFCT SC-4070 

i! ) PL.OTC prFOPP MAX TIMF IS PFACHFD# •) 

STOP 

,* c-Nn ! 

•J t/\cM,T( Cpt ITTM,CPUT!M • USP AS CALL CPUTIM(JTIM) 

T, ( 1) AXRS . WHFRF I TIM IS FLAPSFD CPU T I MF j 

rDUTIM# | A AO « ( ARRAY ) • IN MICROSFCONOS 

1 c-p POTS 

j A AO, APPAV+P? 

yc | , yti Ao,pon 

T r a AO,*0,Xt 1 

j p»yti 

APPAY HFR 77 

rMO 

«FOP,lc e?<s,spp 

cUPRmiTlMP ARROW ( I y 1 « tVl«lXP*IVP«IHl TF ♦ I BASF ) 

Hr f O-mT = I M t T r 
' nAsr = IBAsr 

yR> a 1X7 

YP a tVP 
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r»Y n JYP - IY1 
nx n lyp - ixj 
rf) n t • /WT ( r>X*nX + nY#DY ) 

r act v n nA r ,r * (cn«ov) 
patty n r'A'if’ k (prj*nY) 
y? n yo - HP|THT * ( r n*OX) 
yi yn „ Hr T T.HT * < r .Q» n Y) 

J V/( - (VI 4- FATTY) + » r 

JYA r ( V r » - PATTY) + « r > 

(yr, „ (VI - PATTY) + « n 

|V<~1 r> (YT + PATTY) + * R 

TAIL. L t M r V ( |YA« f YA « IX?. IYP) 
TAII |. I MTV < IX^* !Y r >« !XP» IYP) 
rjpTUnN 

pmh 

»map,u. ppnnyppor, 

(_jn <-y<*«^perrT» 
t H e T> 1 

*rOPO()T top*. ♦LWWAT'A. 
inriiiiMn, 1 ( M*'Ar>, 

• pope- I.HMArn* 

• vot ppnr. 

***# oata pppp #*■** 

tFJM 
t F I M 
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